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CHAPTER 1 
THE MOTIVATION FOR ANALOG COMPUTERS 



Analog Computation, based on the modern electronic 
analog computer, is of fairly recent date. The first 
commercially available general purpose electronic 
analog computers appeared on the market in the 1940's. 
These early machines were an outgrowth of an emerging 
electronics technology and a critical need for automatic 
computing machines that could solve complex dynamic 
problems. Slide Rule and manual equation solving could 
no longer be relied upon as a practical approach to 
seeking engineering solutions. Therefore, analog com- 
puters became important tools in the design of air- 
craft, jet engines, atomic reactors, oil refineries, 
chemical plants, etc. 

Many types of analog computers have evolved over the 
years. The family has included the mechanical differ- 
ential analyzer, electromechanical differential ana- 
lyzer, and most recently, the iterative differential 
analyzer. 

The analog computer has always had several advantages 
compared to a digital computer. These are primarily 
speed, more simulation capability per dollar, an ability 
to integrate, and an excellent man- machine interface. 



The main feature of the analog computer is that it can- 
integrate time- varying voltages. There is no easy 
way to differentiate. Consequently a mathematical 
model of aphysical system which is expressed in terms 
of differential equations cannot be solved with the 
machine directly. It is necessary to reformulate the 
mathematical model in terms of integral equations, 
either implicitly or explicitly. 

The analog computer can integrate only with respect 
to time. Thus, a mathematical model which contains 
partial integrals (corresponding to partial derivatives) 
with respect to several variables must be approximated 
by a set of ordinary integral equations with respect to 
time. Of course, computer time need not correspond 
to time in the physical world, although it usually does. 

Before 1959 the analog computer was a synchronous 
machine; all its integrators operated in unison. In 
1959 the DYSTAC^was introduced. The name was an 
acronym for 'dynamic storage analog computer'. This 
machine was the forerunner of the iterative differential 



DYSTAC is a registered trademark of CSI. 




Two large SD 80 computers are used -in this engineering laboratory to simulate the behavior of 
drone helicopters with different load configurations and various automatic flight control sys- 
tems under a wide range of operating conditions . Simulation saved time and money ^ reduced the 
exposure of personnel and materiel to possible damage when actual flight tests were mad^j and 
it also eliminated the necessity of waiting for special environmental (weather) conditions . 

(Photo courtesy, Gyrodyne Company of America, Inc.) 
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analyzer which appeared in 1960. The iterative differ- 
ential analyzer is an asynchronous computer; the 
integrators need not be controlled in unison. They 
can operate independently either in groups or singly. 

In 1962 the analog computer was augmented with d^ital 
logic. This innovation first appeared in a machine 
called the HYDAC.^ The HYDAC had a very large 
quantity of synchronous digital logic. Since that time 
the inclusion of a smaller complement of asynchronous 
d^ital logic has become accepted practice. 

Digital logic can be used for the implementation of 
logical decisions. These are based on results obtained 
from the analog portion of the computer during the 
solution of the problem. Digital logic also can be used 
for mode control of the analog computer. All of the 
Systron-Donner 10/20 and 40/80 series analog com- 
puters can be operated in the iterative mode and can 
be augmented with digital logic. 

Analog computers have found widespread acceptance 
in virtually every area of scientific investigation. This 
growing interest in analog computers has created a 
need for complete software, specially designed for the 
beginner and less experienced user. Also, the recent 
addition of digital logic control has greatly improved 
the problem- solving capability of analog computers. 
How this new feature can be used in analog computation 
is thoroughly illustrated in this publication. 

It is the purpose of this handbook to provide students 
as well as experienced computer users with compre- 
hensive and up - to - date analog computer software. 
Chapters 2 to 5 develop the basic fundamentals of com- 
puter operation and illustrate the solution of elementary 
problems. Chapter 6 provides a detailed description 
of the operating controls and computer logic of the 
SD 10/20 and 40/80 series computers. This informa- 
tion serves as useful reference material to problem- 
solutions illustrated in the more advanced discussions 
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which relate iterative programming techniques to the 
Systron-Donner computers. Chapters 7 to 22 develop 
more sophisticated programming techniques and appli- 
cations on a progressive basis. 
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On the Apollo Program^ Douglas Aircraft 
coupled an SD 40 computer to a large centri- 
fuge. The computer calculated and integrated 
errors in human performance during Apollo 
lifting body reentry simulation studies . The 
Systron-Donner computer is seen next to the 
programmer y in upper left portion of picture . 

(Photo courtesy, Douglas Aircraft Co. ) 





Analog computers are now standard computing equipment for classroom teaching and research work 
in colleges and universities . Students in the departments of Electrical Engineering , Mechan- 
ical Engineering , Chemistry and Biosciences receive instruction in the use of desk top analog 
computers as basic electronic model builders of dynamic problems . The ease and swiftness of 
presenting a solution on a readout (oscilloscope, XY recorder) , and the ability to vary problem 
parameters and observe immediately their corresponding effects , have made the analog computer 
an important teaching aid. 

(Photo courtesy, University of Santa Clara) 
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CHAPTER 2 
BASIC ANALOG COMPUTING ELEMENTS 



This chapter shows how electronic equipment and 
circuits are used to implement mathematical rela- 
tionships in an analog computer. 

OHM'S LAW 

Ohm's law describes the relationship between the 
current through, and the voltage across a passive 
impedance. A passive impedance is a collection of 
passive elements such as resistors and capacitors 
connected together in an arbitrary way. Such an impe- 
dance is generally considered (in analog computing) 
to be a two- terminal network which can be denoted by 




For a capacitor 



so that Ohm's law is 



or 



where 



C ^= I 
^ dt ^ 



^ - -s' 



E = ZI 



2 = ci 



where Z is the dynamic impedance of the element. 
Ohm's law states 



for 



E = ZI 



Hh 



-H E - 



SERIES and PARALLEL IMPEDANCE 
Impedances are additive in series: 



where I is the time- varying current (in the direction 
indicated) through the passive element generated by 
the impressed time-varying voltage, E (with the 
polarity indicated). 

In order to simplify what follows, transform notation 
will be used where 



+ £■ 



+ E 



Zt = Zi+Z2 



and 



sf(t) =-^f(t) 



-^f = fmat. 



This can be proved in a simple way. The total voltage 
from point a to point b is E-p = Ej + E2 and by 
Ohm's law 



E^ = Z^ I + Zg I = (Zj + Z^) I 



Z^L 



For a resistor Z = R and Ohm's law is 
E = RI 



When impedances are connected in parallel the total 
impedance can be found as the reciprocal of the sum 
of the reciprocals: 
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Zt Z| z^ Z3 



This again can be proved by simple application of 
Ohm's law: 



^1 ^2 ^3 



^T ' E - Z^ Z2 Z3 



The total impedance of two parallel impedances has 
the simple formula: 



1 



V2_ 

z^.z^ 



Typically, Ri ranges between 1 ohm and 100 ohms; 
and R2 between 100, 000 ohms (lOOKfi) and 1, 000, 000 
ohms (IMn). Thus the diode is an approximation to a 
switch for which 

R. = 



R2 - « 
OPERATIONAL AMPLIFIER 

An operational amplifier is of the type that is called a 
d-c amplifier; it amplifies not only time-varying volt- 
ages, but also d-c or constant voltages. It is charac- 
terized by its excellent stability and extremely high 
low-frequency gain (amplification). It is denoted by 




where 



= -A e. 



These voltages are measured with respect to ground 
(zero reference). The gain, A, is usually frequency 
dependent and will decrease with increasing frequency. 
As will be seen presently, this becomes a limiting 
factor in the use of a computer at high speeds. An 
equivalent circuit for the amplifier is 



e: o- 




DIODE 



A diode is a non-linear resistor. Its resistance or 
impedance depends on the direction of the current flow- 
ing through it. It is denoted by 



I 



where 



and generally 



— ^ 

+ E- 

E - ZI 

Z = Rj, E > 



= Rg, E < 



^2 » ^1 



where, generally, Rj is greater than 10 ohms and 
Ro less than 10"^ ohms (closed- loop) at zero fre- 
quency. However, the amplifier is current limited. 
That is, it will only perform satisfactorily if the out- 
put current, Iq, is less than some value. The cur- 
rent limit for the S-D^ amplifier is ±25 milliamperes 
,at ±100 volt output. The amplifier is also voltage 
limited; it will not function satisfactorily if the output 
voltage, eo, is greater, in absolute value, than some 
upper limit. The limit for the S- D solid- state 100 volt 
amplifier is 105 volts. Since the output of these amp- 
lifiers contains unwanted noise whose magnitude typi- 
cally can be 10"2 volts, their effective useful range is 
about three and one-half decades (5xlO~2 to 102 volts). 

GENERATION OF TRANSFER FUNCTIONS 

In what follows, it will be assumed that the input 
impedance of the operational amplifier is infinite, the 
output impedance is zero, and the gain is infiniteii. 



Abbreviation for Systron-Donner. 
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These assumptions introduce negligible error at zero 
frequency and are a good approximation at mid- fre- 
quencies. At high frequency, the assumptions cannot 
be made (particularly for gain). 

Consider an amplifier with input and feedback (from 
output to input) impedances and applied input voltages 
as shown. 



and 
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The simplest case is 
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From Ohm's law 



which can be expressed 



(e. -e.) (e„ -e ) (e„ -e ) 

I = -J: L + _e L + -A L 

T Z, Z„ Z„ 



e. -e 

I = -J 9. 

h Z^ 



(3) 



By a suitable choice of impedances many desired trans- 
fer functions can be generated. (Henceforth, the units 
of megohm and microfarad will be used for resistance 
and capacitance respectively. ) 



•t = If * h 



Now if the input impedance, R., of the amplifier is 
assumed to be infinite 



SUMMER (ADDER) 

Consider the configuration shown below, which is 
called a summer. 



1 



for finite e. . Thus 
] 

^T = h 
and 

(e. -e.) (e„ -e.) (e„ -e.) 

_i L + _^ L + —^ L 

^1 ^2 ^3 



e. -e 



Further, if the gain, A, of the amplifier is assumed 
to be infinite, then 



e9-~AAA^ if 



O^ 



If, in equation (2), the substitutions 



Zj = Rj 



. = ^ . 

This is a reasonable assumption at low frequencies 
since as noted above eQ, at most, will be in the neigh- 
borhood of 100 volts in absolute value. Thus, equation 
(1) becomes 



^2 ~ ^2 



^3 " ^3 



^f - ^f 



are made, then 



1l + ^ + ^ = - 1°. 
^1 ^2 ^3 ^f 



R 



R7 ^1 



^e 
Rg ^2 



ie 
R3 ^3- 
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A set of typical values (in megohms) for these resistors 
in a summer in a computer is 

R^ = R2 = Rj = 1.0 



R„ = 0. 1 



so that 



-e, -eo -lOe, 



In a computer, access is usually provided to the input 
(summing) junction so that additional input resistors 
or feedback components can be added to the summer 
circuit externally. Thus, the summer circuit is 



1.0 

-^AAr 



62- 



1.0 



10 



0.1 



'^T^ 




Substitution in equation (3) of 



Z - R, Zj - ^^ 



leads to the transfer function 
o _ 1 



RCs 



or 



®o " RCs "i 



e. = - 



kf^ 



RC / ^i<'>^' 



K R=l, C=l, then 



/•■ 



(t)dt 



so that the output voltage, e , of the amplifier will be 
the integral with respect toT:ime of the time-varying 
input voltage, e.. 

If a switch is included as shown 



where the gains (multiplying factors) are indicated at 
each input, and the summing junction terminal by J. 
The program symbol for the summer is 




10 




and if 



1) there is an initial charge stored on the capacitor 
which results in the voltage, e (0) 



where the input gains are omitted if they are unity. If 
it is desired to indicate the junction or summing junc- 
tion (high- gain input), the symbol becomes 



10 




INTEGRATOR 

The high-gain input is labeled only if pertinent. Con- 
sider the configuration shown below, which is called 
an integrator 



e; ^AAr 




eo 



2) the switch closes at t=0, 

then 

t 

o 

It is common practice to call this initial voltage the 
'initial condition' for the integrator. The derivation 
of this terminology is obvious from mathematics. 

It is necessary to find a practical way to establish the 
initial condition or initial voltage for the integrator. 
To accomplish this, the circuit below is used. 



-e 



IC 



— ^AAr-^ 
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Here 



With S. closed at, say, t=0, then 



Cs . ^ = RCs . 1 



R 



R 



RCs + 1 
= R 



and substituting in equation (3) 
1 



^0 = <■ 



RCs + 1^ ^~hc^ RCs + 1 ^IC 



t 

%(t) = ^ic - /^Wdt. 
o 



Katt = T, Si is again opened, Bq will stop chan- 
ging and remain at the last value before S- 
opened, namely 



e„(T) 



IC 



/^' 



(t)dt. 



RC 



de 
c 

dt 



+ e = e,„. 
o IC 



The solution of this differential equation is 



% = ^^^ fe) + ej^' 



where A is a constant depending on the initial voltage 
stored on the capacitor before e,„ was applied. Then 



t — CD O 



IC 



For practical cases, it is only necessary that t>10RC, 
since e~10<, 0001, which is compatible with the accu- 
racy of the circuit. The result in simple terms is: 
to guarantee the establishment of the initial condition 
it is necessary to wait lORC seconds after the voltage, 
-e,„, has been applied to the circuit. 

The two previous circuits, combined with appropriate 
switches as shown below, constitute a practical inte- 
grator 



K S2 is closed again while Sj is open, the out- 
put of the integrator will return to 



= e 



IC 



Each integrator in a computer has these switches as 
part of its circuit. They are open or closed depending 
on what the programmer wants the integrator to do. 
The state of these switches is called the integrator 
mode. If all the integrators are controlled in unison, 
the switch states are determined by the main computer 
mode. 

The modes have simple names with obvious interpre- 
tations. During the R, or reset (IC) mode, initial 
voltages are impressed on the integrator capacitors. 
During the compute or C mode, the integrators inte- 
grate input voltages. During the hold or H mode, the 
integrator outputs remain constant at the last value 
achieved before entering the hold mode. Thus, the 
mode permits the programmer to stop the computation 
at any time, enabling him to evaluate what has happened 
thus far in the calculation. The table below shows the 
states of Si andS2 in the various modes. The numeral 1 
indicates the switch is closed (logical 1) while 
indicates the switch is open (logical 0). 




The operation is as follows: 

With S2 closed and Si open, eo = ejc after 
1.0 second (or 0.1 second if C = 0.1, etc.). 



With S2 open, eo still equals ejQ, since the 
current through the capacitor is zero. (This 
must be true since both Si andS2 are open and the 
input impedance of the amplifier is assumed to be 
infinite. ) 



MODE 



R(IC) 



H 



CIO 

Finally, the general integrator circuit is shown below 
together with a definition of its transfer function for 
various modes. 



=IC 



IC 0.05 
o ^AA^ 
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e|0— ^\AA^-i 
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MODE 
R 



TRANSFER FUNCTION 



% = he 



%^^^ = ^IC 



rle^it) + e2(t) + 10e3(t)ldt 



eio 




H 



%(T) = ^IC 



/[• 



ej(t) 



e2(t) 



10 63(1)] 



dt 



Note that the gain (multiplying factor) of the computing 
inputs is noted next to the terminals and is determined 
by the value of the input resistor. The SJ terminal is 
provided to allow the other external input resistors to 
be added. The U terminal provides the ability to gen- 
erate an initial condition which is the sum of several 
voltages. The J terminal permits the connection of 
additional feedback elements around the amplifier (to 
be discussed in chapter 10). Also, these terminals 
provide an external connection for the use of solid- 
state switches in place of the mechanical switches 
(relays) Sj and S2. In this case Sj and Sg are con- 
strained to be open regardless of integrator mode and 
the necessary switching is done by the external solid- 
state switches. 

The program symbol for the integrator is 




The input gains are usually omitted if they are 
unity. 

POTENTIOMETERS 

It is necessary to have a device for entering constant 
parameters in the computer program. This is accomp- 
lished with a potentiometer. The circuit for this ele- 
ment is shown below. 




Then 



so that 



'2 ■= ll * II 



«2 ' «1 ' «L 



and 



Since 



R2 = % 



R, + 



Rl Rl " R2 



1 ' ^ - ^t 

where R, = total potentiometer resistance, 



and 





Rg - R^ - Rj 






1 




^t ^t Ri 


e. 
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e 



«1 \ ] 


^t ^^L "" ^1^ " ^1^ 




^ -^ 


e = 



i e.. 

1 





Now, obviously, it would be time-consuming to deter- 
mine Rj, with the knowle(%e of Rt, Rl for each new a. 
Consequently, in practice, the pot is set with the load 
connected, by reading the output voltage, e , with a 
meter, foraknown input voltage, which is us*&ally 100 
volts. 



Thus the meter reading is equal to 100 a and the pot is 
changed until the desired value for a is obtained. The 
symbol for the potentiometer is 



If the output voltage, Bq, of the potentiometer (pot) is 
applied to an input of another element of the computer, 
then the input impedance, Rl, of the other element is 
connected from the pot output to ground. 




or 
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where H is the notation for the high end, or input. The 
His usually omitted since it is obvious from the com- 
puter program which side is the input to the pot. The 
element discussed above is called a two- terminal pot 
since it has two available terminals (or connections) 
on the computer program board (patchboard). Some- 
times it is desirable to connect the bottom (low) end 
of the pot to some computing- element instead of ground. 
In this case the low end is made available at the patch- 
board. The circuit is 




and the program symbol is 



The principle of operation of the function generator 
depends on the use of a non-linear input impedance 
for an operational amplifier. That is, the impedance 
generates an input current proportional to the fxmction 
which in turn constrains the output voltage to have this 
functional relationship to the input voltage. A simpli- 
fied circuit is shown below. 



1.0 



.lK3i 




-^1 = ^F 



\ = - f (ej) 



e = f'(ei) 
o 1 




or 




where notation for the arm. A, is omitted when it is 
obvious from the computer program. 

ARBITRARY FUNCTION GENERATOR 



The program symbol for an arbitrary function gener- 
ator is 



FIXED FUNCTION GENERATOR 



Many mathematical problems to be solved with a com- 
puter require the generation of an arbitrary function. 
This is accomplished with a device called an arbitrary 
function generator (or frequently, diode function gen- 
erator, since the internal circuitry uses diodes). This 
device allows the programmer to approximate the de- 
sired function with stra^ht line segments. An example 
is shown below 



Fixed function generators are used to generate often- 
used functions such as sine, cosine, log. etc. These 
are similar to arbitrary fimction generators in opera- 
tion, but do not permit the programmer to change the 
parameters within the device. Generally, a fixed 
function generator is more accurate and has better 
frequency and noise specifications than an arbitrary 
function generator. The program symbol for a sine 
generator, for example, is 



f'(e,) 




Here the function f(ej)is approximated by f'{e^) with 
three line segments. In general, each function gen- 
erator, depending on how it is used, will provide either 
10 or 11 line- segments with maximum slope changes 
of 2 or 2. 5:1. Several of these devices can be used to- 
gether if more segments are required. The location 
of the slope discontinuity is called the breakpoint and 
is adjustable. Detailed instructions for the setup of 
this element appear in Chapter 6. 



x(t) 



Sin 



> 



Sin X ( t ) 



MULTIPLIER 

The multiplier used with the S-D computer is called a 
quarter- square multiplier. The name derives from the 
equation 



XY =^ [(X + Y)2 - (X-Y)2] 



The actual equation to be implemented with hardware 
in order to provide multiplication is 



®1 ®2 " 400 



^1 ^ ^2 



(4) 



The multiplier module contains two fixed function gen- 
erators. A block diagram for one of the generators is 
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The transfer fxinction is 



-^I 



(Note that both polarities are required as inputs for both 
e^ and e^. ) The transfer function is 



I = K 



^1 ^ ^2 
400 



I 



K 
400 



1^ ^ ^2l 



J - ^ IC 



I "1 "21 J 100 

A simplified block diagram for the module is 



■e,«- 



■eo»- 



where K is a constant that determines the feedback 
resistor of the output amplifier (to be shown below). 

The block diagram for the other is 



e, •- 



62* 



-ei»- 



62* 



The transfer fimction is 



If the module is connected to an output amplifier 



then 
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The multiplier module block diagram is 



and if K = r^ 
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The program symbol for the multiplier module (with- 
out the output amplifier) is 





\ 

X 

/ 


Z 


/ 

I 

\ 




^1 








^1 






=2 




-e 


2 



With the designations 

M: The module is used as a multiplier 
I: The output is a current 
X, Y: Bipolar inputs 
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If the bipolar inputs are not naturally available from the 
program (e.g. other variables) then summers can be 
used to generate them. This can be done in two ways: 

Case I: 



ex 



ey 







ei 


\ 

X 

/ 


t 


/ 

I 

\ 
















>^ 












62 










e„ e 






f\-e2 


1 = 


^x ^ 
100 






l^ 









CIRCUIT 



Note that the "minus sign" associated with the "M" indi- 
cates a negative current is generated for a positive 
product. 

The program symbol for the multiplier module together 
with an output amplifier is 




eo= - 



^X ^y 
100 



As before, summers can be used to generate the bi- 
polar inputs if they are not otherwise available. 



Case I: 




PROGRAM SYMBOL 



Note that the "plus sign" associated with the "M" indic- 
ates a positive current is generated for a positive pro- 
duct, CxCy. The sign is generally omitted since this 
is the normal condition. 




CIRCUIT 



Case n: 







ei 


\ 

X 

/ 


:>^ 


/ 

I 

\ 














Nsea 




fix * 












-e 

I 


2 

-ex e 


ey ' 




[> 






L— . 




100 



CIRCUIT 




° '00 



PROGRAM SYMBOL 



(Note that here the sign inversion is due to the output 
amplifier.) 

Case II: 





PROGRAM SYMBOL 



CIRCUIT 
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PROGRAM SYMBOL 



I \AAr 



''% 



{> 



■eo 



A multiplier module can be used to make a divider with 
the same output amplifier. (This comiection can be 
made conveniently at the patchboard. ) 



exO~— VVA> f 




Here 



e e e„ 

T - K-9-L - - ^ 
^ ~ '^ 100 R 



Rf<0. Then, I f will be in the direction shown, for 
eo > 0, and If will produce an e j < . Since the gain 
of the amplifier is negative and large, the circuit will 
be unstable. An additional restriction, |exl - l^yl ' ^® 
necessary to ensure that the output amplifier not ex- 
ceed the voltage limit. 

The divider can be used to generate square root in a 
simple way. If e = e^. 



e— vvA/' 




so that for R = jf 



CIRCUIT 



Vy 

^x " " 100 



then 



or 



lOOe 



e = - 

o e. 



or 



e = -100-2- 



The programmer's symbol is 



e = -100 e 
o X 







t 


\ 


ex 


• 


/ 






. 


100 ex 


Cy 




cy 



•eo 



From the circuit it can be seen that the condition ey>0 
must hold in order to avoid instability. This is easy 
to show. The multiplier module can be considered to 
be a non- linear feedback resister, R^. The sign of 
Rf is the same as the polarity of e . Suppose in the 
circuit 



so that 



e = 10 Ae 
o V X 



v°- 



The diode insures that the system will not saturate in 
the wrong direction if e^ goes positive inadvertently. 

This circuit has the program symbol 



ex 



> 



•eo IO>/e7 
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FUNCTION SWITCH 



COMPARATOR 



A function switch is a manual switch which can be oper- 
ated from the control console. It is used to change the 
computer program during execution. The switch has 
three positions: up, down, and center (off). Its circuit 
is shown below. 



Frequently, it is required to determine the sign of the 
sum of two variables. 

A circuit which accomplishes this is 



I o 



D2 



This is the program symbol as well. The terminals 
of the switch are located at the patchboard. 

FUNCTION RELAY 

Afimction relay is usedto make program changes auto- 
matically durii^ execution. It is energized (logical 1) 
for an input equal to or greater than +28 volts and is 
de-energized (logical 0) for an input equal to or less 
than volts. The circuit is shown below 




Ux 




U 






Ux = X if U = 




Ux = X if U = 1 



This is also the program symbol where N is the relay 
number. 

ELECTRONIC SWITCH 



«2 .068 



0.1 
0.1 



ej 




.02 
— o Co 



-100 



D| 



If (ei + e2) > 0, then eQ is limited to approximately 
volts by Di because of its low resistance when it con- 
ducts (logical 1). If (ei + e2) < 0, then eo is limited 
to approximately +28 volts by Dg (logical 0). The 
latter is true since 



1) Ip mustbein the direction shown (e. + e2>0) 
and e^ > e... 

2) eo = 28 volts since ez must be approximately 
volts due to the high gain of the amplifier. 

The program symbol for the comparator is 



=[!> 



U = l if (a+b)<0 
U = if (a+b)>0 



An electronic switch is used in much the same way as 

afunction relay; the only restriction being that it must 

beinseries with the junction of an amplifier. Its speed 

of operation is much greater (10~5 sec) than afunction 

relay (10-3 sec). It requires the same input voltages 

(logic levels) as the function relay (i.e. logical 1: 

28<e. <100, logical 0: -100<e. <0). 
m ' ^ m ' 

The program symbol is 



If both logical comparator outputs are required, then 
Cq is complemented with a biased analog inverter. The 
program symbol is 





N 


-V- 


u 


k. 












r^ 




where N is the switch number. 



and the analog circuit is 




-100 
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CHAPTER 3 
ELEMENTARY ANALOG PROGRAMMING 



All dependent variables in an analog computer are vol- 
tages. Consequently, it is necessary to equate these 
voltages to physical variables. Since the elements of 
the computer are voltage limited, scaling will be re- 
quired. The art of scaling is discussed in Chapter 5. 
In this chapter scaling is ignored since the necessary 
concepts can be developed without reference to scaling. 

As pointed out in Chapter 1, it is necessary to remem- 
ber that an anal(^ computer cannot differentiate easily. 
It is a machine whose main feature is its ability to 
integrate. Consequently, problems which are defined 
by differential equations must be reformulated in terms 
ojf integral equations either explicitly or implicitly. 
Consider the problem 

y = y(x) 

§ . By . 

y(0) = A. 

This can be rewritten as the integral equation 
y + Bfydx = 

y(0) = A. 
For this problem computer time will represent x, and 

y(t) + ByVdt = 0, t - X 

y(0) = A 

An integrator, as shown in Chapter 2, integrates from 
t = onward. It also requires an initial condition: 
namely the value of the integral at t = 0, which in 
this case is A. Thus, the problem is formulated for 
the computer as 

t 



y(t) 



= -jBy{t)dt, 



(1) 



y(0) = A 
The right half of equation (1) can be generated by 



where the constant B and the initial condition A are 
inserted by means of pots. Remember that the inte- 
grator inherently has a sign change . Equation (1) states 
that y(t) is equal to the negative integral. All that is 
needed, to complete the computer program for this 
problem, is to connect the output of the integrator to 
the input of the pot set to B. Hence, the computer 
prc^ram is 



A/100 



-100 




y(t),t = x 



The independent variable, y, can be recorded from the 
output of the integrator. 

In this example the differential equation was reformu- 
lated explicitly in terms of an integral equation. It is 
frequently possible to do this in an implicit way, as 
the next example will show. 

Consider the problem 

y + Ay + By = f(t), y = y (t) (2) 

y(0) = C 

y(0) = D 
This is equivalent to 

y = f(t) - Ay - By 

y(0) - C 

y(0) = D 
First, y can be generated from y by 



A/100 



-100 



+y(t)> 




-/ By (t)dt 




-100 
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(Note that again the pots are used to generate initial 
conditions.) Using this circuit, the complete program 
can be generated by 



D/lOO 



C/IOO 



-100 




where equation (2) is satisfied by applying the correct 
inputs to integrator 1. 



Next, the use of non- linear elements is demonstrated. 
Consider 



y + Ay + ty = cos t, 1 < A < 10, 
y(0) = C, 
y(0) = D. 
The computer program (except for scaling) is 
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As another example, consider 

y + f(t)y = 
^(0) = A 
y(0) = B 

The computer program is 



►+I00 



-100 



-y f (t) 




-100 



The next example shows how the DCU's (Chapter 8) can 
be usedto determine the solution of a differential equa- 
tion for a specified value of the independent variable. 
Suppose it is required to find y (. 17) for 

y + y = 

y(0) = 100. 

From above, the analog program which generates y is 



.-100 




In order to determine y(. 17) the computer is allowed 
to compute for 0. 17 seconds and then put into hold. 
The latter can be accomplished by applying +28 volts 
to the "problem hold" terminal at the patchboard after 
0. 17 seconds. This will halt the entire computer and 
in particular the computation carried out by the above 
prc^ram. The digital program which will do this is 
shown below: 




PROBLEM 
^HOLD 

TERMINAL 



<FRT^ 



The operation is as follows: The clock pulses are not 
applied to DCU 1 until the machine is in compute be- 
cause the output of the RH' bus is a logical in I C 
and this is an input to gate 1. The FF is initially in 
the R state when the machine is in compute having been 
reset by the FRT' Logic output. I^vider 1 reduces the 
clock frequency of 1000 cycles/sec to 100 cycles/sec. 
The output of the count of counter 4 is 10 cycles/sec 
and is connected to the input of counter 5. The 7co\mt 
output of counter 4 and the 1 count output of covmter 5 
are connected to gate 2. Thus, whenever t =» 0. 17 sec 
the output of gate 2 will charge from Oto 1, which will 
change the FF to the S state. At this time the computer 
will go to the hold mode due to the S output of the FF 
being connected to the problem hold terminal at the 
patchboard. Also, the clock and counters will be reset 
because of the logical 1 ^plied to the OV terminal. 
As soon as the machine goes to the reset mode the 
FRT' Ic^ic output will cause the FF to return to the 
R state. Thus, the logic circuit is ready to be usfitd 



again as the computer is manually put first in the reset 
mode and then the compute mode. 

ALGEBRAIC EQUATIONS 

The first examples were intended to illustrate the basic 
approach to programming. The followii^ illustrates 
a more serious application of the analog computer. 
Suppose it is desired to solve a set of simultaneous 
algebraic equations which are expressed in matrix form 
by 



AX = C 



(3) 



where the n x n matrix A and the column matrix C 
are known. For the purpose of illustration it will be 
assumed that A is of rank 2, althoi^h the derivation 
which follows is perfectly general. Consider 



X + AX = C 



(3a) 
3--3 



When the system has reached steady- state (assuming 
there is one) 

X = 

AX = C. 

Thus the solution to equation (3) is obtained from the 
steady- state solution of equation (3a). The latter can 
be written 



^1 ^ ^11^1 ^ ^12^2 ~ ^1 



^2 ^ ^21^1 "^ ^22^2 ~ ^2 



(4) 



and the analog program is 




022 



-100 



-100 



Generally, integrator capacitors are chosen equal to 
0. 001 iiid to decrease the solution time. However, 
there is no guarantee that the system of differential 



equations, (4), is stable. A different approach will 
provide an unconditionally stable system of differential 
equations . If 



it follows that 



AX = C 



A' AX = A'C 



(5) 



where A' is the transpose of A. The solution of equa- 
tion (5) is the same as that for equation (3), Again at 
steady- state the solution of 



X + A' AX - A'C 



(6) 



will be the solution of equation (5) and therefore equa- 
tion (3). The stability of equation (6) is determined by 



X + A' AX = 0. 



(7) 



If equation (7) is multiplied by the row matrix X, then 
the result is 

XX + XA'AX = 0. 

It is well known that XA'AX is positive definite (i.e. 
non- negative for all values of X). Unless X = 0, 

XX< 



and stability is guaranteed. The computer program for 
equation (6) is 




^<>-^-< 



"22 
O 



-100 




100 
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In this program it is easy to change individual coeffi- 
cients because each pot is associated with only one co- 
efficient. The program can be simplified so that it 
requires no more equipment than the previous pro- 
gram if the following coefficient changes are made. 



OLD 



NEW 



f (t) = A(t + a) ,• t > 0: 



df 
dt 


t + a 


f(0) 


= Aa 


Generating Equation 



■^11 



^11 "" ^21 



•^12 ^11^12 ^ ^21^22 



^21 ^11^12 "^ ^21^22 

2 + 2 

•^22 ^22 ^12 



^11^1 ^ ^21^1 



^12^2 ^ ^22^2 



This simplified program saves equipment but is in- 
convenient to use when the solution is required for 
several different values of the coefficients . 



100 



-o-lOO 




-100 



Program 



ARBITRARY CLOSED FUNCTIONS 

Many analytic or closed functions can be generated by 
representing them in terms of their generating dif- 
ferential equations. The following examples illustrate 
the technique. 

f(t) = 1/t, t>a: 

This fvinction has the generating equation 



df 
dt 



= -f 



f(a) = 4 



The program which solves the generating equation 



is 



(lOOa)" 




f(t) = 



-lOOo 



In (t + q), t > 0: 




df 1 




dt t+ a 




f(0) = Ino! 




Generating Equation 






Ina 
100 



-100 




Program 
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f(t) = A cosh(t), t>0: 



This prc^ram is 



f = f 

f(0) = A 

Generatii^ Equation 




-100 



A cosh(t) 



Program 



Mean value: 



The mean value, <f) (t), of a function, f (t), over the 
interval FT. , T- + tl is defined by 



T +t 
<t>it) =4/^ ^(^)^^' *^°- 
"^1 



(8) 



The generatii^ equation for ^ can be obtained by dif- 
ferentiation: 



|^.4{,(Tj*t)- ,(t)} 
The program is 



(9) 




*(t) 



-100 



Suppose the integrators 1, 2 are in the initial condition 
mode imtil t = T^, at which time their mode is changed 
to compute. Then, theoretically, the program will 
generate the solution of equation (8). However, the 
program will not work because the right side of equa- 
tion (9) is indeterminate for t = 0. Thus the divider, 
in the program will have an unstable output. Practi- 
cal limitations require the Y- input to be at least 3 
volts for useful operation. Consequently, it is neces- 
sary to choose a different initial condition: 



1 1 
(0 'T f ^(^)^^' 



■f (f) -*- 




"^ 



As before, when t' = T^, the mode of integrators 1, 
2 is changed from initial condition to compute . Initially 
SI is closed and S2 is open. The output of integrator 
1 will not be <f) (t) until t = « . When t = « , the cor- 
rect initial condition, namely ^(e), will have been es- 
tablished for integrator 1 . At this time the comparator 
will be activated, which will open SI and close S2. 
Also the output of integrator 2 will be € . Thus, the 
generating equation is implemented by the program 
starting when t = «' . The output of integrator 1 will 
be <p (t), t>«. 

COORDINATE TRANSFORMATION 

Suppose two coordinate systems have the same origin 
and are displaced by an angular rotation, e : 




An X- coordinate point, sayx, will have the coordinates 

X' = X cos e 

Y* = -X sin e 

in the x', yV coordinate system. Similarly, for a 
y- coordinate, say y, 

x' = y sin e 

y = y cos e. 
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Thus, a vector, R, in the X, Y- system with components 
X, y, will have the components 



positive rotations 0, 4>,\ respectively about the X, Y, 
Z axes correspond to the transformations 



X' = X cos e + y sin e 
y/ = -X sin d + y cos e 

in the X', Y' -system. In matrix notation 

R' = TR 

where R, R' are the representations of the vector in 
the \mprimed and primed coordinate systems respec- 
tively, and where 



r cos e sin di 
[jsin 9 cos ej. 



In the three-dimensional coordinate system 

Z 




^v 



R^ = 



R-^ - 



K- 



1 

COS 6 sin 6 
-sin e cos e 

cos </) -sin (p 

1 
sin (p cos <j 

cos X sin A 

-sin X cos \ 

1 



Any arbitrary three-dimensional rotation can be repre- 
sented as a product of these three matrices. The super- 
scripts refer to the axis about which the coordinate 
system is rotated and the subscripts to the name of 
the angle of rotation. 

The transformation of a vector, V, by R is 
V = RV 



V' = V 

X X 



V' = V cos + V sin 

y y z 



V = -V sine + V cos d. 
z y z 



The program is 




Note that 9 is limited by the range of the sine and cosine function generators. 
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The program symbol used frequently is 



V X X 

, Y y' 



The transformations R^, R have similar programs, 
MATRICES 

A matrix, A, is a collection of elements: 
i = 1, 



^ =Kj}j = i;.'.*,'m 



(in what follows it is assumed that n = m). Thus, if 
n = 3 



A = 



^11 
^21 



^12 



13 



'22 '^23 



31 32 



33 



Multiplication by a constant, a, follows the rule 



QiA 



"^1 


\2 


^3" 




"^^1 


aa^2 


^13" 


^21 


^22 


^23 


= 


^^21 


^22 


"=^23 


.^31 


^32 


^33_ 




."^31 


°^32 


"^33. 



which also can be denoted 



q;A = aa. 



aA 



afa..] = raa.-l 



when it is understood that the a. . are the elements of A. 
Addition follows the rule 

A + B = a.. + b.. = (a + b).. 



or for, say, n = 2 
A + B 



hi ^2 



^21 t>22 



''ll ^2 
^^21 ^^22 



hi ^^11 h2^^2 
^21 ■''^21 ^22 ^^^22 



Multiplication: 



Thus for, say, n 



AB = 



a. b . 



V a. D . 



AB = 



hl^ll ^ h2^21 
^21*^11 "" ^22^21 



hl'^12 ^ h2^22 
^21^12 ^ ^22^22 



Note that AB ^ BA. The unit matrix, I, is 

I, i = J 
0, i /^ j 



I = a.. = I, i = j 



Thus for, say, n = 3 



I = 



1 











1 











1 



I, has the property 



lA = AI = A 



Differentiation: 



dx dx 

or for, say, n = 2 



d . d (a..) da.. 

— :-A = -TT-^ ir = — — n 



-^A = ^ 
dx dx 



hi h2 
^21 ^22 



dx 



^hl ^^12 



dx dx 

^^21 ^^22 



dx dx 

Transpose: The transpose, A', of A is 



or for, say, n = 2 
A' 



a.. 



hi ^21 
h2 ^22 



Submatrices: The submatrix, Sy, corresponding to 
the element ay in A is obtained by 
deleting the i-tn row and j-th column 
of A. Thus for, say, n = 3 



'13 



"^11 
^21 



'12 *13 



"^22 *23 



•^31 32 



'33 



^21 h2 



^31 h2 



Determinant: The determinant, |a|, of A is a scalar 
quantity and can be found by induction. 
By definition 



N 
|a|= Z (-l^^a |S |, for any i 
j=l •* •• 



(10) 



(Note that some particular choice is 
made for i). 



For n = 1 



A = a 



11* 
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Thus, for n = 2, 
2 



l^h .^/-i) Sjl^ijh^Ai - ^ir 



12 



- ajja22 - ^^2^21 



or 



2 
|a|= I (-D^'^a |s I = -a2iS2^.a22S 



j=l 



22 



■^21^12''^22^11 



and either calculation leads to the same result. For 
n - 3, 



Cofactor: Thecofactor, C, corresponding to the ele- 
ment a... in A is 



C. = (-i)^''J(s..|, 

13 ^ ' iJ ' 



or for, say, n = 3 



13 



and 



'23 



^21 ^22 
^31 ^32 



^11 ^12 
^31 ^32 



and in particular, any of the equations 
3 



|A| . I^(-l)'*'ay|S,.| 



3 

I 
j=l 

3 

3 



1^1 = ^,^-'^^\\^i\ 



|A| = I M)^*'a3.|S3.| 



Inverse: The inverse, A' , of matrix, A, is 

A-' 
For example, if n = 2 



C. 

31 

TAI 



-1 



^11^22'^12^21 



^22 "^21 



■^12 ^11 



Orthogonality: A matrix, A, is said to be orthogonal 
if 



will lead to the same lA . If i = 1 



A| - ajjS^^-a^2Sl2"*'^13^13 



- aj^(a22ag3-a23a22) - ^i2^^21^33"^23^31^ 
^ h3^^21^32'^22^31^ 



^11^22^33"^11^23^32"^12^21^33 "^ ^12^23^31 



+ aj3a2ja32-a^ga22a3j 



Alternatively, | A can be determined from 

|a| = I (-D^-'Ja |S I, for any j. (11) 

i=l ^J ^^ 

A determinant can be multiplied by a constant, a, in 
a simple way: 

|qA| = a|A|. 

From equations (10) and (11) it can be seen that if a 
matrix B can be generated by m adjacent row inter- 
changes and p adjacent column interchanges starting 
with matrix A, then 



-1 



A' = A 



Note that R , R , R are orthogonal 
matrices. The geometric interpreta- 
tion of orthogonality is rotation. 

Logical Matrix Operations: 

(AB)' = B'A' 

(AB)"-^ = B'V-^ 

DECOMPOSITION OF ORTHOGONAL MATRICES 

Frequently the programmer is presented with a com- 
posite orthogonal (rotation) matrix which is a fimction 
of several angular variables. Usually these at^ular 
variables represent successive rotations about coordi- 
nate axes. The simplest way to generate the computer 
program is to decompose the matrix into a product of 
the RX, rY^ rZ matrices and then generate subpro- 
grams for each of these matrices . This can be accom- 
plished by following a simple set of rules. These are 
stated without proof and examples are given. R^, RY, 
rZ are said to be tjrpe 1, 2, 3 matrices respectively. 
Suppose the composite matrix, M, can be represented 



b| = (-1) 



m+p 



A. 



M = AjA2. 



^, 



3-9 



The outer matrices are A., A^ and the inner matrix 
Mj, is defined by i jn 



Mj = A2A3 



^.- 



Now suppose the simplest element (lowest degree homo- 
geneity) of M is located in row i and column j . Then A- 
is of type i andAj^is of type j. The angular variables 
corresponding to Ai, A^ will appear in the products 
of columns j, i respectively. The sense of the rotation, 
Ai, (positive or negative) can be foimd by setting all 
the variables corresponding to A2..., An equal to 
zero. The sense of the rotation, A^, is similarly 
found by setting the variables corresponding to A, . . . , 
An-1 equal to zero. If the simplest element of M is 
zero, then M has the representation 

M = Aj A2 

and there is no inner matrix. K the simplest element 
of M is not zero, then there is an inner matrix and Mj 
can be found by setting the variables corresponding to 
Ai, Ajj equal to zero. The further decomposition of 
Mj is carried out with the same procedure as for M. 



Example: 



Cos * 


Sin* 





Cos e Sin * 


Cos e Cos * 


Sin 9 


Sin d Sin * 


-Sin 9 Cos * 


Cos ( 



M = 



The simplest element is zero and is in row 1 and col- 
umn 3. Therefore there is no inner matrix and A^, 
A2, are of types 1, 3 respectively. The products of 
column 3 contain the variable 9 and the products of 
row 1 contain the variable for *. Therefore, the var- 
iable for Ai is 9 and the variable for A2 is «'. If * = 



M = 



1 








Cos 9 





-Sine 



If 9 = 



M = 




-Sin 9 
Cos 9 



Cos * Sin * 
■Sin * Cos * 



Finally, 



Example: 



M = 



M = R^R^ 

6 * 



Cos<^ Cos* + Sine Sin0 Sin* 

-Cose Sin 
Cos(f> Sin* + Sine Sin^ Cos* 



The simplest element of Mis Sin e and is in row 2 and 
columns. The products of column 3 contain the vari- 
able * and those of row 2 the variable <p . Therefore 
Aj, A3 are of types 2, 3 with variables *, cp, respec- 
tively. 



K e = = 0, then 



^1 = 



If * = e = 0, then 
If (^ = * = 0, then 



A2 = Mj = 



Cos * 


Sin * 



Cos 
-Sin 




-Sin * 

1 

Cos * 



Sin 

Cos (j) 

1 



Cos e Sin e 
-Sin e Cos e 



Finally, 



M = rY R^R^ 

-* e ^ 



As can be seen from the previous example, it is not 
necessary to find all the elements of Ai to determine 
the sense of *. Since it is known that Ai is of type 2 
it is only required to evaluate the element in row 1 and 
column 3 with 9 = <(> = 0. Similarly, if the element 
in row 1 and column 2 is found for * = e = 0, then 
the sense of ^ is determined. 

Example : (Here we use the abbreviations : C = Cosa, 
S = Sina) " 



M = 



c^V^-c^c^s^s^-s^c^c^ 



s^c c. +s s s s -c s c c^s,-s s,c 
<f) 9 ■^ 0e\* ex e * e x * 

S^SX -S.C„S,S,+C.C„C, S^S 4C s c 

(t> 9 ■^ (p 9 \ yf/ <p e \ 9 <i/ 9 \ -i 



-s c s -c s 

* X. * $ X. 



X * 



The simplest element of.M is C^C which is in row 3 
and column 3. The products of colinnn 3 contain e and 
those of row 3 contain . Thus, Ai is type 3 and Aj^ 
is also type 3. 9 , (^ are the variables for Aj, Aj^ 
respectively. Consider the element in row 1 and 
column 2. For X = * = ^ = its value is Sine and 
consequently 



\ - «-% 



Sin0 Cos* - Sine Cos0 Sin* -Cose Sin* 

Cose Cos<^ ,Sine 

Sin0 Sin* - Sine Cos<^ Cos* Cose Cos* 



For 9 = \ = 



* = its value is Sini^ and consequently 



\ ' 



= R 
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To determine the innermatrix, set 6 



so that 



Finally, 






s* 


^x 


V* 


-\ 


f^xS, 



c s 



The simplest element of Mjis a zero in row 1 and col- 
umn 2. Also for * = the element in row 2 and col- 
unm 3 is SinX. The products of column 2 and row 1 
contain the variables X , *, respectively. Thus 



M, = Rf r: 

I X * 



A = R^R^RyR^ 

-0 X * 



COORDINATE TRANSFORMATIONS WITH UNRE- 
STRICTED ANGLES 

The coordinate transformation program is the same as 
before, except for the use of Sine and Cosine function 
generators. For unrestricted range we use the vari- 
able 6 and the equations 



Sine = /eCose dt 
Cose ^-SeSmQ dt 



The program is, 



-100 





100 sin© 







J3 



+ 100 



The lower two multipliers and summer generate 

2 2 

E = -(100 - Sin - Cos 6). 

This error is used to change the damping of the inte- 
grators to ensure orthonormality of Sine, Cose. 






100 
00% 9 
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CHAPTER 4 
BLOCK PROGRAMMING FOR PHYSICAL SYSTEMS 

A Valuable Aid to Simulation of Lumped -Parameter Physical Systems 



The analog computer, or differential analyzer, is a 
powerful tool for obtaining dynamic solutions to differ- 
ential equations. Generally the equations to be sol- 
ved relate to a physical system of particular interest 
to the computer programmer, in which case the pro- 
grammer wishes to simulate his physical system by 
means of an electrical network whose defining equa- 
tions are analogous (hence the name analog computer) 
to the subject physical system. It is desirable then, 
to arrange the computer mechanization to have a one- 
to-one correspondence between the problem board and 
the physical system, so that the programmer "sees" 
his system rather than an abstract electrical network. 

The following discussion will illustrate a method of 
programming in which the physical elements are 
treated as blocks rather than as sets of differential 
equations. Certain useful basic building blocks will 
be developed, and methods for combining the basic 
blocks to form any desired physical system (simula- 
tion) will be discussed. Simulations generated by the 
methods to be described will have the advantage that 
any or all physical parameters appear as single con- 
trols in the computer so that any parameter may be 
changed without affecting any other setting. In addi- 
tion, any or all variables appear as computer outputs 
and may be individually monitored. 

A lumped-parameter system, by definition, consists 
of a number of interconnected discrete elements whose 
individual or collective responses to impressedforces 
or stimuli are of analytical interest. The forces and 
responses are related by a mathematical operator. 
Block programming starts with the selection of basic 
definition and rules. 

MATHEMATICAL OPERATOR (O) 

The operator, or combination of operators, describes 
the functional relationship between quantities. The 
most common operators are the following: 

Summer 

Constant Multiplier (potentiometer, commonly 
referred to as "pot") 

Inverter - - 1 Sign Changer 

Integrator - — (La Place Transform Notation) 

Differentiator - s (This operator is never used in 
an analog simulation if it can 
be avoided, as it usually can). 

Variable Multiplier 

Arbitrary Function - f(x) 

Analytical Function - Trigonometric, Log, 

Hyperbolic, etc. 

Inverse Operator - (O) s, — are inverses 



ELEMENT 

The element is the basic unit of the block representation 
of a system. It has a pair of terminals or nodes, with 
a value associated with each node. It also has a trans- 
ference quantity between nodes, or "through" the ele- 
ment which is functionally related by some operator to 
the node-pair value. The transference quantity may 
also be referred to as the branch transference or trans- 
mission. 

The elements required to depict the passive lumped- 
parameters of most physical systems are very few in 
number. In particular, five basic elements will sviffice 
for many linear electrical, mechanical, and thermal 
systems . These are the Summer, Constant Multiplier, 
Sign Changer, Integrator, and Differentiator. 

BLOCK 

A block is the computer mechanization related to a 
physical element. The input/output variables of a com- 
puter block are voltages and currents, but the operator 
relationship is the same as the simulated element. 

Table I gives examples of single elements and their 
corresponding blocks. Table II gives afew commonly 
occurring two-element combinations. Combinations 
of more than two elements can be generated as shown 
in the following examples . Useful tables of more com- 
plicated blocks or transfer-function simulations appear 
in many publications. ^ 

DRIVING SOURCE 

A driving source is a source of energy or power for a 
network of elements. 

SYSTEM STABILITY 

A stable system is one whose responses are bounded 
(finite amplitude limit) for any finite input. An un- 
stable system is one whose responses are not bounded. 

NETWORK 

A network is an interconnected set of elements and 
sources. Elements may be joined together at their 
nodes, providing the connected nodes have common 
dimensions and value, and provided the transference 
quantities have common dimensions. 

This discussion concerns elements and networks in 
which the following two rules apply: (Kirchoff's laws) 

1 . The algebraic sum of the node values around 
any closed path (loop) in a network is zero. 

2 . The algebraic sum of the transmission quan- 
tities at any node is zero. 



Korn and Huskey, "Computer Handbook' 
McGraw-Hill, New York, 1962. 
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TABLE I. Single- element Computing Blocks 



PHYSICAL 
ELEMENT 



SCHEMATIC 
SYMBOL 



MATHEMATICAL 
RELATIONSHIP 



PROGRAM BLOCK 



Summer 1 
Y = - (X^ + X2) 






D> 



Constant 
Multiplier 

(POT) 
Y = KX 



X (k) Y 



Inverter 
Y = -X 



Integrator 

s 



D> 



H> 



-Y 



Resistor 



Linear 
Spring 



Torsional 
Spring 



Heat 
Conductor 



Viscous 
Damper 



^I 



,-^^rr^( 



Differentiator 
Y =sX 



I = 



(Ej - Eg) 



R 



F = K(Xj - Xg] 



F = K(0^ - e^) 



K 
T,_J\/\/\/L_T2 
^Q 



Q =K(Tj - Tg) 



V| 



F=B(v^ - Vg) 







-Y 



-X, 




4-2 



TABLE I (Continued) 



PHYSICAL 
ELEMENT 



SCHEMATIC 
SYMBOL 



MATHEMATICAL 
RELATIONSHIP 



PROGRAM BLOCK 



Electrical 
Capacitor 



Viscous 
Damper 



Mechanical 
Mass 



Thermal 
Mass 



^^ 



(Ei-S'=ic 



H] 



B 



(Xj - Xg) - gg 



10—^2 



(^1 - ^2^ = iM 



— F 



|0-7-'"2 



^Q 



(Ti-T2)=^^ 




-(Xj-Xg) 




Electrical 
Inductor 



Mechanical 
Mass 




•"1 



1 = 



(El - Eg) 
sL 



i®— ^ 



X = 



(vi - Vg) 



-^X 




Y (k) S 



(Xj-Xg) 



COMPUTER PROGRAM MECHANIZATION 
COMMENTS 



1 . Operational amplifiers operate with the junc- 
tion at virtually zero potential; the currents 
into the junction through all the input and 
feedback paths must sum to zero; and there 
is as^ inversion between the amplifier in- 
put and output quantities. 

2 . Any c losed- loop in a problem mechanization 
simulating a stable physical system will gen- 
erally have an odd number of sign inversions. 

3. In mechanizing an equation, set the highest 
order derivative only on the left of the equals 
sign to obviate the necessity for differen- 
tiation. 

4. The mechanization is accomplished by: 



a) assigning aposition in the mechanization 
diagram to all necessary variables and 
operators. 

b) accomplishing the necessary intercon- 
nections so that each output has its pro- 
per inputs. 

c) checking to see that all input/output rela- 
tionships in the mechanization are satis- 
fied. 

EXAMPLES OF NETWORK PROGRAMMING 

The following examples illustrate some of the tech- 
niques used in applying the block programming meth- 
ods. Examples are given for linear elements only. 
Non-linearities can be included simply by replacing 
the potentiometer corresponding to the non- linear ele- 
ment with a multiplier, function generator, or other 
suitable computing element. Backlash, hysteresis, 
stiction, etc can be included as required. 
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TABLE n. Two- Element Computing Blocks 



SCHEMATIC 



MATHEMATICAL 
RELATIONSHIP 



PROGRAM BLOCK 



as 



(x^ - x^) 



Y - K (X^ - Xg) 
a- s 



— Y 



Hh 



-"2 
— I 



2 
•F 



H} 



M 



(E^ - Eg) 



^-r(^1-^2^ 
sC 



(Xj - Xg) = 



F - K (E^ - Eg) 
^B 



(Vj - Vg) = 



F - B (v^ - Vg) 

iM 



r€> 



(^-n> 



^»_^(X,-X2) 



K 



l/a-s 



Y = 



— Y 

R L 



(X^ - Xg) - KY 
a • s 



-*I 



M 



dJ 



1 = 



F = 



--F 



(E^ - Eg) - RI 
sL 



^^1 - ^2^ - B 
sM 




X, — a- 



l/bs 



Y = 



■*"Y 



^^1 - ^2) -a- s 
bs 



-^I 



M 



I = 



V = 



sL 



(^1 - ^2' - f 

sM 




-»v 
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MASS - SPRING - DAMPER SYSTEM 



r— X2,V2 




^P MOTION IN X 
^•^ DIRECTION ONLY 



V ///////////////// //////// 



Element - force equations 

Spring K^: F^^ = K^ (X^ - 0) 

Spring K2: Fg2 = Kg (Xg - X^) 
Damper B: F = B (Vg - v ) 



Initial values at t = 0. 
X^ (0), X2 (0), Vj (0), vg (0) 



Forces acting on masses 



^Ml = -^31" ^s2 -^ ^D 



F = -F - F + F 

*M2 ^s2 ^D 



Steps to mechanize 

1. Assume X J, -Xg, -v^^, Vg and F are available 

2. Draw blocks for M^ M^ 

3. Draw blocks for K., Kg, B 

4. Make input connections to mass blocks to produce X^ -Xg, -v^, Vg 




4-5 



I 



The above mechanization makes each parameter avail- 
able as an individual adjustment. The prc^ramming 
is easily accomplished by generating the forces im- 
posed on the elements and then operating on the forces 
to obtain velocity and displacement. Any or all of the 
initial conditions, Xj(0), X2(0), vi(0) and V2(0), may 
be zero, or may also be variable. 

Note that the block representing the damper contains 



two summers. The first summer produces (vj - V2), 
and the second summer acts only as a sign changer. 
The two- amplifier, one-pot combination can be reduced 
to a more simple two-pot combination by connecting 
Vj and V2 through B pots to the mass block inputs. 
The less complex alternate is shown below. This con- 
figuration requires fewer amplifiers, but variations 
inM^, M„, or B necessitate two or three adjustments. 




* Do not connect to -v- to get pot setting of B instead 

B 
of -TT- . A coefficient pot output cannot be connected 

to another coefficient pot input. This is a practical 
limitation of the computer hardware. 



Rigid Bar Supported on Springs 



— / ^ 






C.G 



UA 



T 



X2 



777777777777777777777777777777 



Ri^id bar of mass M, moment of inertia about C.G. of 
I. Assume small such that «sin B. Vertical motion 
only. Yq is displacement of C.G. 



* ~ si sM 



F, = -K, Yj 



Yl = Y,-.^a 



^2 " "^2 ^2 



Y2 = Yo+lgO 
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I~ LINEAR 
MASS 



SPRING I i SPRING zTrOTATIONAL I 

MASS 0-9 




One -Dimensional Heat Transfer 




K- THERMAL CONDUCTIVITY 

M-MASS 

C -SPECIFIC HEAT 



Ql = ^1 C^l - '^2^ ^2 = ^2 ^"^2 - '^3^ 



Q. 



• • 



sT, = 



1 MjC. 



sT„ = t. , ^ sT„ = ^ 



2 MgCg 



3 M3C3 



T 



-T, (0) 




-T, 



v^M.C. 



M,C, 



A 





K| 
-Q, 



V^MoOsI 



M2C2J 
61-02! 

A 





Om 



K2 
62 



3^3 



A 



L_ 



THERMAL 
MASS 



THERMAL! 
CONDUCT- 
ANCE 



THERMAL 
MASS 



THERMAL , 
CONDUCT- 
ANCE 



THERMAL 
MASS 
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Electrical R-L-C Networks 



II ^ R| 

-^O ^AA^ 



R2 

AAAr- 
-^13 



It^' t 



"1 



R3 



■l = 


^1-^2 


«1 


h ' 


E2-E3 


«2 


T 


E3-E4 


's 


SL 


^4 = 


'5«3 



^2 ^ ^1 - I3 



^4 = I3 - ^5 



'2 sC. 



'3 sC. 




L 



RESISTOR CAPACITOR RESISTOR CAPACITOR INDUCTOR RESISTOR 



C2 E3 




E1-E2 



sC. 



^2" S 



sC. 



^3 = ^3^3 



I„ = 



^2 
sL 



In = 



Il-^2 
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E3/R3 




TWIN-T NETWORK 



9z 



o— *— V\Ar"~* 



^l2 
C, 



♦ R2 



E2 



•R3 






2 " R. ' 3 sC„ 



"3 R„ 



(is-^V 



(^4 - ^2^ = sC 



11-^2 



E4 ^(Il-l2^^3^V^3 

E^ = (E^ - E4) + E^; E2 = E4 - (E4 - Eg) 




!.•> 
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CHAPTER 5 
SCALING 



The analog computer is voltage limited to a practical 
range of operation of ±100 volts as a maximum and ±1 
volt or zero as a minimum. It is usually not possible 
to numerically equate analog voltages to physical vari- 
ables. That is, the magnitude of physical variables 
must be scaled to fall within the useful range of the 
machine. This procedure is called amplitude scaling. 

The analog computer is also limited as to the speed with 
which it will solve a problem. Practical bounds on 
solution time are 100 seconds as a maximum and 100 
milliseconds as a minimum. Events in the physical 
world usually occur in time intervals which fall out- 
side these limits. Thus, the simulation of real world 
phenomena with the analog computer is ordinarily 
faster or slower. The procedure for relating computer 
time to physical time is called time scaling. 

Correct scaling is an important factor in reducing sim- 
ulation errors. An error analysis of any particular 
program is about as complex as the problem which is 
being solved by the program. Consequently, except 
for unusual cases, an analysis of error is not made. 
An estimate of the accuracy of the program can be 
found from check cases for which the answer is already 
known. Error is reduced by scaling the program so 
that all potentiometer values and amplifier gains are 
reasonable. Frequently, in the process of scaling a 
problem, it will be found that some parts of the model 
are not significant and can be eliminated. Scaling is 
a good check on reasonability. 

Scaling is an art, not a science. A good deal of exper- 
ience is required to become proficient. Scaling can be 
done in several ways . The techniques discussed in this 
chapter are those in most prevalent use among analog 
programmers. 

If X is a physical variable corresponding to computer 
voltage, V, then a scale factor, a, is chosen so that 



The general procedure for scaling is: 

1. Generate an unsealed program having a math- 
ematical structure which agrees with the pro- 
blem (model) to be solved. 

2. Identify the location of all voltages corres- 
ponding to physical variables. 

3. Associate a scale factor with each of these 
voltages (e.g. V^ " \^^* 

4. Label the program with the scaled physical 
variables (i.e. a X) rather than the voltages. 

5. Estimate maximum values of the physical 
variables . 

6. Choose the amplitude scale factors and the 
time scale factor, N, so that all pot settings 
and gains are reasonable. 

In order to accomplish step (4) it is necessary to know 
how input and output scale factors are related for 
various computer elements. 

POTENTIOMETER 



Physical equation: Y = bX 



Scale factors 
Scaled equation 

Program 



V = a X, V 
X X ' y 



a Y 

y 



V = b-^ V 
y a^ X 



a X 

X 



O 



a Y 

y 



b a 



V = a • X 

will fall within the practical range of the computer. It 
is good practice to make V as large as possible. If 
tp is physical (real) time corresponding to t^, the com- 
puting time, then these two can be related by a scale 
factor, N: 

t = Nt , dt = N dt 
c p' c p 



It follows that 



/(F)dt = N/(F)dt 



(1) 



N is chosen so that the problem solution time on the 
analog computer falls within reasonable limits . Prac- 
tically, this is accomplished by a choice of N which 
will not result in prohibitive amplifier gains. 



Note that a scale change can be made withapot. Assum- 
ing either a^ or a^ has been established previously, 
the other is chosen so that bax/ay is a reasonable pot 
setting. The pot can be used solely for a scale change. 
In that case b = 1. 

SUMMER 

Physical 

equation: Z = X + Y 

Scale 

factors : V = a Z, V = a X, V„ = a Y 
z z ' X X ' y y 



Scaled a^ a^ 

equation: v^ =— V^ + — V 
X y -^ 
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Now, if only one -gain inputs for the summer are used, 
then it must be true that a = a = a and the pro- 
gram is 



MULTIPLIER 



Physical 

equation: XY = Z 




a^Z 



Thus, it is seen that input scale factors must be the 
same for a one-gain summer and that a change of 
scale cannot be made through a summer. However, 
a different input scale factor can be used foralO-gain 
input. The program for this case is. 



-QvX 




a,Z 



Scale 

factors: V = a X, V = a Y, V^ = a^Z 

X X y y ' z z 

Scaled 

equation: V V 100a 

^ y ^ = \T 



Program: 



100 ' a a 
X y 




a^Z 



where it must be true that a 



X 



10a 



It must be true that 



INTEGRATOR 

As with the summer, the input scale factors (except 
for the 10-gain input) must be the same. However, a 
change of scale can be made with an integrator. Both 
amplitude and time scaling can be done. Scaling for 
an integrator is derived below, where C is the value 
of the capacitor in ^fd, and R is in megohms. 



100a 

2 

a a 
X y 



= 1. 



From this requirement it is seen that the relationship 
among the scale factors is the same as that for input and 
output multiplier voltages: 

a a 

X y 

a = 

^ 100 



Physical 

equation: Y =JXdt 

Scale 

factors : V = a Y, V = a X, t 



Nt 



Scaled a 

equation: ^ ^ '^ 
^ X 



^f 



V dt 
X c 



Thus, the integrator gain, ^^7 must be 



RC a^ 



(2) 



DIVIDER 



Physical ^ 
equation: — = Z 



Scale 

factors : V = a X, V„ = a^ Y, V^ = a Z 

X X ' y y ' z z 



Scaled 
equation: 100 



a a 

z y 

lOOa 



Scale 

factor a 

constraint: a = 100— 

z a 



and the program is 



-a«X. 




Program: 






• 


UxX 


• 




OyY 




a,Z 
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SQUARE ROOT 



Physical 

equation: Y = \/x 



Scale 

factors : V = aY, V =aX 
y y ' X X 

Scaled a 

equation: V = -^ . 10 /v 

^ 10 /i~ 



Scale 
factor 

constraint: a ==10 Va 

y ^ X 

Program: 



a,X. 



y~ 




FUNCTION GENERATOR 



QyY 



Estimates: 



\ p/max 



8 X 10 ; Xmax. 



The unsealed program is: 



\> 



-jdX 



M 



o 



M 



O 



i-X(O) 




<0 



The scaled program is obtained by assigning literal 
scale factors, determining their numerical values, and 
calculating the coefficient potentiometer settings. 



Physical 

equation: Y = f(x) 



V = aY, V =aX 

y y X X 



Scale 
factors : 

Scaled V 

equation: V = a f(— ^) 

y y a^ ' 



EXAMPLE; 

Mass- spring -damper system 



M 



dX. 
dtn 



HF 



K 



Physical 

equations: Spring Force, F = KX 



Damper Force, F _ = B 



dX 
dt 



d^x -(I's^^d) 



Mass 

Acceleration, ^^—^ - ,, 

dt 2 ^ 

P 

M = 1; B = 4 X 10"'*; K = 6. 4 x lO"^ 

^ (0) = 0, X(0)=.95 
P 



c .d£x_ 
N2 dto2 



I 




b d><, 
N dtn 



cB dX 



(sy 




N^Mdtp 



cB 
NbM 



■d)- 



cKX 



cK 



The —5 and TT factors appear 

N 
due to the inherent gain of N 
through each integrator, from 
Eq. (1). 



In general, the determination of the numerical values 
involves a trial and error approach. Any convenient 
parameter can be picked for the starting point. If a 
selected value leads to unreasonable potentiometer 
settings or vinreasonable amplifier gains, then new 
values may be required. Reasonable pot gains are 
between .1 and 1 . Reasonable amplifier gains are 
between 1 and 20. Amplifier gains should be integer 
values. Non- integer values can be obtained from a 
potentiometer, in combination with gain-of-10 inputs, 
if necessary. 

1. Determination of a. In order to utilize the 
full dynamic range of integrator 2, set 
aX max =100 (volts). 

a = 100 
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2. Potentiometer 1 setting. 



4. Potentiometer 2 setting; 



PI - aX(0) _ 



(3) 



At this time, a trial 



3. Determination of ^r. 

N 

value for N maybe selected. Analog solution 
times of . 1 to 100 seconds, or computer nat- 
ural frequencies from . 1 to 100 radians per 
second are the preferable operating ranges. 
An estimate of the physical time or frequency 
ranges can be used to select a value for N. 

For this example, a value of N = .01 appears 
reasonable. 



_b ^ 100 (volts) 

N (dX/dt )max. 
P 



(4) 



Looking ahead, 
require a gain 



note that integrator 2 will 

The gain of N is in- 

1. There- 



. aN 
°^ b" 



herent in the integrator from Eq. 

fore, the p^ value for the integrator will be 

-g-. It is highly desirable that this factor be 

an integer, preferably either 1 or 10. There- 
fore, b = 100 will be used. The effect of this 
is to force integrator 1 to utilize something 
less than its full dynamic range, but the effect 
is not serious in this example. As a general 
rule, it is not possible to scale so that all 
amplifiers work over their full range. 



5. 



P2 = -y^ = (6.4x 10"^)c 

N aM 



(5) 



The value of c is chosen to assure that a rea- 
sonable value for P2 is obtained and to assure 

that— is an integer. 

c = 100 will satisfy both conditions. 

P2 = .64 
Potentiometer 3 setting: 

cB 



P3 = 



NbM 



= .04 



This setting for potentiometer 3 is below the 
desired minimum potentiometer setting value 
of . 1. In this example, however, it is not pos- 
sible to have a large setting for potentiometer 3. 
This follows from a consideration of the phy- 
sical problem. The damper, B, was chosen 
to have a very slight effect on the mass-spring 
system, so it necessarily follows that the 
portion of the analog program related to the 
damper will have only a slight effect. It should 
be noted that abnormally low potentiometer 
settings will frequently arise in simulations of 
systems containing negligibly small elements. 
It should be recognized that the low settings 
stem from the actual system characteristics. 
The programmer, therefore, need not spend 
time in futile attempts to improve the scaling. 



The final program with numerical values is: 



-100 




X= 



V2_ 
100 
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GENERAL COMMENTS ON TIME SCALING: 

The following points concerning time scaling are of 
particular importance: 

1. If all integrator gains are simultaneously 

changed by a factor k, the time scale changes 

N 
from N to -^, and all computer frequencies 

increase by k. The amplitude scaling, how- 
ever, does not change. All numerical scale 
factors, all potentiometer settings, and all 
computer problem voltages remain the same. 
This property of amplitude invariance under 
time scale change allows the programmer to 
easily speed up or slow down the computer 
solution. 



The most convenient way of changing all inte- 
grator gains simultaneously is by changing 
the capacitors by a fixed factor. Integrator 
capacitors can be simultaneously changed by 
means of the computer Time Scale control^. 

An event which occurs in a fixed real-time 
interval can occur in several different com- 
puter time intervals depending on the choice 
of integrator capacitors, that is, depending 
on the time scale factor. 

Events of identical character which occur in 
different real-time intervals can all be made 
to occur in the same computer time interval 
by appropriate time scaling. 



On the SD 10/20 and SD 40/80 analog computers, a master Time Scale switch is a standard feature. This per- 
mits fast and convenient switching of capacitors that are patched -- up to 1000:1 time scale (xlO,xlOO, xlOOO, 
depending on patch panel connections). 
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CHAPTER 6 
COMPUTER OPERATION 



This chapter is based on materials from the Systr on- Donne r 10/20 and 
40/80 Computer Operating Manuals. The descriptions of the computer 
mode controls, computer logic, and various computing modules permit 
the reader to relate all applications and programming examples given in 
this handbook to an existing computer. For example, Chapter 8 offers 
detailed information on the use of digital logic. The modules needed to 
implement the applications described in Chapter 8 are shown in this 
chapter. 

Not all available S-D computing components are described here, only 
those which illustrate and serve to explain a basic function. For the 
latest, up-to-date availability of all computing elements and modules in 
the SD 10/20 and 40/80 series, please consult Systron-Donner or your 
nearest SD engineering sales representative. 
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SD 10/20 Computer 



The Model SD 10/20 is an all solid state analog com- 
puter with an operating range of ±100V and expansion 
capabilities for up to 24 operational amplifiers. The 
principle features of the computer are full ilOOV com- 
puting range, all solid-state design, compact modular 
construction of all computing components, removable 
problem board which mates directly with computing 
modules to eliminate all problem board cabling, a 
control center with readout selectors for amplifiers 
and potentiometers, high-speed reed relays and solid- 
state switches, complete logic capability for repetitive 
and iterative operations, 0.01% accuracy of computing 
resistor and capacitor networks, a panel voltmeter 
having a null mode of measurement with 0.02% (of full 
scale) resolution. The amplifiers are protected against 
accidental grounding of their outputs as are all power 
supplies. 



All computing components in the Model SD 10/20 are 
modular plug-in units. They are identical to those 
used in the larger SD 80. This high degree of standard- 
ization between SD 10/20 and 40/80 computers maxi- 
mizes usefulness of computing modules. Since all 
modules are standard plug-in units, the computer user 
is offered a wide choice of computing components. 
As new modules become aval loble, the SD 10/20-40/80 
series computers can thus be constantly up-graded to 
include the new modules. Furthermore, a user may 
start with an SD 10/20 and later use the same 10/20 
modules in an SD 80. This common design feature 
results in important cost savings. As the need for 
more computer modules increases, the expansion cost 
from an SD 10/20 to an SD 80 is minimized. 




m Panel Meter Ranges: 1, 3, 10, 30, 100, and 300 volts, and 
± null. Full scale accuracy: 3%. Null position provides 
0.02% F.S. resolution with reference potentiometer having a 
±0.05% linearity at 25°C. 

m Compute Time Selector — compute time continuously vari- 
able from 5 msec to 10 sec. Reset time varies from 5 msec 
to 5 sec, depending upon coarse steps of compute range. 

m Meter Range Selector — with positions for 300 v, 100 v, 30 v, 
10 v, 3 v, and 1 v. Serves also as sensitivity adjustment for 
± null. 

m Function Selector — for rapid choice of: + null, — null, 
Meter, External {connects selected bus to external jack). 

m Address Selector — address of all amplifiers and potentio-^ 
meters. 

m Null Reference Potentiometer — provides high accuracy 
readout using null method with 0.02% F.S. resolution. Lin- 
earity is ±0.05% at 25°C. 

[t] Function Switches — provide manual switching flexibility 
in problem solutions. 



Ij] Mode Selection (lighted pushbuttons): 

Hold — places problem solution on all integrators into 

hold position. 
Compute — applies problem voltages to all integrators. 
Reset — applies initial condition voltages to integrators. 
Rep-Op — places integrators into a repetitive operation 

cycle. Compute time variable from 5 msec to 10 sec. 
Bal/Pot Set — disconnects junction and grounds the 

input resistor summing junctions of all amplifiers. 

Each amplifier is converted to a gain of 2500 for 

precision monitoring of junction offset. 
Pwer On/ Off — energizes and de-energizes computer. 
Oven — indicates + 28-volt oven power is on to maintain 

constant temperature of computing capacitors. 
OL Hold, OL Reset — lights up when any amplifier is 

overloaded. When depressed, computer goes into 

Hold; when released, normal operation is resumed. 
Time Scale — activates relays in each integrator module 

to change computing capacitor, (x 10, x 100, x 1000, 

depending on patchpanel connections.) 
Slave — permits operation of computer control circuitry 

from a second console. 

m Coefficient Potentiometers — up to 24, available in groups 
of 6, featuring 10-turn wire-wound with iockable counting 
dials. 
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SD 40/80 Computer 



The same compact design and high performance com- 
puting capability found in the small SD 10/20 computer 
is also available in the large SD 80 computer. With an 
expansion capability to 84 amplifiers (126 special 
order), 125 potentiometers, and a full set of digital 
logic control modules, the SD 80 provides a problem- 
solving power capable of handling large, complex 
problems normally associated with the simulation of 
dynamic electrical, mechanical, thermal and similar 
systems. 

Yet, the large-capacity SD 80 fully retains the 
versatility of a desk-top computer. Mounted on a desk 
or table equipped with lockable casters, it can be 
easily moved to new locations. In addition to creating 
a highly compact and rugged computer, the all solid 
state design increases reliability, simplifies mainten- 
ance, and eliminates the need to operate the computer 
in an air conditioned laboratory environment. 

Programming a problem on a fully expanded SD 80 
is virtually as simple as operating a small SD 10/20. 
Visual computer circuits in symbols everyone under- 
stands show at a glance how a circuit is patched. 
This feature saves time, reduces patching errors, and 
makes patching a task that the non-expert can perform. 

To allow maximum operator control over a computer 
with a capacity as large as that of the SD 80, all com- 
puter controls are centralized in the left wing which 



is movable to any convenient position to suit the 
operator. Major SD 80 controls that speed and simplify 
problem checking and solving include: 

4 A 199-point digital address pushbutton selector 
system for rapid access to amplifier outputs and 
potentiometer arms. 

A Static and Dynamic checks by simple pushbutton 
control to verify patching and test time scaling 
rapidly. 

4 Instant alarm controls, both visual and audible, 
with automatic transfer to HOLD at the occurrence 
of an overload. 

4 Individual control of compute and reset times for 

precise adjustment in Iterative and REP-OP modes. 

* Time scaling the problem with a single switch on 
the Control Wing. 

4 One pushbutton to check any amplifier's balance 
without repatching and to set any potentiometer. 

The problem board is removable, just as it is with 
the SD 10/20. There is no restriction as to grouping 
of the modular patch panels in the overall problem 
board layout. In both the SD 10/20 and 40/80 compu- 
ters, digital logic control is provided by modular plug- 
in units that slip into the central patchbay. They thus 
form an integral part of the computer and the main 
problem board, not a separate added-on element. 




1. Digital Voltmeter — 5 digit readout, automatic ranging and 
polarity. 

2. Address Selector Panel — with an amplifier and potentio- 
meter address capability of 199 points. 

3. Control Panel — for selection of: 

Op. Modes: Hold, Compute, Reset and Rep-Op. 

Master Time Scale Switch. 

Compute Time: 5 msec to 5 sec. 

Reset Time: 5 msec to 5 sec. 

Slave: for remote control. 

Oven and Power pushbuttons. 

Check Mode switches for Static and Dynamic checks. 

Overload Controls for: 1) audio alarm, 2) transfer to hold 
mode, 3) removal from Hold mode. 

Amplifier Balance and Pot Set: changes all amplifiers to 
gain of 2500 for simplified monitoring of junction off- 
set; and for proper loading of potentiometers. 



4. Function Switches 

5. Coefficient Potentiometers, No. 1-5 

6. Null Meter Reference Potentiometer 



7. Removable Patchboard. — An SD 80 board can hold up to 

42 individual computing module panels. All modules 
are mounted in a universal patchbay directly behind 
the removable problem board. Nomenclature on patch 
panels matches textbook circuit diagrams. Operator 
can thus "see" the circuit while patching a problem. 

8. Built-in Digital Logic Modules, consisting of Flip-flop, 

Gate, and Time/Event Control Modules. 

9. Potentiometer Wing can hold up to 6 panels for a total of 

125 potentiometers (first 5 pots are on Control Wing.) 
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Computer Logic 



Each operating mode switch energizes a logic buss 
which is then distributed to the logic relays in the com- 
puting modules (primarily to the integrators). The logic 
levels in the Analog sections of the computer are: 

+28V = True = logical 1 
OV = False = logical 

A logical one energizes a buss and places the com- 
puter in that mode of operation. All contradictory modes 
have their switches interlocked to prevent two modes of 
operation from being energized at the same time. 

The table t+iat follows outlines the logic busses and 
operating modes. 

In most cases the nomenclature for each logic buss 
is self-explanatory. The following explanations should 
clarify the remainder of the logic nomenclature. 



Logic Buss 



Description 



F. R. (R. 0.) Forward Reset (Repetitive Operation). 
This is the normal Reset buss and also 
is the Reset buss energized in the reset 
portion of Repetitive operation. 

R. R. (R. 0.) Reverse Reset (Repetitive Operation). 
This buss is the logical complement of 
forward reset. It is used for complemen- 
tory integrators and track/store circuits 
in iterative computation. Energized in 
the Compute mode of operation and in 
the Compute portion of Rep-Op. 

F. R. (R. T.) Forward Reset (Real Time). This buss 
is used to reset non-repetitive operation 
integrators. It is energized only when 
the reset switch is engaged. 

F. H. Forward Hold is the normal Hold buss. 

It is energized in the Hold mode of 
operation or if an overload occurs when 
the Overload Hold switch is engaged. 

R. H. Reverse Hold. Reverse Hold is the logi- 

cal complement of F.R.R.T. It is used for 
complementary integrators and Track/ 
Store circuits in iterative computation. 
Energized in the compute mode of opera- 
tion and in Rep-Op. 



Table 1. CONTROL LOGIC TABLE 





COMPUTER OPERATING MODE 








LOGIC 
BUS 


POT-SET 
BALANCE 


RESET 


COMPUTE 


RtfP-OP 


MANUAL 
HOLD 


PROBLEM 
HOLD 


O.L. 
HOLD 


STATIC 
CHECK 


DYNAMIC 
CHECK 


TIME SCALE 
CHANGE 


Bai 


+28 
























Pot -Set 


+28 
















+28 








F R (R ) 





+28 





0, +28 

















R.R.(R.O.) 








+28 


+28, 

















F.R (R.T.) 





+28 





















F.H. Bus 














+28 


+28 










R.H. Bus 
Static Ck. 








+28 


+28 









+28 
+28 



+28 




Dyn. Ck. 
Time Scale 




















+28 



(1) Blank indicates either 1 or +28 condition can exist, two values indicate cyclic operation. 

(2) Manual Hold is actuated by pushbutton switch on Logic Control Panel. p_uu„ Board 

(3) Problem Hold is actuated by application of +28V or greater Problem Hold Trunklme on Problem Board... 

fNote- Will not be damaged by any voltage in computer range) . „„„lifi-, in 

(4) O.L Hold is actuated by operating "CL. Hold" switch with an O.L. condition (either temporary or permanent) or any amphf.er m 

the computer. 
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Use of the Mode Controls 

The mode controls are of the push button type with 
interlocks to prevent two contradictory modes from being 
energized at the same time. 

All control in the computer is accomplished in the 
integrators. Thus, when a mode is energized it opens 
or closes a relay in each integrator. Reference to the 
illustration, Integrator Block Diagram, on following page 
(6-6), while reading the following explanations will 
assist in understanding the control system. 

Reset 

The Reset or Initial Condition (IC) mode enables 
the integrators to be set to their initial values before 
starting computation. When the Reset relay is energized, 
an input and a feedback resistor are associated with the 
integrator amplifier making it, in effect, an inverter. 
Applying a voltage to the input resistor (in this cose the 
IC terminal) will cause that input to appear at the output 
and to simultaneously charge the integrator capacitor 
to the same voltage. 
Compute 

In the compute mode there are no energized relays. 
The logic rs such that if the computer rs not in Hold and 
not in Reset, it rs in compute. 

Hold 

In the Hold mode a normally closed relay contact is 
operted which disconnects the input, signal from the 
amplifier. As there is no discharge path for the inte- 
grator capacitors the output will remain at its last com- 
puted value. In thrs manner outputs may be read or re- 
corded and the problem resumed at will by returning to 
Compute. 

Rep-Op 

In the Rep-Op (Repetitive Operation) mode the com- 
puter is automatically switched between the Reset and 
Compute modes of operation. This mode enables the user 
to obtain repetitive solutions which are easily monitored 
on an oscilloscope. 

The amount of time the computer stays in the compute 
mode is determined by the setting of the Compute Time 
switch located at the lower left of the control center. 
The outer dial of the switch may be positioned in any of 
three ranges: .005 to .1 sec, .05 to 1 sec, and .5 to 10 
sec or in the External (EXT) position. The inner dial 
adjusts the compute time over the range selected. 

The compute and reset times may be adjusted arbi- 
trarily by observing the solution to the problem and ad- 
justing the times until the desired display is obtained^. 
If more precise compute and reset times ore required the 
Compute-Reset cycle may be observed on an oscillo- 
scope. Connect the oscilloscope to any of the FOR- 
WARD RESET terminals on the integrators or to Test 
Point 6 at the rear of the computer. A rectangular wave- 
form will be observed. The positive (+28V) portion is 
Reset, the negative (OV) portion Compute. The Reset 



and Compute time controls may then be adjusted to give 
the desired times. 

If the Compute Time switch is placed in the External 
(EXT) position external circuits may be used to generate 
the Rep-Op cycle. The drive signals must be compatible 
with the logic used in the computer (+28V = True, OV = 
False) and should be connected to the.R (reset) and 
C (Compute) terminals located on the problem Jjoard. 

OL Hold - OL Reset 

When engaged, the Overload Hold portion of this switch 
will cause the computer to automatically switch to the 
Hold mode at the occurrence of on overload in any ampli- 
fier. When the overload is found and eliminated, the 
switch is re-engaged toreset theoverload hold circuitry. 

If an overload occurs, the light within the OL Hold 
/OL Reset will begin to blink on and off and will con- 
tinue to do so until the overload is removed. 

Time Scale 

The Time Scale switch determines which of the 
integrator capacitors connected to the A and B terminals 
on the integrator modules is in the circuit. With the Time 
Scale switch not engaged, the capacitor connected to 
the A terminal is in the circuit. When the Time Scale 
switch is engaged the capacitor connected to the B ter- 
minal is in the circuit. Any two of the three available 
capacitors may be patched to the A and B terminals. If 
time scaling is not required the capacitors may be patch- 
ed directly to the output. 

Special Trunklines 

There are six special trunklines located on the pro- 
blem board of the Computer. These ore located on the 
top left row of modules. The first module contains the 
M, X and Y special trunk lines and the third module 
contains the H, C and R special trunklines. The special 
trunklines are used as follows: 

M; This line is connected to the panel meter or the EXT 
jack as determined by the position of the Meter Select 
Switch. If the Meter Select Switch is in the METER 
function the M line is connected to the EXT jack; if 
the Meter Select Switch is in the EXT function, the 
M line is connected to the panel meter. 

X & Y: The X and Y lines are connected directly to 
the X and Y terminals on the left side panel of the 
computer. These lines are intended for use with an 
X-Y plotter but may be used to connect signals to 
any external device or they may be used to connect 
external inputs into the computer. 

H, C & R: The H (Hold), C (Compute), and R (Reset) 
lines are used as external input lines to the Control 
Logic circuitry. A logical 1 applied to any of the 
lines will place the computer in that mode of opera- 
tion. To use these lines the computer must be in 
REP-OP and the Compute Time switch must be in 
EXT. 
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ntegrator Modules 








STAT. 9-"° ^ 

CHECK * 



IM 



Integrator Block diagram illustrating flexibility 
of circuit and patching convenience. 
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Model 3320 

al integrator with two uncom- 
mitted operational amplifiers. 



Moc/e/ 3320, Dual Integrator 

The Model 3320 is divided into three sections: the 
top and bottom integrators and a logic section common 
to both. Basic patching of the integrator requires only 
a patch from the end of the integrating capacitor to the 
A terminal of the time scale relay. By patching the end 
of another capacitor to the B terminal a second time 
scale is set up and is controllable by the TIME SCALE 
switch on the control panel. 

Initial conditions (IC) are applied at the IC terminal. 
The source of the IC voltage may be a pot, amplifier 
output, the reference voltage, or from an external source. 

The logic section consists of a 7-terminal matrix 
near the center of the patching block. The nomenclature 
is the same as that explained in the Logic Control sec- 
tion of this manual. The two center terminals of the matrix 
are connected to the coils of the Reset and Hold relays. 
For normal and repetitive operation, patch from the F 
(forward) terminals to the relay coils. For reverse and 
iterative operation, patch from the R (reverse) terminals 
to the relay coils. For nonrepetitive operation, patch 
from the F(RT) terminal to the reset relay coil. (In this 
mode the integrators will reset ONLY when the reset 
button on the control panel is pushed.) Note that the 
logic patching controls both integrators in this module. 

To use the integrator as a summer, patch from the 
input side of one of the resistors to the output of the 
integrator. This will connect a resistive feedback around 
the amplifier; the other resistors may be used as summer 
inputs. 



Model 3329 

Quad integrator, dual multiplier, 
operational relay. Contains four 
integrating networks, two four 
quadrant multipliers and one 
operational relay. Also avai lable, 
less multipliers. Model 3329A. 




Mode/ 3329 - Quad Integrator, 
Dual Multiplier, Function Relay 

The Model 3329 has four integrating networks without 
operational amplifiers. The operational amplifiers neces- 
sary for integration are normally patched into the Model 
3329 from an adjacent Model 3325. The input resistor 
network of the Model 3325 then becomes the input to 
the integrator. Relay logic connections are in the lower 
left hand section of the Model 3329 patchboard. Nomen- 
clature is the same as that used on the Model 3320, Dual 
Integrator. Each Integrator Network has a Hold (H) and 
Reset (R) terminal associated with it, these terminals 
must be individually patched to the Relay Logic termi- 
nals. The Squid (multiple) patch cords should be used 
when more than one integrating network is to be used 
at the same time. 

The Dual Multiplier section of the Model 3329 may 
be patched for division, squaring and square root func- 
tions as well as for multiplication. For multiplication 
or division operation, three operational amplifiers are 
required. In the division and square root functions arbi- 
trary limitations as to polarity and relative amplitude of 
the inputs are imposed. For example, in division the 
Y input must always be positive and equal to or greater 
than the X input. Note that in square root an external 
diode must be patched into the circuit. 

Patching and operation of the Function Relay in the 
Model 3329 is identical to that of the Function Relays 
in the Model 3322A. The relay is energized by a nominal 
+28 volt signal and will not be damaged by a signal of 
up to ±100 volts. 
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Summer Modules 
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Model 3321: Dual summer with 
two uncommitted operational 
amplifiers. 



Model 3325: Quad summer, four 
summing networks, with four 
uncommitted operational 
am pi if ier s. 




Model 3321 - Dual Summer 



Model 3325 - Quad Summer 



The Model 3321 contains two operational amplifiers and 
the components necessary for summing operations. Basic 
patching consists of connecting one end of the feedback 
resistor to the output. The top summer has three input 
resistors, the bottom summer two input resistors. There 
are four spare summing resistors located in the center 
section of the module which may be used as either input 
or feedback resistors with either of the summers. 

The module also has two 20 kohm resistors and 
two uncommitted diodes associated with each summer. 
These may be used to patch the summer os a limiter- 

comparator or for any of the common diode circuits used 

in computation. 

A function switch termination is provided on each 
summer module. The function switch is a single pole- 
double throw type switch and is physically located near 
the control center. Normally only the first three summing 
modules will have switches terminated in them. 

The center section of the summing module has ter- 
mination for a diode function generator. Use of these 
terminations and of the diode function generator is ex- 
plained in that section of this manual. 



The Model 3325, Quad Summer is identical in per- 
formance to the Model 3321, Dual Summer. As illustrated, 
patching of the Model 33^5 is, however, considerably 
different than that of the Model 3321. Basic patching for 
any of the four summers in the module consists of patch- 
ing the SJ' terminal to the J terminal and patching the 
right end of the feedback resistor to any of the output 
terminals. Any of the resistor input networks may be 
used with any amplifier. Thus summers with six, nine 
or twelve inputs may be patched. 

The Model 3325 also contains four spare resistors, 
four diodes and four trunk lines. In addition, it has 
terminations for a Function Switch and for a Diode 
Function Generator. 

The Model 3325 is also used as a source of amplifiers 
for the Model 3329. When used in this manner, the Model 
3325 furnishes the amplifier and input resistors for the 
integrator circuit. 
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Model 3322A - Dual Inverter, Dual Funcfion Relay 
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Mod«l 3322A: Dual inverter, dual operational relay. Two 
uncommitted operational amplifiers with 
two operational relays. 



The Dual Inverter contains two operational amplifiers 
for use as inverters or high gain amplifiers in comparator 
and limiter circuits. Operation as an inverter requires 
two patches: (1) From the summing junction to the junc- 
tion. (2) From the input side of one of the resistors to 
the amplifier output. The inverter may be used as a 
summer by patching additional resistors into the sum- 
ming junction. (Note: Potentiometers used as inputs 
to the Model 3322A should be set in the Reset mode of 
operation.) 

The Model 3322A also contains two function relays. 
Each is a two form C (2FC) relay corresponding to a 
double pole-double throw (DPDT) switch. The relays may 
be energized by any voltage from +28 volts to + 100 volts. 
The energizing source is typically the output of an amp- 
lifier or a logic signal, but any source within the +28 
volt to +100 volt range will energize the relays. Any 
voltage in the range plus and minus 100 may be applied 
without damage. 

As in the Dual Summer, the Dual Inverter contains 
uncommitted diodes, limiter resistors, terminations for 
a Function Switch and Diode Function Generator ter- 
minations. The use and operation of these elements is 
the same as that described in the section on the Model 
3321. 
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Moc/e/ 3323 ~ Dual Inverfr, Dual kiultipliar 
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Moc/e/ 3323: Dual multiplier, two four quadrant quarter 
square type multipliers with two uncom- 
mitted operational amplifiers. 

I^ulti plication Operation. The multiplier /divider is 

based upon the quarter-square principle, in which the 
following relationship holds: 

1/4 [ (X + Y)2 - (X-Y)2] = XY 

In ±100 volt computers, the multiplier output voltage is 
always scaled to -0.01 XY. 

The circuitry consists mainly of two diode squaring 
networks, one based upon positive polarities and the 
other upon negative. For four-quadrant operation, volt- 
ages representing X, -X, Y, and -Y must be present at 
the inputs to the squaring networks regardless of the 
polarity of input voltages. 

The center section of the module contains the termi- 
nations for the multipliers. The section is divided into 
two identical halves, one for each multiplier. The multi- 
plier will also perform division, squaring, and square 
root operations 



'"'"'^"O 




Division Operation. Division is accomplished by the 

terminal connections shown in the figure below. The 
concept is illustrated by the circuit diagram in the same 
illustration. If the output of amplifier 1 is designated Z, 
the output of the multiplier must be+.OlYZ. By reference 
to the basic discussion on operational amplifiers (para- 
graph 2.8), it follows that X = -.OlYZ, or Z = -lOOX/Y. 
Since +.01YZ and Z must always be opposite in polarity 
to X, it also follows that Y must always be positive in 
order to maintain thrs relationship. 

If Y occurs in the problem only as a negative voltage, 
apply it to the -Y multiplier terminal and patch the input 
and output terminals of the inverting amplifier to the -Y 
and Y terminals, respectively. Z still equals -lOOX/Y. 
The absolute value of Y must always be larger than X; 
otherwise, the output would tend to be larger than 100 
volts, causing an overload. 



-0.01 xz 
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O Y (ALWAYS +) 
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INSERT 




CONNECTIONS FOR MULTIPLICATION 



CHANGE ONLY THESE 
CONNECTIONS TO CON- 
VERT FROM MULTI- 
PLICATION TO DIVISION 



CIRCUIT FOR DIVISION OPERATION WITH INSERT SHOWING 
HOW TO CONVERT FROM MULTIPLIER PATCHING. 



Squaring. Squaring is performed as a special case 

of multiplication where X = Y. Only one input inverter 
amplifier Is required. Patch together the R and Y input 
terminals and the -R and -Y terminals. Leave the other 
connections as for multiplication. The output voltage 
taken at the amplifier output terminal represents -0.01X2 

Squareroot Operation. The squareroot operation is a 

special case of division where Z = Y. The relationship 
Z = -lOOX/Y of division becomes Z2 = -lOOX for -100 

<x<o, z= loV^, foro<x <+ioo, z = -io>/5r 

If X is always negative, patch together the R and Y 
terminals and the -R and -Y. If X is always positive, 
patch R to -Y and -R to Y. Leave all other connections 
OS for division. 

The Inverter sections of the Model 3323 are similar to 
those described for the Model 3322A. The two basic 
patches required for inverter operation ore: From the 
summing junction to the junction and from the input 
side of one of the resistors to the output of the amplifier. 
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Moc/e/ 3324 - Dual Inverter, Quad Electronic Switch 




Model 3324: Electronic switch, four high speed elec- 
tronic switches with two uncommitted 
operational amplifiers. 



4100 -100 




The center section of the Model 3324 contains four 
electronic switches. There are two normally open (NO) 
and two normally closed (NC) switches. The normal volt- 
age required to excite the switch is a nominal +28 volts; 
however, voltages of up to ±100 volts will not damage 
the switch. The current required to excite the switch is 
less than three milliamperes at 28 volts. 

The Electronic switches were designed and are nor- 
mally used to replace the reed type relays in the inte- 
grator modules. The Electronic switches may be used as 
normal switches as long as the input voltage and current 
through the switch limitations are observed 

Patching to the Model 3320, Dual Integrator is shown 
in the accompanying illustration. Patching to the Model 
3329, Dual Integrator is similar. Also included in the 
Model 3324 are two .001 mfd. capacitors and four .01%, 
10 K ohm resistors for use in high speed Sample-Hold 
circuits. 

The inverter sections of the Model 3324 are identical In 
operation and performance to those of the Model 3322A. 
Patching is identical to that of the Model 3322A and the 
instructions for that Model may be used. 



IN O— 



ELECTRONIC 
SWITCHES 




INTECRATOR 



Integrator with Electronic Switches 



Electronic Switch Block Diagram 
(Connections shown for comparator amplifier) 
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Digital Logic Modules 




Model 3326 

Twelve RST Flip-flops with 
four relay drivers. 




Model 3327 

Twelve three-input AND 
Gates with four relay drivers. 



Digital control, the new way of multiplying the efficiency 
of an analog computer, can be included in the SD 40/80 as 
well as the small 10/20 computers. The SD hybrid computer 
expansion system is comprised of three types of compact plug- 
in modules. 

The advantages gained by digital control in an analog com- 
puter are of far-reaching significance. Here are some important 
new advantages made possible by SD's Digital Logic Control: 

1. Track and hold operation by individual integrators. 

2. Sub-routines can be flexibly programmed at different speeds de- 
pending on decisions made by logical equations. 

3. Program statements can be arranged into a flow chart quite sim- 
ilar to those used in digital computation. 

The flexibility gained through this interplay of analog/ 
digital equipment results in: 

1. Better and greater problem-solving capacity. 

2. Ability to solve a wide range of problems that before could not 
easily be handled by an analog computer. 

3. Speed. Problem solving time is greatly reduced. 

Through digital logic control, sub-routines start and termi- 
nate when the corresponding binary control variables change 
state as logical functions of: 

1. External control (switches, relays controlled by external devices). 

2. The states of timers or sub-routine counters. 

3. Analog-comparator decisions. 

The interplay of binary control variables and analog com- 
putation results in a special hybrid analog-digital structure. 
Relays or electronic switches implement analog sub-routine 
changes under control of digital (binary) control variables 
and constitute the digital-to-analog interface of the computer. 
Analog solutions, in turn, can modify digital control. 




• ••••• 

• ••••• 

• ••••• 

• ••••• 



• • 



• • • 



## 



Model 3326A 

Six RST Flip-flops, six AND 
Gates and four relay drivers. 
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Model 3328 

Time-event Control. 1 kHz 

clock with three decade 

dividers, three decade counters 

and one gated relay driver. 



Additionally, in the SD system it is possible to insert a 
delay of the Hold-mode command to the Complementary 
Integrators. This delay is important because it enables the 
C Integrators to store the final values of their inputs in a 
problem solution where the C Integrators are tracking rapidly 
changing problem variables. 

The Reset-Compute-Hold modes are, of course, settable to 
any desired intervals, and it is this new degree of freedom 
imparted to integrators that illustrates the importance of 
digital logic control. 

The combination of the three SD digital logic modules 
(Flip^Flops, Gates, Time/ Event Control) results in a most 
flexible hybrid analog-digital structure. The operator can 
easily set all Reset-Compute-Hold intervals of integrators with 
the Time/ Event Control Module. 

A typical example that illustrates the use of SD's digital 
logic is mode control of iterative integrators. In iterative 
operation (10), results obtained during or at the end of one 
solution of the problem are used to change parameters or the 
circuit configuration (switching) for the next solution. The 
following figure shows how in lO, integrators are paired into 
normal and complementary (opposite) logic to implement 
iterative solutions: 



Normal 



Complementary 



Mode Duty Cycle for Integrators in Iterative Operation 




th( 



* When programming Involves the use of digital control, logic levels for the digital logic elements ore 
opposite of those for the mode control and function relays. 

* Refer to Chapter 8, "Basic Operation of Digital Logic Elements," for complete operating descriptions 
of the Systron-Donner digital logic control modules. 
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Function Generator Operation 



Principle of Operation: The Function 

Generator approximates the curve of an arbitrary function 
with a series of connected straight-line segments as 
shown in Figure lA. This graph represents the sum of 
the individual line segments shown in Figure IB: Each 
line segment corresponds to the output voltage of a diode 
network in the function generator. The point along the X 
axis where each line begins corresponds to the selected 
bias or "breakpoint" voltage. The diode networks are 
arranged so that six conduct when the X input is positive 
(0 to +100 volts) and six when the X input is negative 
(0 to -100 volts). For better accuracy, the input voltage 
may be biased to use all 12 networks above or below 
zero. Two function generator units can be operated as 
a single-channel of 24 segments for more accurate 
simulation. 

The slope of each line segment corresponds to the 
ratio of output voltage per volt-change of input and can 
be varied between the values of +2 and -2 by adjusting 
the current flow through the diode network. For steeper 
changes of line slope at any point, two or more networks 
can be operated in parallel by shifting to the same 
breakpoint. 

When the individual outputs of all diode networks are 
added together by an operational summing amplifier, the 
total output over the entire range of X-inpuf will repro- 
duce the original- straight-line approximation. No output 
is produced when X = 0; therefore, a constant voltage, 
Yo/ must be added at the output summing amplifier, when 
Y is not equal to at X = 0. 

Set Up Procedure: A function may be set up on the 
function generator either by the preparation of a table of 
values to which the individual slope and breakpoint 
potentiometers are adjusted, or by adjusting the potentio- 



meters so that a trace of the output voltage on an XY 
plotter will match a prepared graph as closely as pos- 
sible (visual method). The second method may be faster 
but is limited in accuracy by the read-out instrument. 




A. Straight Line Approximation of a Curve. 




B. Outputs of individuol Diode Bridge Networks, the 
of which produces the Straight Line Approximatior 
shown in Figure 1 A. 



Figure 1. Simulation of a Function with Diode Bridge Net- 
work Outputs 



TABLE 2. 



DATA FOR SETTING UP FUNCTION Y = COSX 



X (DEGREES) 


Y (COSINE) 


X (VOLTS) 


Y (VOLTS) 


CONTROL* 
BREAK SLOPE 


POINT 


-90 





-90 





-6 


-6 


-60 


+0.500 


-60 


+50,0 


-6 -5 


-5 


-45 


+0.707 


-45 


+70,7 


-5 -4 


-4 


-30 


+0.866 


-30 


+86.6 


-4 -3 


-3 


-20 


+0.940 


-20 


+94.0 


-3 ^^^^ -2 


-2 


-10 



+0,985 
+ 1.00 


10 


-I-9R S 


9 7^ . 1 


1 





+ 100 






+ 10 


+0.985 


+ 10 


+98.5 


+ ? Zll^^^^ + l 


+ 1 


+ 20 


+0.940 


+20 


+94.0 


+3 " :» +2 


+2 


+30 


+0.866 


+ 30 


+86,6 


+4 +3 


+ 3 


+45 


+0.707 


+45 


+70.7 


+5 +4 


+4 


+60 


+0.500 


+60 


+50,0 


+6 +5 


+5 


+90 





+90 





+6 


+6 


* Arrows indicate 


sequence of adjustme 


nts. 
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Y (VOLTS) 




BREAKPOINT: 



-6 -S -4 -J -2 11 >2 +3 '■* *S t« 



SCALING: X: 1 VOLT = 1* (INPUT) 

Y: 100 VOLTS = 1.00 (OUTPUT COSINE VALUE) 



Figure 2 Typical Graph prepared for setting up the Function 

Y = COSX. 

Tabular Mef/iot/ (for Model 3352) 

General: The maximum slope change for all segments, 
except the first, Kx, segment, is ±2.5 : 1. The first, 
Kx, segment has a maximum slope change of ±5:1. 
Breakpoints may be stacked to achieve greater slope 
changes. 

Step 1: Set all Break-point/slope change (polarity) 
switches to the correct combination of polarities as de- 
termined by the function to be simulated. (See previous 
discussion on tabulating a function.) 

Step 2: Set all Slope potentiometers to approximately 
center position. Set all Break-point potentiometers to 
one end as determined by the arrows associated with 
the polarity switches. (E.g., if the polarity switch is in 
the +BP, +AS or in the -BP, -AS position the break- 
point potentiometers should be set to the counter-clock- 
wise end.) During this procedure the output amplifier 
may go into overload. Although no harm is being done, 
the computer may be placed in the Pot Set/Bal mode to 
prevent overload. The computer must be returned to 
Reset before proceeding with the adjustments. 
Step 3: Set the Kx/Yo switch as determined by the 
function and by the following instructions: 

a. The Yo potentiometer sets the value of Y when X 
equals zero. This value may be set anywhere within 
the range of plus to minus 100 volts. 

b. The Kx potentiometer sets the slope of the first line 
segment. The Kx slope is adjustable within the range 
of plus to minus five (5). 

c. If the Kx/Yo switch is placed in the Off position 
both Yo and Kx will be zero. 



d. If the Kx/Yo switch is placed in the Kx position Yo 
will be zero, Kx will be adjustable over its range. 

e. If the Kx/Yo switch is placed in the Yo position, 
Kx will be zero, Yo will be adjustable over its range. 

f. If the Kx/Yo switch is placed in the Kx & Yo 
position both Kx and Yo are adjustable over their 
ranges. 

When the switch has been set, the Yo adjustment 
should be made. With X = 0, the potentiometer may be 
adjusted for the proper value of Y^ at the output of the 
Model 3352 (02)- The Null voltmeter or a precision ex- 
ternal meter should be used for this and all succeeding 
voltage adjustments in this procedure. 



Step 4: The Kx adjustment should be mode next. 
Apply the value of X at the first break-point (BP) to the 
input of the Model 3352. (If a potentiometer is used to 
obtain the X value its output should be applied through 
an inverter to the Model 3352. This will prevent errors 
due to the diode function generator loading the potentio- 
meter.) Adjust the Kx potentiometer until the proper value 
of Y is obtained at the output of the Model 3352. 

Step 5: The first break-point (BP) is now adjusted. 
Adjust the first break-point potentiometer until the Kx 
value previously set is offset by approximately 150milli- 
volts. (0.15 volts) This offset compensates for the diode 
characteristics and will help produce a more accurate 
function. 



Step 6: Set the input to the value of X at the second 
break point. Adjust the slope potentiometer associated 
with the first break-point until the output is equal to the 
value of Y at the second break-point. 



Step 7: Adjust the B.P. potentiometer associated with 
the next segment until the value of Y set in the preced- 
ing specification is offset approximately 150 millivolts. 
(See discussion in Step 5). 

The remaining slope and break-point potentiometers 
are set alternately as indicated by the large arrows on 
the panel, by repeating Steps 6 and 7 until all segments 
have been adjusted. 
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CHAPTER 7 
LOGICAL ALGEBRA 



It is necessary to have a rudimentary understanding of 
logical algebra in order to use logic elements. The 
fundamentals of logical algebra are presented for this 
reason. 



LOGICAL VARIABLES 

A logical variable, L, has two values: 

1 (true) 

(false) 



COMPLEMENT 

The cornplement of L is denoted by L. Loosely speak- 
ing, L is the inverse of L. That is, whenever L = 1 
or 0, then L = or 1, respectively. In manipulating 
logical formulas, the two obvious relations below are 
often useful. 

L + L = 1 



COMMUTATIVE LAWS 



ADDITION (OR) 

The logical sum, Lo = L^ + L^, is defined by the fol- 
lowing table: 



Li 


Lo 


Lo 


1 


2 


3 














1 


1 


1 





1 


1 


1 


1 



This is commonly called the OR function and has the 
alternative symbology 



L^U L2 



^1 ^ ^2 = h "- h 



^2 - h ' h 



ASSOCIATIVE LAWS 



h " (4 " L3) = (L, . L^) . L3 



Ll • (L2 • L3) = (L^ • L^) 



DISTRIBUTIVE LAW 



Ll • (L2 - L3) 



COMPLEMENTED SUM 



^1 ■ ^2 ^ -""l ' ^3 



MULTIPLICATION (AND) 



The logical product, L„ = L^ • Lg is defined by the 
following table: 



h 


^2 


^3 














1 





1 








1 


1 


1 



This is commonly called the AND function and has the 
alternative symbology 

L3 = Lj n Lj . 



h* h' h 



This relationship can be seen from the following table: 



Ll 


4 


^1^4 


L^+L2 


Ll 


4 


^1-^2 











1 


1 


1 


1 





1 


1 





1 








1 





1 








1 





1 


1 


1 















COMPLEMENTED PRODUCT 



h = h ' ^2 
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DE MORGAN FORMULAS 





+ 


4 = 

^2 " 








h 


= ^1 




h 


h 


= h 


+ 


\ 




T,. K u^ • ^ u 1 .• xu . In Other words, the exclusive sum is the sum L, + L„ 

These can be obtamed by complementmg the comple- minus the product Li • Lo. A more convenient ex? 

mented sum and product formulas. pression for © can be derived. 



EXCLUSIVE SUM 



7-2 



Lj © L2 = (Li+L2)- (Lj-L2) 



The exclusive sum, ©, of logical variables is defined 12 12' 

by _ _ _ 



h ® 4 = (h ^ V • (Li • L2). ,L2-L2 

This is also called the exclusive OR. L^ © L„ is the = I .T" t -T 

shaded area of the "set" diagram below. 1 2 2 1" 



CHAPTER 8 
BASIC DIGITAL LOGIC ELEMENTS 



Digital logic can be included in an analog computer. 
It is used for making program decisions based on events 
occurring during the analog simulation. It is also used 
for program mode control. The application of these 
elements is discussed in other chapters of this hand- 
book. This chapter describes the basic operation of 
the logic elements. 

OPERATION OF DIGITAL LOGIC IN AN ANALOG 
COMPUTER. 

Logic levels for the elements are 

+28 volts = false (logical zero) 

volts = true (logical one) 

Most of the elements have a logically complemented 
output in addition to the normal output. All of the ele- 
ments except the drivers and the dividers have status 
lights. When the light is on, the element is in the 
"true" state. 

All of the elements, with the exception of the drivers, 
are designed so that their outputs can be connected 
together without damage. This feature is very impor- 
tant because it permits an OR operation by connecting 
two or more outputs to the same point on the patch- 
board. (See driver exception below. ) 

However, when two outputs are connected together in 
this way, they cannot be used as separate logical var- 
iables. To avoid this difficulty diodes can be used. 
Thus, if Li, L2 are two logical outputs and Li + L2, 
Li, L2 are all required for logical inputs elsewhere 
in the program, the circuit below can be used 



r^^' 



L, + L: 



i 



Each of the logic elements will drive any reasonable 
number of others. However, the logic elements are 
not capable of activating anal(^ control relays or fvmc- 
tion relays. For this purpose drivers are supplied. 
The driver acts as both a power buffer and lexical in- 
verter. It should be noted that the logic levels for 



mode control and function relays are the opposite or 
complement of those for digital logic. That is 

volts = false (relay de-energized) 

+28 volts = true (relay energized) 

Thus, the driver not only delivers the current neces- 
sary to operate several relays but also provides the 
logic level inversion required. To use a driver, the 
output of the logic element is connected to the driver 
and the output of the driver to the relay input. These 
connections are made readily at the patchboard. It 
should be noted that the driver is a special element 
and its output cannot be used as an input to the flip- 
flops or DCU's. It can be used as a gate input. 

Important: The driver outputs must never be connected 
to any other digital logic output. Con- 
necting a driver output to ground (logical 1) 
will destroy the output transistor, that is, 
if connected to a gate, it should not be 
OR-ed with any other element. 

The remainder of this section describes the basic oper- 
ation of the digital logic elements . 

DRIVER 

The driver is used as a power buffer and logical in- 
verter from digital logic elements to analog mode con- 
trol, fimction relays, and electronic switches. Its 
logic levels are 



INPUT 


OUTPUT 


(false) 


volts 


1 (true) 


+28 volts 



The program symbol is 



^ 



GATE 

The lexical gates are of the coincidence type. With the 
exception of input number 1, they operate as AND gates 
only on the inputs connected . That is, if no logical 
variable is connected to one of the inputs at the patch- 
board, then this input is essentially de- activated. Note 
that this is the same as connecting a logical 1 to the 
input. 
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Input number 1 is different because it must be connec- 
ted; otherwise the gate will not operate. However, if 
only input 1 is connected, the gate will act as an in- 
verter if its complemented output is used. 

The program symbol for the gate is 



The gate is shown on the patchboard as 




or 




(NORMAL) 
(INVERTED) 



where the input numbers are shown. 

The truth tables below will clarify the operation of the 
gate. 



INPUTS 




NORMAL 


INVERTED 


1 2 


_3. 


OUTPUT 


OUTPUT 


none any any 





1 


INPUTS 




NORMAL 


INVERTED 


1 2 


^ 


OUTPUT 


OUTPUT 


none none 





1 


1 none none 


1 





INPUTS 




NORMAL 


INVERTED 


1 2 


_3_ 


OUTPUT 


OUTPUT 


none 





1 


1 







1 


1 







1 


1 1 




1 





INPUTS 




NORMAL 


INVERTED 


1 2 


3_ 


OUTPUT 


OUTPUT 


none 








1 





1 





1 


1 








1 


1 


1 


1 





INPUTS 




NORMAL 


INVERTED 


1 2 


3 


OUTPUT 


OUTPUT 
















1 







1 










1 


1 







1 










1 


1 







1 1 










1 1 


1 


1 






20- 



J ^_Y,.'LIGHT 



36 o- 



The true state response to inputs 1 and 2 is delayed 
by 50 usee, while input 3 is not. As will be seen in 
Chapter 9, covering logic applications, this permits 
the gates to be used in connection with flip-flops to 
form a shift register. Also, this feature frequently 
obviates the necessity of using one-shots for delays 
in a logic program. 

Each gate has three inputs on the patchboard. If the 
logicalproductof more than three variables is desired, 
the additional elements can be connected in through 
externally patched diodes, as shown below. 




The additional elements can be connected to any of the 
three patchboard inputs, except number 1. 

The gate can also be usedfor pulse-shaping. It is ON 
(logical 1 at the output) when no input voltage is greater 
than 4 volts and OFF when any input is greater than 
6.5 volts. When the net input is between 4 and 6.5 
volts, the output corresponds to the last state of the gate. 

GATE-DRIVER 



This is a special element found in the TIME/EVENT 
CONTROL MODULE. It is a power buffered, logically 
inverted AND gate which performs the combined fvmc- 
tion of an AND gate and driver. The program symbol 
is 




where Li, L2 and L3 and Cq are respectively logical 
input variables and the output voltage. This gate is 
the same in logical operation as the AND gate except 
there is no delay introduced at any input. When Li* 
L2* L3 = 1, the output of the gate is 28 volts; other- 
wise it is zero. 
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The primary use of this gate is for implementing ana- 
log mode control changes based on logic outputs. For 
example, in this way a multiple decade counter can be 
used either to stop the analog computer after a pre- 
scribed number of runs or in conjunction with a clock 
to specify precisely the length of time the analog com- 
puter is in one of several modes . This is accomplished 
by patching the output of a gate- driver to either the C, 
R or H inputs on the patchboard. 

DECADE COUNTING UNIT (DCU) 

The DCU is either a decade divider or a decade counter. 
It contains ten internal binary states which correspond 
to integers from to 9. Note that this is different than 
a binary coxmter which counts in the binary number 
system. There are two types of DCU's in the comput- 
er. Both of these are in the TIME/EVENT CONTROL 
MODULE (also referred to as Coxmter Module). One 
has all the integer binary states available as a "ten 
line" output. This element is called a decade cotmter. 
A status light is associated with each of the integer 
binary states. The other type has only the decade out- 
put available and is called a decade divider. 

There are two master resets for each group of 3 decade 
covinters together with 3 decade dividers. These are 
the OV and +V terminals. The OV input is directly 
coupled to the clock. As long as a logical one is applied 
to this terminal the clock will be stopped. A state 
change from zero to one will cause all the covmters 
and dividers to be reset. As long as a logical one is 
maintained at this terminal, all other reset signals are 
inhibited . 

A state change from logical 1 to at the +V terminal 
will cause all the counters and dividers to be reset. 
(This is called a trailing edge trigger. ) The +V input 
does not inhibit other reset signals, and it does not 
turn off the clock. 

All counter or divider inputs (including reset inputs) 
require a change from a logical to a logical 1 for 
operation. Thus, the application of a logical 1 to an 
input will cause it to operate once. It will not operate 
again imtil the logical one is removed and reapplied. 
This is true also for the master resets for each count- 
er/divider group. 

The program symbol for the decade counter is 




The counter outputs are the ten integer counts from 
to 9. Whenever a logical 1 is applied to the reset ter- 
minal, the counter will return to count 0. Whenever 
a logical 1 is applied to the input terminal, the counter 
will advance to the next higher count. The counter is 
cyclic: when it is in count 9 and receives an input 
pulse, it returns to count 0. 



Note that the counters can be used in series. This is 
accomplished by connecting any output from the first 
counter to the input of the second counter. Then, when- 
ever the transition from alogical to a logical 1 occurs 
at the output of the first counter, this transition from 
to 1 is applied to the input of the second counter. 
This causes the second counter to advance one count. 
Normally, the second counter receives its input from 
the "0" digit of the first. Then, when a transition from 
count to count 1 occurs in the first counter, the out- 
put from count goes from a logical 1 to a logical 
and the second counter will not be advanced in count: 
its input is activated only with alogical change of state 
from zero to one. Note that when the first counter 
goes from count 9 to count 0, then the count output 
goes from to 1, and the second counter is advanced. 



The program symbol for the decade divider is 



DIVIDER 



OUT 



The divider has ten internal integer counts in the same 
way as the counter. However, these are not available 
as logical outputs. An input to the reset terminal will 
cause the divider to return to count internally. The 
output will only return to state if it is in state 1. If the 
logical variable applied to the input terminal changes 
state from to 1, then the divider will count inter- 
nally. The divider is cyclic: count follows count 9. 

With a sequence of input pulses, the output is alogical 
1 except when the internal counter is in count 8 or count 
9. Thus, when the internal transition from count 9 to 
count occurs, there is a transition from a logical 
to 1 at the output which can be used as input to other 
logic elements. However, it should be remembered 
that if the divider is reset durir^ the time when inter- 
nally it is in count 8 or count 9, the transition from 
to 1 will occur at the output then. 



CLOCK 

Each TIME/EVENT CONTROL MODULE provides a 
lOOO cycle clock output at the patchboard. 
Maintaining a logical 1 at the OV master reset termi- 
nal will stop the clock. 

The program symbol for the clock is 




I kc 



The logic^ cycle of the clock is repetitively a logical 
1 for 800 /x sec and logical for 200 n sec. 
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FLIP-FLOP 

A flip-flop (FF) is a one -bit memory. It has a set (S) 
and a reset (R) state. The set state corresponds 
to true (logical 1) and the reset state to false . The 
program symbol is 




There is a status light connected to the S output so 
that when the FF is in the S state the light will be on. 
The FF can be placed in the S state by applying a logi- 
cal 1 to the S input. It can be placed in the R state by 
applying a logical 1 to the R input. Applying a logical 
state change from to 1 to the T input will cause the 
FFto change to the opposite state. The flip-flop R and 
S inputs are direct-coupled. Each can be a-c coupled 
by inserting a 200 pF capacitor in series with the input 
so that these inputs will be activated by only logic state 
changes from to 1 and not by logic levels. 

Internally, the application of a logical 1 to the S input 
turns off the R output, and a lexical 1 at the R input 
turns off the S output. The R and S inputs are level 
sensitive in contrast to the T input which is activated 
by state changes. Consequently, if a logical 1 is ap- 
plied to both the R and S inputs simultaneously, both 
the R and S outputs will be a lexical 0, and a state 
change at the T input will have no effect on the state of 
the FF. If both R and S inputs are a logical 1 and then 
both are changed to a logical simultaneously, the 
final state of the FF will be either R or S with nearly 
equal probability. If the change to does not occur 
simultaneoiisly, the input which sees the logical 1 last 
will control the state of the FF. 



If either the R or S input is energized by a logical 1 and 
if a pulse or logical variable state change is applied to 
the T input during this time, the FF will act like a one- 
shot mono- stable multivibrator. That is, it will make 
a transition to the opposite state for about 5 M sec. 

Two master reset terminals are located at the patch- 
board for each FF module. One will put all the FFs 
in the R state for a transition from to 1. The other 
reset is level sensitive. As long as a logical one is 
applied all the FFs are constrained to be in the R state. 

COMPUTER CONTROL LOGIC OUTPUTS 

The analog control logic outputs in the integrator mod- 
ule must not be used as an input to any logic element 
except a gate. Therefore, special F-RT and R-H out- 
puts are located in the FF and gate modules. These 
can be used as an input to all of the logical elements. 
They are identified as F-RT' and R-H'. The logic 
levels for both F-RT' and R-H' are consistent with the 
digital logic elements. There is a lexical inversion 
between F-RT and F-RT' as well as R-H and R-H'. 
For example, when F-RT' is connected to a flip-flop 
R input, the flip-flop will be in the reset state when 
The F-RT bus is energized. When the F-RT bus is 
de- energized the FF is free to respond to signals at 
its T and S inputs. 



The F-RT and R-H busses are logical complements, 
which is the reason for their selection in logic control. 

LOGIC SWITCH 

There is a manual logic switch in the counter module 
which is convenient for single -stepping programs . The 
output of this Ic^ic switch (ungrounded terminal) can- 
not be connected to any other logic output. The reason 
for this is because the output of the logic switch is 
shimted internally by a large c^acitor. Consequently, 
if this element is connected to another logical output 
it will cause degradation of the rise-time of that ele- 
ment. 
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CHAPTER 9 
CIRCUITS FOR SIMPLE LOGICAL FUNCTIONS 



Several fundamental circuits which are used frequently 
for logical programming are presented below. In all 
cases the trigger is a negative- going voltage. 



This is due to the 50 /usee delay at input 1 of the gate. 
Its behavior for multiple triggering is the same as for 
the 5 /Ltsec one- shot. 



ONE- SHOT 



Variable: 



A one-shot circuit will produce a logical pulse of pre- 
determined length whenever its input is triggered. 
That is, before the trigger, its output will be a logical 
zero. When the trigger occurs its output changes to 
a logical one for a fixed time interval. A one-shot can 
be constructed with logic elements in several ways. 

5 /isec: 

The circuit below will provide a 5 /xsec logical pulse 
when triggered. 



^ 



^ 



If a second trigger occurs during the pulse, then the 
pulse will last for 5 /xsec after the second trigger. In 
other words, after the one- shot is turned on by the 
trigger, it will stay on until 5 /xsec after the last pulse. 
The sensitivity to multiple triggering can be eliminated 
with a gate: 




The logical pulse from this circuit will be delayed 50 
/xsec due to the delayed input of the gate used for L^^. 

50 /xsec: 

The circuit below will provide a 50 ^sec pulse when 
triggered. 



^ 



Ll 



^ 




^ 



X 



R 
>AAA— 1 



vfy 



This one-shot will produce a logical pulse of fixed 
duration. The pulse duration depends on the choice 
of R, C. Its behavior for multiple triggering is some- 
times erratic. The sensitivity to multiple triggering 
can be eliminated by gating. 




This one- shot will produce a logical pulse whose dur- 
ation is determined by both the clock frequency and the 
output line number from the counter. The pulse dura- 
tion for the circuit shown is n times the clock period. 
Obviously several dividers may be inserted in series 
with the clock if desired and several counters may be 
used. This circuit is insensitive to multiple triggering. 

HALF-ADDER 

A half- adder is a circuit for implementing an exclusive 
sum, ® : 




■Lo=L,.L2+L,.L2 



9-1 



Although digital inverters are shown in this circuit, 
they are not usually necessary because all logical ele- 
ments except for the counters have both complemented 
and non- complemented outputs, 

SHIFT REGISTER 

A shift register is a binary register which is shifted 
one bit each time a trigger is applied. The triggered 
half- adder is needed for this circuit shown at right. 

If Li© L2 = 1, then whenever a trigger, T, is ap- 
plied, Lq will be a logical pulse whose length is equal 
to the duration of T. Note that T is applied to inputs. 

If the duration is less than 50 /isec then changes of L^, 
L2 after the application of T will have no influence on 
Lq due to the built-in delay for gate inputs 1 and 2. 




The circuit for a shift register is 




Notice that the inverters for the triggered half- adder 
are not required since complemented inputs are avail- 
able. An initial setting for this register can be entered 
by resetting the flip-flops individually, or with the com- 
mon reset, and subsequently gating a logical one into 
each S input as required. This circuitry is not shown 
because it requires little explanation. Whenever a 
trigger, T, is applied which is less than 50 /xsec (a 5 
^sec one- shot is a convenient source) the following 
state transfers will be made 



Thus, when each trigger occurs, the current value of 
Li is entered into the first flip- flop. The last flip-flop 
can be connected to the first flip-flop, instead of L^ , 
to obtain a register with an end- around shift. 

UP-COUNTER (BINARY) 



The circuit for the binary up-counter is (shown for 
three bits) 



L. -FFl 
FF -FF2 



FF . — FFn 
n-1 



FFI 

T 



v!y 



FF2 



vly 



FF3 



^ 
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The counter is initialized with the common reset. Al- 
ternatively, an initial value can be entered in the same 
manner as for the shift register. This circuit is acti- 
vated by a logically positive change of state. It will 
covint the number of positive changes of state in the 
binary system. Each time FF n-1 returns to the R 
state, the state of FFj^ is changed. The table below 
shows the states of the FF's for successive logically 
positive changes of state of L. 



FF NUMBER 
3 2 1 



1 

10 



Number of 
State Changes 
for L- 


1 
2 



FF NUMBER 
3 2 1 





1 


1 


1 








1 





1 


1 


1 





1 


1 


1 












Number of 
State Changes 
for Li 



DOWN- COUNTER (BINARY) 

The circuit for the binary down-counter is (shown for 
three bits) 



FFI 

r:r\ 



^ 



INITIALIZE' 



FF2 



FF3 



^ 



^ 



2°-bit 



2 '-bit 



2-bit 



INITIAL VALUE 



This circuit is activated in the same way as the up- 
counter. The table below shows how it coimts. (It is 
assumed that the initial value is the binary number 111.) 



FF NUMBER 
3 2 1 


INPUT 
STATE CHANGES 


1 1 







1 1 





1 


1 




2 


1 





3 


1 




4 


1 





5 







6 








7 


1 1 




8 
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TIME INTERVAL MEASUREMENTS USING DIGITAL 
LOGIC 



The digital logic modules, which are part of the Systron- 
Donner Analog Computer complement, can be used to 
obtain measurements of time intervals with a resolution 
of up to one millisecond. The overall timing accuracy 
depends on the number of digits displ«iyed. It can be 
as good as . 01%, which is a significant improvement 
over the typical 3% accuracy obtained through the use 
of an oscilloscope. Also, this particular capability 
of using digital logic modules to make accurate time 
interval measurements can eliminate the need for a 
frequency coimter, which is generally an extra cost 
accessory. 

The basic time interval measurii^ program is given 
in Figure lA, and the timing relationships are shown 
in Figure 2. This program measures the duration of 
a logic signal, L. Initially, all flip-flops and DCUs 
(Decimal Co\mting Units) are in the Reset state. When, 
and only when L chaises from to 1, Flip- Flop 1 is 
toggled into the Set condition. Gate 3 is opened and 
the clock timing signals are sent to the Time Interval 
Indicator. This is a cascade of DCUs of siofficient 
quantity to measure the total time interval. At the end 
of the time_interval being measured, L returns to the 
state. L from the inverted output of Gate 1 resets 
F-Fl. F-Fl, in turn, toggles F-F2 into the Set state 
and turns off Gate 1. The Time Interval DCUs stop 
coimting and hold the elapsed time display. The Set 
output of F- F2 inhibits the output of Gate 2 through 
Diode 1, preventing any further changes in L from 
toggling F-Fl. F-F2*s Set output also opens Gate 4, 
and clock signals are sent into the Display Time DCUs. 

The Display Time Generator is also a cascade of DCUs. 
Outputs from the Display Time Generator are logically 
summed in Gate 5. When the count reaches the value 
corresponding to the desired display time. Gate 5 is 
changed to the 1 state and toggles F-F3. F-F3 is 
connectedas a 5 microsecond one-shot. That is, F-F3 
goes into the Set state for 5 microseconds before 
returning to the Reset state. The Set output of F-F3 
resets all flip-flops and DCUs, and the circuit is ready 
to measure another time interval. The automatic 
Display Time Generator can be dispensed with, and 
the system can be manually reset by using the logic 
switch on the clock module (Model 3328). 

If L changes state before the end of the display time, 
nothing happens due to the Inhibiting action of F- F2 and 
Dl. If the reset action at the end of the display time 
occurs while L is in the 1 state, the system remains 
in Reset. Interval measurement will not commence 
again until L changes from to 1. This is illustrated 
in the Timing Sequence diagram, 3rd interval time. 

To summarize, this program takes no action until L 
changes from to 1. It then measures the current 
time interval, holds the display imtil automatically or 
manually reset, and then waits for the next state change 
from to 1 in L. 

An alternate program is shown in Figure IB for meas- 
uring the time interval between two independent logic 



signals. Lg^ initiates the timing action. Gate 2 acts 
as a buffer so that L^ is not affected by the inhibiting 
action through Dl . L^ terminates the time interval 
when it goes into the 1 state. The interval terminates 
at that time independently of the state of L^. If L^ 
causes trouble by being in the 1 state during part of the 
interval to be timed, it can be connected to a one-shot 
similar to F-F3 and then into the Reset input of F-Fl. 

The program of Figure 3 allows measuring of the Nth 
time interval after reset. It is identical to the program 
of Figure 1 except that a pre-set event counter is in- 
serted between L and the toggle input of F-Fl. The 
event counter can be made of DCUs, or a cascade of 
flip-flops . 

The period of N intervals can also be_measured, by 
inserting the pre-set circuit between L and the Reset 
input of F-Fl. This is useful for measuring equal- 
period intervals with an effective resolution of 1/N 
milliseconds, or for averaging readings over several 
cycles. 

Figure 4 shows how to convert analog signals to logic 
levels. The logic gates can be operated from analog 
inputs; the output will conform to logic states. How- 
ever, the gates have hysteresis. That is, the gate is 
off for inputs greater than about + 8 volts . Once opened, 
the gate will not turn on imtil the input becomes less 
than about +4 volts. Negative inputs to the gates cause 
no action. Also, the gates are not damaged by input 
voltages within the ±100 volt computing range. 

Gate B is turned on at the start of the time interval 
period, thereby toggling F-F A into the Set state. 
Gate B turns off during the interval, but this has no 
effect on F-F A. When Gate B turns on again, F-F A 
resumes its Reset state and the timing is terminated. 

Figure 5 indicates the method used to convert inte- 
grator mode logic bus signals to logic levels. Gate 
input 1 must be grounded or connected to a lexical 1. 
When the mode bus is energized, it turns the gate off. 
L then is taken from the inverted output of the gate 
since this is a logical 1 when the mode logic bus is 
energized. 

Display of the less significant digits of a time interval 
(when the more significant digits are known), can 
^pear on the Model 3328 which allows only a 3-digit 
display. A 3-digit display of the less significant digits 
will suffice when the more significant digits are known. 
For example, to read an 11 second interval to a reso- 
lution of 1 millisecond requires a 5-digit display. If 
the first two digits are already known, the last 3 can 
be indicated on the clock. The coimters are cyclic, 
so overflows of the more significant-digit counters do 
not affect the count in the less significant digits. 

As another example, the period of a 1 radian per second 
sine wave can be measured to 1 millisecond resolution, 
and . 015% accuracy, without displaying the most signif- 
icant digit since the period is known to be 6. + seconds. 
The display would indicate the . 283 part of the 6. 283 
second period. 
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TO FLIP-FLOP 
MASTER RESET 

AND 
DCU OV RESET 
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R 
3T 



1^-^ 
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DCUS 
DISPLAY TIME 
GENERATOR 
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I PRESET TIME PERIOD 



I GENERATOR | 




DCUS 
TIME INTERVAL 
INDICATOR 



Figure lA. Time Interval Measurement Program 





DCUS 

TIME INTERVAL 

INDICATOR 



TO FLIP-FLOP 
MASTER RESET 

AND 
DCU OV RESET 
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DCUS 
DISPLAY TIME 
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PRESET TIME PERIOD 
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Figure IB. Alternate for Interval Between Two Logical Events 
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u 



COUNT 

DISPLAY 

RESET 

COUNT 
RESET 



DISPLAY 
TIME 



INTERVAL 



INTERVAL 
2 



I 

I I 

I I 

I I 



INTERVAL 
3 



I I 



TIME 

INTERVAL 

INDICATOR 



n_?' 



INTERVAL 
4 



ISPLAY 
ME 
GENERATOR 



Figure 2. Timing Sequence 



{> 



DCU'S OR 

FLIP-FLOPS 

EVENT COUNTER 



n 



/g^ 



^ 



^>\— 



I P?|SEJ JO_N I 



Figure 3. Alternate To Measure Nth Interval After Reset 



t^^ 



ANALOG 
SIGNAL, S^ 



T A 



S/^ MIN.< +4 VOLTS 
Sa MAX.> + 8 volts 



Figure 4. Program to Convert Analog Signal to Logic Signal 



(LOGICAL I) 



INTEGRATOR 
MODE LOGIC 
BUS ( + 28 V) 



FL 



■C^L 



Figure 5. Program to Convert Integrator Mode Logic (+28V) Bus to Logic Signal 
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CHAPTER 10 
CIRCUITS FOR SIMPLE LINEAR AND NON-LINEAR FUNCTIONS 



Several simple operations are used frequently for pro- 
gramming. They include division by a constant, dif- 
ferentiation, time delay, and a variety of sectionally 
linear transfer functions. These are described below. 
(It should be remembered that all electronic switches 
in this section must be used as J inputs in series with 
the input resistor. ) 

DIVISION BY A CONSTANT 

The circuit 



e, — vVV^ 



O.IM 
eg — AAA 




The symbol 




DIFFERENTIATION 

Programs are usually written to avoid differentiation. 
This is done because differentiation inherently stresses 
noise and because the circuits used may lead to insta- 
bility. However, in some instances it is difficult to 
avoid forming derivatives. This is especially true for 
on-line data processing or editing applications. For 
this reason, it is helpful to be familiar with available 
differentiation techniques. 

A single amplifier program which can be used is 



-AAA- 



R 
-AAA- 



O 



♦— e© 



which has the transfer function 



— = RCs • -^ = 

e. rCs + 1 



r is usually about lOK producing reduced gain at fre- 
quencies above (2TTrC)"-'^ cps. Since C is generally 
large, the circuit tends to load unduly the preceding 
computer element. Also, a large C may lead to sta- 
bility problems in the amplifier driving the differenti- 
ator. The small resistor, r, helps prevent drive cir- 
cuit instability. 



The more sophisticated program 



satisfies the condition 



ae^ + e^ + lOeg = 0. 



(The input gain labels are usually omitted if they are 
unity. ) The transfer function is 



1 10 



This circuit is useful for dividing a sum of variables 
by the same constant, a . Sometimes it is more con- 
venient to set a pot to a, if a is a natural parameter 
in the program. Otherwise the pot would be put between 
the amplifier output and the following input and woUld 
require the calculation of a~^ for each new setting. 

This configuration is also useful for obtaining gains ol 
greater than 10 without requiring more input resistors 
or more amplifiers. 




has the transfer function 



'o s 



1 + Ts 



e, K,K2 • 



1 + 



m 
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where T = CRin» Hiii' the input resistance at ej, is 
usually 1 Mn. Thus, expressing C in microfarads, 
T(sec.) = C(/xf). 



Ideally, 



= 1 - 



K^KgC 



but in practice 'fx is adjusted to a slightly different value 
due to amplifier and potentiometer phase shift. The 
error of the circuit is strongly dependent on a causing 
its setting to be critical. C must be selected so that 



a is large but not close to unity. This is due to the 
stability of the loop composed of amplifiers 1 and 2, 
and pot 3. A good choice for C is 



C = 



0.3 



K, 



^2 



The circuit will give good performance for 10 < K1K2 
^ 1000. The lower limit is necessary in order to limit 
the size of C and thus avoid stability problems. It fol- 
lows that the setting for a is roughly 0. 7. 

The test circuit for adjusting a is (the block denotes 
the differentiator circuit). 



OSCILLATOR 




DIFFERENTIATOR 
K| K2 



+S 



K.K 



|r\2 



KiKzp V^ 



■ERROR 



The oscillator is set to the highest frequency of interest 
and is adjusted until the error is minimum, using a 
scope. In practice a> 0. 7. 

An excellent low-frequency differentiation circuit is 



This circuit works well for 



K, 



1 < — ^ < 100. If desired, 
^2 




pot 1 can be adjusted precisely with the previous test 
circuit, in which the integrator input is set to 






TIME DELAY (Continuous) 

Sometimes a program for time delay is needed for 
solving a problem. The time- domain definition of time 
delay is 



It has the transfer function 



e^(t) = e^(t-T) 



o 



^ = -K.S f^^-^ 

^1 ^ \. OOls"^ + .01s + Kg 



and at low frequency 
e 






This has the transfer function 



= e 



•sT 



1 



There is no convenient way of generating a continuous 
time delay with only analog computer elements. How- 
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ever, various approximations can be made; most of which are due to Fade. The simplest (first order) is 



e^ 1 - -s-s 

O i 



1 1 + ys 

which works reasonably well for a delay of 4>2 0. 6 where 4> is expressed in radians. Since 4> = wT, where w is the 
frequency, 

rr, 0.6 

0) 

That is, at very low frequencies, the approximation will provide relatively large time delays with reasonable 
accuracy. The program is 



ei -t 



^?t^ 




^> 



In practice, the highest significant frequency in the simulation is estimated in order to determine an upper 
bovind for T. 

A second-order approximation which is reasonably accurate for 4> £ 1. 8 is 

Ts T^s^ 

e. Ts. T^- 

^ ^ ^ 2 12 



The program is 



Q>-<!^^><) 



D> 



<] 



•eo 



A fourth-order approximation which is reasonably accurate for <t) s 7. 5 is 



1 - ^(Ts) + ^(Ts)2 - -A(Ts)3 ^ ^(Ts) 



28 



84^ 



1 1 +i(Ts) + ^(Tsr + oi(Ts)'' + 



28 



84^ 



1680 



(Ts) 
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DEAD- ZONE 



The dead-zone operation relates e-, e by the sectionally linear function 



-b 




slope * -m 



The analog program is 



+100 




e.o " 



-100 



In the following circuits note that mode relays and electronic switches are operated by logical zeroes. The 
circuits have been shown with electronic switches for mode control. Function relays can be used instead of 
electronic switches, if high-speed operation is not required. 
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Pot 1 is set so that its diode will not conduct until e. > a. Pot 2 is set so that its diode will not conduct until 
e. < -b. Thus, for -b < e. < a, neither diode will conduct and e = 0. Also if pot 3 is set so the total gain is m 
then 

e^ = -m (e^ - a), e^> a 



= 



, -b< e. ::£ a 



= -m (e. + b), e^^ ^ -b. 

Note that pot 3 cannot be set to m exactly, since the source impedance due to the pot-diode networks is not zero. 
If digital logic is used, the program is 



ma 
100 

-100 (7) 



mb 
100 



+100 



-<D 



ELECTRONIC 
SWITCHES 




DEAD ZONE 

4'. 



a 



m (slope) 



ei 



-,11 N.O. =Normally Open | ^ ,, ^ 
^^z ^whenU = 



N.C.=Normally Closed 



me^ 


Ll 



1 


L2 

1 


1 


^2 



^3 
1 


ES 


-mb < me^ £ ma 


o 


+, > ma 


1 





1 





1 


1 


c 


-, < -mb 





1 














c 



ES, 



ES. 



O 

c 
o 



o 

o 
c 



■m(e^ - a) 
■m(e. +b) 



Pots 2 and 3 are used to set the breakpoints for the 
fxmctions. Their outputs are used both to set the com- 
parator references and to provide a biased input to the 
amplifier as required. The amplifier output will be 
zero when all three switches are open. Logical vari- 
ables Lp L„ are defined by 

h^ = sgn (e- - a) (output of CI) 

Lg = sgn (e^ + b) (output of C2) 



ES 



is open when L. = L, = 1, which is equivalent to 



-b < e. :£ a. 



Since pot 1 is set to m, eQ will have the desired func- 
tional relationship to e.. Normally closed and nor- 
mally open contacts of the electronic switches are used 
to provide the logic level inversion required. 
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SATURATION 

The saturation function is 



eo 



-m 



a 



The analog program is 



m 



^>-=t> 



+ 100 





i\ e< 



-100 



Pots 1 and 2 are set to limit at -a and b respectively. This circuit gives a "soft" limit. That is e« does not 
remam constant at the saturation level, but continues to increase with increasing ei due to the nin-Sero feedback 
resistance when the diodes are conducting. ^ leeaoacK 

The program with digital logic is 



SATURATION 
FUNCTION 



+ 100 



-100 




^1 




^1 



^1 
1 


1 


1 


^2 

1 





ESj 


ES2 


ES3 


b 

m 1 


a 
m 


C 








J- m 







1 

















c 





-, e- < — 

1 m 




1 





1 





1 


1 








c 



-me. 
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The operation of the electronic switches is governed by 
the logical equations 



The analog program is 



ESj = 


W 


^2 


^m 


" ^ ^ m^ 


ESj = 


4 




<«i 


ni 


ES3 = 


h 




(e^ 


m 



The normally closed or open contacts of the electronic 
switches are used to provide the logic level inversions 
where required. The logic program will provide hard- 
limiting; i.e., the slope after saturation will be zero. 
This is not true for the purely analog program. 

BINARY (bang- bang) 

The binary fvinction is 



— 



CiO 



-100 



+ 100 




6-100 



Pots 1 and 2 are set to limit at -a and b respectively. 
As in the "soft" saturation function, the output in the 
limiting region has a non-zero slope. 

The program with logic is 




-100 



The NC (normally closed) contacts of the switch are used. 
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ABSOLUTE VALUE (variable slope) 

The absolute value function with variable slope is 



eo 



-m, 




m, 



The analog program is 



The program with logic is 



■2'('^2''"''''l' 

■o- 



O 



o 



m, 



[> 





■^(mg+m,) 



BACKLASH 



The backlash function is 
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The analog program is 




^ 



♦f- 



[E> 



<l> 



>-IOO 



Pots is set to ^ where ni is the desired slope. 



Pots 



1 and 2 are set to limit at -b and a respectively. The 
limiting network and amplifier 3 comprise a dead-zone 
circuit whose output is an integrator rather than a sum- 
mer. Integrator drift will be about 50fiv/sec/[j.fd, so 
it is necessary to consider two factors when selecting 
the integrator capacitor: 

1 - Sufficiently small IC time constant to permit 
adequate tracking 



2 - Sufficiently large C to ensure that the inte- 
grator drift is permissible. 

The logic program for this function is obtained by re- 
placing the analog dead- zone circuit with the equivalent 
logic program. 

BISTABLE FLIP-FLOP 

The bistable flip-flop function is 





/ 


•o 




— a 




e 


2 








V 














\ 


/ 




/ 


\ 


A 


/ 


\ ^^^ 


^^ 




-c 


\ 






c 


*l 




^v 


/ 




y^ 






-a 











The analog program is 



DITHER 



The dither function is 




-A 



Pot 1 adjusts the value of e- which results in a change 
of state, namely C. 
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The analog program is 





BANG -BANG 






r^ 










^>i 








-^ 


_J1>^- 






A 
100 ^ 


V 








-MWW^ 1 / 







In the bang-bang circuit the limits are set to A, -A. The logic prc^ram is 



k. 



1 



-100 shA/St 



+100 AAA/ 




100 



In both the analc^ and the logic programs, the gain of integrator 1 determines the dither frequency. 
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CHAPTER 11 
SIMULATION OF CONSTANT COEFFICIENT TRANSFER FUNCTIONS 



The simulation of control systems requires programs 
for constant coefficient transfer functions of the type 



where G(s) is the forward path transfer function and 
H(s) is the feedback path transfer fxmction. 



T(s) 



P(s) 



where P, Q, have polynomial form: 



m 



P(s) =X;ajs', Q(s) £ b.s^ 
i = j = 




Kn > m differentiation is implied. Differentiators are 
treated in Chapter 10. Here, it is assumed that n < m. 
For physically realizable systems the ai, bj are all 
real numbers. It is well known that any polynomial 
with real coefficients can be represented as aproduct 
of quadratic factors. Further, in many practical 
problems the physical model is composed of dominant 
second or first order transfer functions. For these 
reasons only cases for which m < 2 are considered. 
Programs for specific higher order transfer functions 
can be derived easily with the approach outlined in 
Chapter 3. The material here is intended to comprise 
a dictionary for second or first order programs . With 
this information, generating a program for a complex 
system diagram with many blocks is easy. The pro- 
gram is written for each block and the system is rep- 
resented by making appropriate interconnections be- 
tween blocks . 

In developii^ transfer fvmctions using only pots, am- 
plifiers, and integrators, the followii^fimdamental 
characteristics of feedback or closed- loop systems 
are of interest. 

General Closed- Loop System 



The transfer functions of this system are 



^0 ^1 S 



El 1 - HGi Gg 



^ S ± / M 

Eg - l-HGjGg = G^\eJ 



That is, the response to inputs other than Ei is equal 
to the normal response divided by the transfer fimction 
preceding the actual input. 

Specific transfer function programs are given below. 



T(s) = a: 



E,«- 



O 



Eo 



-1^^?) — SUMMING POINT SYMBOL 



-^^g) 1 G(S) I 1 

I H(S) 



The transfer function of this system is 
^0 _ G 



E^ 1 - HG 



T(s) =^: 



T(s) = -4-: 
^ ' 8 + a 

E| 



1^ 



;fl> 



c ^0 " 



a 



■Eo 



Eo.jl_ 
El *s+a 



U-1 



T(«) = siVi 




T(s) = i^ 



E| 



/8 
OH 



^ 



L-O 



•En :r- 



Eo, -)9 



E, S + a 



T(s) = ^^ : 
^ ' s + a 



/3 
E.-0 



;i> 




Eq -iSS 



E, S + a 



T(s) = ^ + " 

^^®^ S + /3 



a>0 E, — f 



a-i8 



o 




^ 




i9 



Eq S-ha 
^0 E|"s+)3 



a< 



y9 



E, — • 








^ 



^0 . S+g 
E| "s+iS 



T(s) = ^ 



S + as + /3 



iS>0 E 
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T(s) 



s + as +/8 



fi'O 



<i> 



ffi> 



■O 



B> 



E| 



Eo 



= T(S) 



(2) 



)9<0 E, 



v^ 



vqs 



1I>T^^^ 



o 



b^ = T{S) 
E| 



(3) 



For a, i3 > 0, a^ > 4^, another program is preferable. 
het6=f, \^ = 6^ -p. Then 

y 
"^^^^ " (s + 6 + x) (s + a -X) 



and the program is 



E, 6-— []/> 



•-o- 



s + x 



^^ 



* Eo 



|S = T(S) 




Note that the condition a, /3 > Censures (S+X)^ (8-X)>0. 



T(s) = -2 



US + V 



s + as + 



The program for this function uses types (1), (2), & 



(3)ofT(s) =—2 



shown above, with the 



s + as +/8 
addition of two pots and one or two summers. In each 

of (1), (2), (3), pot 1 is set with 7=1 
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iS>o r 



ii> 



'-^> 



11' 



(I) 



^ 



r 

o- 



o 



t> 



En iri=T(s) 



/9-0 r 




En -2=T(S) 



/9<o r 



Mv^ 




T(s) = -^ 



S + liS + V 



+ as + 



For convenience T(s) is simplified 



T(s) = 1 + i^^- oc)s ^ {.-^) 

s + as + 6 



and the programs are (where againextra circuitry has 
been added to (1), (2), (3) above and in each pot 1 is set 
with y = 1). 
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E,— f- 






li'l 



c> 






O ^ 



O 



Eo |2.T(S) 



iS.O 




i9<0 




Eo — »T(S) 
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CHAPTER 12 
CONTROL SYSTEM SIMULATION 



A model of a control system is usually represented by 
a block diagram The block diagram contains constant 
coefficient transfer functions, non-linearities, and 
frequently, special differential equations. The easiest 
way to generate the complete computer program for 
the model is to 



1 - derive a program for each block 

2 - interconnect the blocks. 

This method iscalledblockprogrammingandis simple 
using information developed in Chapter 4. Some illus- 
trative examples are shown below. 



MODEL: 



KS 



S+ o 



S+j8 



s + a 



1 

's + 




<3 



O— ' 

/9 



PROGRAM: 






<3 



<v 
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MODEL: 



R +xr-K + 



'dS^^ 



Ks 



s + s+ I 



s + l 
s+3 



vTcM 



7 



1 



Ks 



s + s + 1 



tt> 



<3 



O 



\!> 



[> 



s + 1 
s + 3 



p> a: 



0.2 



O 



;^ 



10 



3> 



•-o- 



0.3 



x. 



\ 



i+ioo 
— * 



•-I00 



{^^ 



12-2 



I I - 



t> 



t> 



> 



PROGRAM: 



ABSOLUTE VALUE 
CIRCUIT 



O— 1 




Note that no pot is required for the parameter, B, be- 
cause this parameter is used to translate the dead- 
zone fvmction. 
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CHAPTER 13 
FUNDAMENTALS OF VECTOR ANALYSIS 



The development of vector operators given below is 
within the framework of Euclidean (x, y, or x, y, z) 
space. Due to the invariance property of vectors and 
vector operators under transformation to other coor- 
dinate systems of engineering interest, it is reasonable 
to proceed in this manner. Sometimes we will show 
the representation of the vector operator under discus- 
sion in other coordinate systems for the purpose of 
illustration. Indeed, the main motivation for using 
vector notation is its non- dependence on a particular 
coordinate system. 

DEFINITIONS 

(A • B) is the scalar (dot) product of two vectors. 

A . B = |a|. |b|Cos (A, B) 

(A X B) is the vector (cross) product of two vectors. 

If 

C = AxB 

then 

|C| = |a|- |b| Sin (A, B). 

Note that 

(C x (A X B)) = 

Vis called the vector operator del, 



V ^^ 



8_ _d_ 

ay' dz 



(Here ordered n- tuples are used instead of the cus- 



tomary i j k unit vector notation, i. e. 

a .. 



d_^ d_ ^ d_. 

ax ' ay ax 



If (^ = ^ (x, y, z) is a scalar- valued function, then 

^ ax ' ay az 
is a vector-valued function also known as grad <t>. 
Vhas different forms in other coordinate systems. 
Example: In cylindrical coordinates 

ar' rae' az 



For a vector V = Vx, Vy, Vz; 



, av av av 

„ -, A z y X 

y z z 



aV aV aV 
z y _ X 

a. ■ 



X X 

is a vector-valued fvmction known as curl V. 



V V = V 



Laplacian. 



2 2 2 

82^82^92 
X y z 



is known as the 



SUMMARY OF OPERATOR ALGEBRA (greek letters 
are scalar-valued function^, capital letters are vec- 
tor-valued fxmctions) 



V- (aV) = a(v- V) + V- (Va) 
Vx(aV) = a(VxV) + (Va) xV 

V • (VxU) = U • (^V) - V • {"TkU) 

V(V • U) = V X (VxU) + (V • V) U + Ux (^V) + 
(U • V) V 

Vx (V X U) = (U • V) V - (V • V) U - U (V • V) + 
V(V- U) 

Vx (Va) = o 

V . (v X V) = o 

Vx (Vx A) = V (v • A) - v2 A (this defines v 
which is not a true vector operation). 

AxB = - BxA 

A • (B X C) = (AxB)' C (triple scalar product) 

DIFFERENTIAL OPERATIONS 

The remainder of what follows will be derived for two 
dimensions but may be readily extended to three. 

Divergence: 

If f (x, y) is a scalar-valued function then the family 
f = c of curves in the plane are level lines for f . K 
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f is thought of as a potential then these curves are 
equipotential lines. 



Outflux 
unit area 



V Jx+h, y + f) - V^(x, y + f) 



V (xH, y+k) - V (x + i', y) 

+ — i sL 



Vf • (dx, dy) = 

where (dx, dy) is a vector element of arc along the 
curve f = c. (This is so because we constrain the 
values of x, y so that f = c). Thus we see that 

Vf _L (dx, dy) , 

and the family of curves obtained from the differential 
tangent field, Vf, are orthogonal to the family of level 

lines and represent a flux (or flow), vf • N = ^ ^^ 

called the directional derivative of f in the N direction 
(N is a unit vector) at the point (x, y). The meaning 

av av 
is obvious. V • V = t— ^ + r— ^ is a scalar- valued 

9 d 

X y 

fimction known as the divergence of the vector-valued 
f\uiction,V. It is commonly written as div V. 

Consider a small rectangle in the x, y-plane 



«(x+h,y+k) 



(x,y) 



In the limit as h, k-^0 we have 



Outflux 



av. 



vmit area dx 



av 

—1 = 
9y 



V- V = div V 



Thus if (^ is a conservative potential field 



V (f) =OorV^<^ =0 



At a point x, y; V <^ 



rt • 1- 1 .source . ., ^ 
implies a local . , m the 
^ sink 



case of incompressible flow. Note that with compres- 
sible steady flow v^ ^ ?= 0, in general, without the 
presence of sources or sinks. In transient physical 
processes v^ is a measure of the change of energy 
from one form to another. For example, in heat trans- 
fer processes heat is stored locally by the distributed 
heat capacity of the medium in which the process occurs. 
Thus, if T is the temperature at (x, y), the local con- 
ductive heat transfer is k VT where k is the thermal 
conductivity, and the local heat stored by the heat ca- 
pacity of the medium is cAT where c is distributed 
heat capacity. The loss of thermal flux (heat flow) 
from the conductive process equals the heat flux re- 
quired for heat storage and this time- varying relation- 
ship is given by 



We think of V as a flux, and determine the net emer- 
gence of flux from the rectangle. 



or 



kVT = c 



a_T 
at 



Vy(x+C,y+k) 



-Vx(x,y+C') 



-■ f{x+h,y+k) 



(x,y)* 



Vx(x+h,y+C) 



■Vy(x + C',y) 



^ ^ k at 



1 91 

K at 



(Diffusion 
equation). 



where K is the thermal diffusivity constant. 

In the case of the small oscillations of a uniform mem- 
brane, if w(x, y) is defined to be deflection, then the 
potential energy stored by deformation of the membrane 
is a function of Vw while the local kinetic energy of 
motion due to the mass of the membrane is a function 

aw 
of ^r- • Thus the time- varying transfer of energy from 

at 

one state to the other is described by 



a w 



T,2 1 

^ ^ ^ -2 
c at 



— K- (wave equation) 



We have 

CXitflow = 



[v^(x.h, 



y + t) - V^(x, y + f 



')] k 



+ rVy(x+^y+k)- Vy(x + r, y)l h 



where c is the wave propagation velocity. 

Curl: 

If we think of the vector-valued fvinction, V, as a ve- 
locity field, we can interpret curl V in a convenient 
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way. Consider the tangential velocity around the sides 
of a small rectangle. 



-Vx(x+C,y+k) 

* • ](x + h,y + k) 




Vy{x + h,y+C) 



Vx(x + c;y) 



The ,^ (spin) around the boundary is given by 



unit area 



Vy(x+h, y + g) - Vy(x, y + r) 



V^(x + |', y) - V^(x + | , y+k) 



In the limit as h, k — *-0 we have 



, .. av av 

circ ulation/ . v y x „ ,, 

unit area ^"P"^) = -^^^ - — = VxV=curlV. 



To visualize the result in three-dimensions, we think 
of an infinitesimal paddle wheel located at the point 
(x, y, z). This paddle wheel is angularly oriented so 
that its shaft has the greatest rate of spin. Using the 



righthand rule, the shaft then points in the direction 
of vxV and its angular velocity is proportional to the 
magnitude of VxV. 

If VxV = 0, V is said to be an irrotational field. 
If V = VxA, V is said to be a solenoidal field. 
If V = V- 0then V is irrotational. (VxV- ^ = 0). 
If V = VxAthenVhas zero divergence (v- VxA = 0). 

VECTOR TRANSFORMATION THEOREMS 

Theor 1: (Divergence or Gauss theorem) 

Let R be a region over which a vector valued function 
V is defined, then (differentiability assumed) 



R aR 



V 



N 



where n is a unit exterior normal to aR (boimdary of R). 

Theor 2: (Stokes' theorem) 

Let R be a 3-dimensional region over which a vector 
valued function V is defined, and let (t be a 2- dimen- 
sional subregion of R, then (differentiability assumed) 



/n 



(VxV) 



da- 



ds 



where N is a imit positive (sense to be determined by aa) 
normal to a and ds is a direction vector which lies in aa. 
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CHAPTER 14 
PARTIAL DIFFERENTIAL EQUATIONS - PART 1 



This section treats the class of partial differential 
equations which can be represented 

V- KV<^ = S<^ <^ = 0(xj, •••, x^, t) 
where the operation, E, is defined by one of 



Helmholtz: 



.001 






r 

"L 



a4> = 

H0 = k 

H0 = k^ 

E*^ - C(X) 



at 



E0 = M(X)^ 

8t^ 



(Laplace equation) 
(Poisson equation) 
(Helmholtz equation) 
(Diffusion equation) 

(Wave equation) 



?i 



The analog computer is used to integrate E^, in other 

words to generate E . This will be discussed before 
considering the generation of V • K v^. The block pro- 
grams for these cases are 

Laplace: 



.001 



- w 

The capacitor is added to prevent algebraic loop insta- 
bility. 

Diffusion: (N is the time scale factor) 




C(x)N 




- B 4> 



Wave: (N is the time scale factor) 



- H * 



The high-gain amplifier forces E^ = when connected 
to the remainder of the program. The capacitor is used 
to reduce high-frequency gain which may result in 
instability. 

Poisson: 



001 




M(x)N' 




- B <t> ' + 100 

The high- gain amplifier ensures E0 - k =0. 



B i> 



With the blocks for E«^ available, all that remains is 
to find how to represent V • Kv<^; the Laplacian of <p. 
The form of V will depend on the coordinate system 
(see Chapter 13). To illustrate the concept we will 
derive the Laplacian for one dimension first and then 
develop it for both two-dimensional Euclidean and Rei- 
mannian manifolds. The extension to three dimensions 
will be obvious . For simplicity it will be assumed that 
K is not a function of time. For the case where K is a 
fvmctionof time, thepotsettoKis replaced by a multi- 
plier. 
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ONE DIMENSION 

VKV* .|j[k(x) |^*(K,t)]. 

The derivatives will be replaced by finite differences: 

Since 

^ f (x t) = lim f(x-^Ax,t) - f(x,t) 
nx,i; iim ^ 

Ax — 
then, for small Ax 



9_ f r^ + M f N - f(x+Ax, t) - f(x, t) 
9x ^^^ 2 ' ^' AX 



Adopting the notations 



1 X 1 + X 

X . 1 n + 1 n 
^^2 = 2 



and 



Using the same approximation 



[v. KV.^], 



[KV0] 1 -[KV<^]^ 1 
X 1 - X 1 



= H</)„ 



Therefore, by substituting the different equations pre- 
viously obtained, we have 



\4 [Vl-*n] Kn4 [^n-Vl] 



X --X 

n+1 n 



X -X . 
n n-1 



/X 1 - X l\ S<i> 
i n+2 "-2J ^« 



The pr<^ram which generates H^ is 



f (x , t) = f for any n 



we have, for finite differences 



Af 



/f2lfllM 



f(x2, t) - f(xj, t) 



Ax 



^2-^1 



For any n 



Ax 



An 







Xn*l-X„.l 



Af 1 
Ax 



f(x -, t) - f(x , t) If X ^- - X *= Ax, for all n, then the program reduces to 



X .1 - X 

n+1 n 



Vi - \ 



Thus 



[KV<^]^ 1 =K(x 1). 



</<X^^jt)-^X^,t) 



^^2 



2 Vl - \ 



Vl ■ "^n 



= K 1 • — 

n+Tf X - - X 

2 n+1 n 



dis- 




similarly 



[KV0] 



<t - <b 
^n ^' 



n-l 



. = K 1 

X. 1 n-^ 

n-j 2 



X - X , 

n n-l 



K_ 1 



K_.l 



p = £ p = £ p = i_ 

^1 AX ' *^2 Ax ' *^3 AX 
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Finally, we have the program for all five types of 
equations above. Two are given below as examples. 

Helmholtz: 




Wave: 



-<^n+l 




ONE-DIMENSIONAL BOUNDARY CONDITIONS 

Boundary conditions of two types are commonly en- 
countered: 



<p(a,t) = f(t) 



[Type 1] 



K(a)V0(a,t) = f(t) [Type 2] 

Sometimes both types are specified simultaneously as 
a mixed boimdary condition. These are approximated 
by 

\ = f(t), x^ = a 



K 1 V<^ ^1 = f(t), x^^l = a 



The programs are shown below for the Wave equation, 
but are typical for all H^. 



Type 1: 
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Type 2: 



TWO-DIMENSIONAL EUCLIDEAN MANIFOLD 




MpN^AX 



A two-dimensional Euclidean manifold is commonly 
called an x, y-coordinate system. Here 

^■^^*-{k'^)- (K(x,y)||,K(x,y)|^) 



and 



^ ax ax ay ay ^ 



As for the one -dimensional case, finite difference 
approximations are made for the derivatives : we will 
adopt the notation 



^, 



n, m 



f n' ^m) 



Mixed boimdary conditions are represented by 




f2(t) 
AX 



This discussion can be extended to higher- dimension 
cases by treating each of the boundary points in the 
same fashion. 



Kl /0i -<^ \k1, m/0 -01 ^ 

n+p, m l^n+1, m ^n, ml n-^ I n, m ^n-1, m 

X 1-x A ^n+l"^n / x 1-x iVWl / 



^■^2 ^"2 



K ^1 

n, m+2 



^-^2 ^"2 



i 



n, m+1 ^n, ml 



y 1-y 1\ ^m+l"^m 
^m+2 m-2^ 



n, m-o /^ -(h ii 

' 2 rn, m ^n, m- 1 

y 1-y l\ym " V-1 J 
m+2 m-gX / 



n, m 



If X i - X = Ax, y 1-y = Ay, for all n, m 
n+1 n ' ■'m+l "'m 

then the equation reduces to 



K 1 

n+g, m 



K 1 

n-rt, m 



2' / \ 2' / \ 

Ax rn+l, m n, ml Ax yn-l, m n, ml 



K ^1 

n, m+5 



K 1 

n, m-K 



Ay 



rn, m+1 n, ml Ay rn, m-1 n, m J 



Ax AyH<^ 



n, m. 
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The program is 




TWO-DIMENSIONAL RIEMANN MANIFOLD 



Note that 



Let the metric for the manifold be given by (ds = arc 
length) 

ds = gij dx dy . 

In what follows we will use summation notation (i. e. , we 
sum on repeated indices). 



g = |gij| = ^ia^^ ^"° ^^^ °" ^^ 



so that G is the cof actor of g. . in the determinant I g. . I . 

1] ij 



g g^ 



ik 



It can be shown, that if ^ is a scalar function, 



V" (f) 



1/9 

Vg" 



VaxV 9x 



Thus, we have a general expression for V in non- 
orthogonal curvilinear coordinates. We could treat the 
non- orthogonal case. However, that treatment is much 
more complicated than for the orthogonal coordinates. 
Since the latter is of greater practical interest, we will 
restrict our attention to 



We define 



iJ _ G^ 



g^j = 0, i ^ j. 



rl3 _ 



- 0, i ^ j 
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Thus, in three dimensions, 



When approximated with finite differences 



ds^ = gjj(dx^)2 + g22^^^^^ ^ SSS^^^)^ 



S - Sll ^22 S33 



— =- I — ll""^ 1 2 3 

8X \ ViJJ 9X / Ax AX Ax 



^)- 



ii 1 / .X 

g = — (no sum on 1) 



We denote 



'11 



(^x^x*, X ); x" = A ,a=i, j,k, A are constant 



Consider the diffusion equation 

V • K V = a<p 

We will need to replace the space derivatives with 
finite differences. 



Ax 



^Ax^ /KV^2^3' 



Ax 



vf7i 



Vl,j,k"\j,k 



i--2' h k 



11 / i+l/2,k,j_ 



Ax^Ax^ /KVg^2^3' 



Ax 



Vs7- 



11 



h-l, j, k'^^i, j, k 



VKV.^ =-^ -^ (K V^g^j^ ) 
Vg" 9x ax" 

= E<^ K=K(x\x^x^) 



and similarly for the other terms. Let us consider 
the one-dimensional projection to determine the sig- 
nificance of these terms. From the definition of the 
metric, clearly 



Now 



and 



■^(KVfgii ^ ) = Vg- E<^ 



ax 



ax^ 



g^^ = 


0, i'^ 


Vgg 


Vg "^^22^33 

^11 viTi 


.r- 22 
Vgg 


^eiig33 


Vgg 


^^11^22 



^ 



ax 



^^ 



KVg^gVg^S a<^ 



y/s^i 






Kv^iVg^a ^ 



VgZ 



22 



ax 



ax 



V ^3 ^/ 



AX Vgj^ 
AX^Vi^ 

AX^Vi^ 

12 3 
are lengths in the x , x , x directions respectively. 

12 3 12 

We denote these by 1 , 1 , 1 respectively. Ax Ax 

3 12 3 

Ax V^ is the volume V = 111. 



In general, the value of the lengths and the volume for 

12 3 
increments in x , x , x depends on location. The one- 
dimensional projection of the problem is 






\j,k **^i, j, k 
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i' ^*V't 



('^'')i-i,M 



<^ DEFINED HERE 



\ i + l,j,k 




(^DEFINED HERE 



X' 
COORDINATE 



(f2^^)i*^J.K 



(CROSS -SECTIONAL AREA) 
<^ DEFINED HERE 



Thus, for example, we see that, in the Diffusion 
equation 

\ iW i-|,j,k \ 1^ / i+l/2,j,k 

are lumped conductivities defined at appropriate loca- 
tions and V. . . S is a lumped capacity, also defined at 

i> ]) k 
an appropriate location. The extension to the non-pro- 
jected problem in three-dimensions is obvious. It is 
interesting to note that the theoretical result agrees 
with our intuition. 

ASYMMETRIC LATTICE 

The set of space points over which the operator 



be found such that the perpendiculars in each polygon 
enclosed by a branch loop intersect at only one point. 




4>.^ 



is defined is called a lattice. In arbitrary regions one 
cannot find a metric in any simple way. Thus, the ap- 
proach developed for the Riemann manifold is npt ap- 
plicable. If a Euclidean system is used, then the bound- 
ary must be composed of a collection of orthogonal line 
segments. This has two disadvantages: the lattice 
density must be high on theboxmdary in order to repre- 
sent the boxmdary. The lattice density on thebovmdary 
determines the lattice density in the interior of the 
region, leading to an exorbitant equipment requirement 
for the program. The asymmetric lattice representa- 
tion of V • KV<^ will allow one to avoid both of these 
difficulties. 

Consider a lattice whose only restriction is that a set 
of perpendiculars to the branches (one for each) can 



Referring to the figure above, 

1. Let a point fimction ^ be defined over the 
lattice 

2. Consider one point of the lattice where we 
arbitrarily choose cjjq and denote by the <^'s 

at branch connected neighbors 

3. Define i. to be the length of the branch between 
lattice points associated with <^q, <I>^. 

4. Define r. to be the length of the _L to i^ 

r. 

5. Define the symmetry ratio Y. = — 
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Now consider the representation of V^<^ at the lattice 
point given by 

a 



Where a is the index set of the neighbors. We will 
choose the Y. and A such as to minimize error. Now 
for the contiguous medium (where we denote the loca- 
tion of latticepoint i by S.)each (^(S.)is given in terms 

of 0(S ) and the derivatives of at S by 



n=l 

We have assigned a Cartesian metric to the domain so 
that i. is the pair (x., y.). Thus the representation 

0{(j) ) is given by 



Define 



Px 'TEVi^' Py4EVl^' exy=27iVl 

a •' cc •'a 



Then can be expressed 



^ MD ^x ^xxo ^y yyo ^xy xyo 



+ higher order terms. 



We want to reduce this to canonical form. The result 
above is independent of an orthogonal transformation 
on (x, y). Note that there exists an orthogonal trans- 
formation such that 6 =0. Henceforth, without loss 

^xy 

of generality, we will assume r = 0. 
Thus, 



or 



Thus 



E-, v^ 1 / 9 a sn 



' n=l 



n! i 9x 



i 9y' 



n=l a 



0{(f>) = E Y.x. <^ + y; Y.y. 



yo 



H:iEY.x.2^^^ -TEVi'^ 



yyo 



+ VY.x.y. d> + higher order terms 
^ 1 ri ^xyo ^ 



0(d) ) = B <P + B <t> + higher order terms. 



\Vq/ V 2 ' ^^xxo ^yyo-^ 



+ ( — jH4(0 - ) + higher order terms. 
^ 2 '^^xxo ^yyo' ^ 



We see that if p is significantly different in magnitude 

from p , then is a poor approximation of v . Thus 

we require the polygon formed by the perpendiculars 
to have relatively good symmetry properties. The 
higher order terms can be reduced by increasing the 
lattice density. One can get an intuitive feeling for the 
symmetry properties of the polygon by observing that 
the system has a physical interpretation in terms of 
moments of inertia as shown below. 



Now we require 

EY.X. = 0, EY.y. = (1) 

a a 

(1), above, is the x and y projection of E Y.i . . 

a 

EY.i. =Er. 

a a 

E^. projected on x, y is zero since the polygon formed. 
a 

by the r is closed. Thus, (1) is automatically satis- 
fied. Evidently our choice of Y. was suitable. 




Y. = Mass at i. 
1 1 
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All that remains is to find the normalizing constant 
A . We choose 



A ^ "" "y 

o 2 



o 4 ^ V 1 •'i ' 2^ 1 1 

a a 



Interpretation: 

If i .. is a branch from <^ to 0., then the lumpedflux 

J ■'•J 

from polygon i to polygon j corresponding to V0 is 



r.. 
_il 
i 



^ U - 0.1 so that Y. . = -^ is a lumped admit- 
ij L J ij 

tance. A. is the area associated with lattice point j. 



o 4 -^ 1 1 

a 

Geometrically, 



MacNeal, R.H., "An Asymmetrical Finite Differ- 
ence Network", Quarterly Appl Math,October 1953) 
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CHAPTER 15 
BASIC ITERATIVE PROGRAMMING 



REPETITIVE OPERATION 

In Repetitive Operation (RO)the entire problem on the 
computer is solved repetitively. This is useful for 
determining the influence of paranieters on the solu- 
tion. Potentiometers and switches are operated man- 
ually while the computer solves the problem repeatedly. 
The RO mode cycle is R, C, R, C, etc. The length 
of the R (initial condition or reset) and C (compute) 



intervals can be set approximately by the controls at 
the console. Usually the C interval is adjusted so that 
the problem solution time is long enoi^h for one to 
observe the significant part of the solution. The R 
interval is set so that all integrators return to correct 
initial values. The RO controls are calibrated only 
roughly. For a more precise setting of the IC inter- 
vals, the program below can be used with the RO 
selector switch at the console set to EXTERNAL: 




<R-Hl 



DIVIDER 



COUNTER 



mi 



1— COUNTER 



m2 



n2 



COUNTER 



m3 



COMP TIME 



n3 




OV 

MASTER 

RESET 



<F-RT' 
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The output from the F-RT' terminal resets the divider, 
counters, and clock, because the logical function 
(GATE 1) + F-RT' is connected to the OV reset ter- 
minal. Note that RH=F-RT'. The clock is off as long 
as the F-RT' terminal is a logical 1. The F-RT' ter- 
minal is a logical 1 when the main machine mode is R. 
The FF is reset whenever 

(GATE 1) + F-RT' = 1 

When the FF is reset all the integrators will be in the 
R (IC) mode since the R output of the FF is connected 
through a driver to the R terminal. Each time the FF 
goes to the S state, the divider and counters will be 
reset momentarily due to the connection of the in- 
verted output of Gate 2 to the +V terminal. The FF 
will be set when Gate 2 = 1. The logical equations 
for the gates are, 



The mode duty cycle for integrators in lO is 



NORMAL^H 




Gl = m, 



m. 



COMPLEMENTARY 



The insertion of the hold interval avoids timing prob- 
lems between mode switches for normal and com- 
plementary integrators . It also has another important 
purpose as shown below. Suppose a complementary 
integrator is tracking a problem variable, x, 



G2 = n^ . n2 • ng. 



Thus, the reset and compute intervals are, 
R = 100m„ + 10m„ + m. 



C = lOOUg + lOng + n^ 




times the period of the divider output. 

Note: The connection to patchboard C terminal turns 
on mode switch compute light during C portion of cycle. 



ITERATIVE OPERATION 

In iterative operation (lO), results obtained during or 
at the end of one solution of the problem are used to 
change parameters or the circuit configuration (switch- 
ing) for the next solution. The term lO has an obvious 
derivation from mathematics. In 10, integrators and 
D AC'S are used for memory. Thus, in particular, there 
is a need for integrators with opposite or complemen- 
tary logic. These are called complementary integra- 
tors. They must be in the R(IC)-mode when the pro- 
blem is being solved and in the Hold- mode while the 
normal integrator is in the R-mode, The C-mode can 
be used instead of the H-mode if the integrator has no 
compute input, but only an initial condition input. 




The chart below shows the relationship between inte- 
grator modes 



Normal 
Integrator 

IC 
H 
C 



Complementary 
Integrator 

C 
H 
IC 



and X is rapidly changing. Then if the C- integrator 
goes to compute when the N- integrator goes to initial 
condition, the C- integrator will not store the final 
value of the N- integrator due to the time constant of 
the IC- circuit: 



N-OUTPUT 
C-OUTPUT 




COMPUTE 



N-MODE 



This problem can be solved if the C-mode command to 
the C- integrator is delayed: 



N-OUTPUT 



C-OUTPUT 




N-MODE 



COMPLEMENTARY 
C-COMMAND 



Without the insertion of the hold interval this cannot 
be accomplished. In some problems, both the N-, and 
C- integrators are used for computation, in sort of a 
bang- bang fashion. For example, if there are two sub- 
simulations, each of which depends on the final values 
of the other, then both types of integrators are used 
for equation solution, not merely memory. In this 
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case delayed commands for both types are necessary. 
If all the variables are sufficiently slow, then the need 
for delayed commands is obviated. The implementation 
of mode duty cycles and delayed commands will be dis- 
cussed presently. 

The programmer needs a switch at the console to start 
10. This is the RO mode switch when the RO selector 
switch is in external. It follows that there must be a 
"safe" mode when not in lO. This is the standby or 
S-mode and is activated by the reset switch for itera- 
tive integrators. Consider the circuit 




It is plain that the N- integrator must be in IC in the 
S-mode. The C- integrator cannot also be in IC in the 
S-mode because the two integrators would form a posi- 
tive feedback system resulting in instability. The C- 
integrator cannot be in compute since it may have an 
initial compute input and consequently would integrate 
to overload. Thus, the C- integrator must be in hold 
during the S-mode. Finally, there may be non-IO 
integrators in the program and these must remain in 



Compute or Hold during the iterative cycle. The fore- 
going can be summarized by 



NON-IO ^ 

INTEGRATOR H 

(NO HOLD) IC 



NON-IO 
INTEGRATOR 
(WITH HOLD) 



N-INTEGRATOR 



C-INTEGRATOR 



MACHINE 
MODE 




t = o 



The patching diagram for 10 integrator control is 



^IRT 

RR^viRO F 
o ^o— o 



RESET 



HOLD 



NON-IO WITH HOLD 



^IRT 

Rr^oro^ 

o— — -^o— — — o 



RESET 



R F 
o o— o 

HOLD 



NON-IO WITHOUT 
HOLD 




RESET 



HOLD 



NORMAL 10 



/iRT 
RESET 




HOLD 



COMPLEMENTARY 10 



The F-RT terminal is energized when the machine 
mode is Reset, putting N- integrators in IC and C-inte- 
gratorsinHold. The F- Hold terminal will be energized 
when the Hold button is pushed, putting both N- integra- 
tors and C- integrators in Hold. There is no conflict 
with the RO C-, R-terminals since these terminals are 
disconnected from the patchboard in any mode except 
RO. When the machine mode is RO, the problem Hold 
or H-terminal on the patchboard will be connected to 



the F-Hold terminal. Consequently, to put both the 
N- and C-integrators in Hold, it is necessary that 
the H-terminal be energized and that neither of the 
RO C-, R-terminals be energized by C, R inputs 
from the patchboard respectively. The diagram above 
shows patching for various integrators with mode 
relays. If electronic mode control is used, then the 
relays are replaced with electronic switches as shown 
in Chapter 6. 
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The logic program for basic lO control is shown below. 




in 



DH1 



HI 



CI 



Rl 



in 



DH2 



H2 



C2 



R2 



FSl 



DOWN A UP 



1^- 



FS2 

► < 

DOWN i UP 



R 

T 
S 



-^DH 



R 

T 
S 




a>' 



R 

T 
S 



fC> 



[> 



[>• 



TO RESET COILS 
OF DELAYED, 
COMPL./ 



> 



> 



+ V 



^ 



HOLD 



> 
^ 



COMR 



RESET 



OV 



A 

F-RT' 



Fig. 1. Iterative Operation Logic Control 



Assume FSl (Function Switch No. 1) is up and FS2 is 
down. When the machine is in the Reset mode, the out- 
put of the F-RT' terminal will be a logical 1 . This out- 
put is used to reset the divider and counters and to 
place both FFl and FF2 in the R- state. The clock is 
off. FFl and FF2 form atwo-bit binary counter. (The 
lO cycle H, C, H, IC corresponds to the counter out- 
put 00, 01, 10,11.) Gate H is ready. When the RO 
mode switch is activated the decade counters will start 
and both N, C- integrators will be in Hold. At lOHg + H^ 



the coincidence at Gate H will reset the decade counters 
and advance the binary counter to 01. (The use of the 
output of Gate DH will be discussed later.) The output 
of the binary counter will ready Gate C, thereby putting 
the N- integrators in Compute, and the C- integrators 
in IC by energizing the C-terminal through its driver. 
At IOC, + C^ the coincidence at Gate C will reset the 

decade counter and advance the binary counter to 10. 
The output of the binary coimter will ready gate 1 and 



15-4 



put all integrators in Hold. At lOHg + H^^ the coinci- 
dence at Gate H will reset the decade counters and 
advance the binary counter to 11. The output of the 
binary counter will ready Gate R, put the C- integrators 
in Compute, and put the N- integrators in Reset by 
energizing the R-terminal through its driver. At 
10R„ + R. the coincidence at Gate R will reset the 

decade counters and advance the binary counter to 00. 
The circuit is now ready to go through the same cycle 
again. 



If FS2isup, then the lO cycle will start with C (i.e. 
will skip the first and only the first H interval). 



it 



Delayed Hold commands are provided for the delayed 
complementary integrators by gate DH and the associ- 
ated FF. The logical equation for gate DH is, 



DH 



H • DH. 



DH, 



Note that 



lODHg + DH^ < 



IOH2-. 



H, 



so that the FF will be toggled before the end of each 
Hold interval. This causes the Reset period of the de- 
layed complementary integrators to bracket the Com- 
pute interval of the normal integrators . 

Delayed Hold can be obtained by other means. Delayed 
Hold can be generated for both N-, and C-integrators 
by delaying the deactivation of the C, R inputs, res- 
pectively. The circuits for obtaining these delays are 
shown in the following diagrams. The inverters are 
those associated with the preceding gates . 





^ 



»• C 



K FSl is down, the logic switch found on the counter 
module can be used to single- step the computer through 
the iterative program. 

If the machine Hold switch is depressed, all integra- 
tors, 10 and non-IO, will go into Hold. This happens 
because the R, C inputs at the patchboard are discon- 
nected from the F-RO and R-RO terminals in the inte- 
grator modules and the main machine mode takes pre- 
cedence. The clock input to the divider is also stopped 
because the RH terminal output is a logical 0. This 
saves the state of the divider, counters, and two-bit 
binary counter. Consequently, the iterative program 
can be put in Hold at any time while preserving its 
status . To return to the iterative program after being 
in Hold, it is only necessary to activate the RO mode 
switch. 

Sometimes it is desirable to control either the C or the 
IC- interval or both by an event in the program. This 
is accomplished by using comparator outputs instead 
of decade counter outputs as inputs to the appropriate 
gates . 



A simpler but less flexible lO logical control program 
is 




frt' > 
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This program will permit only an lO cycle composed 
of 10 time intervals. As before, the RH' output is 
used to gate the clock. 

MODE CONTROL NOTATION 

When delayed commands are to be used for an integra- 
tor it will be denoted by: 



ACCUMULATOR 



r|ND> OR ncD> 




This program will compute the sum of n successive 
values of f : 



depending on its type 

Integrators can be controlled separately by pairs, not 
under general lO control. 

In this case the following notation is used: 



1 " 
, ^ J_ y f (t.) 

'o ^ i-i 1 



The pot is for scaling. The output, e^, as a function 
of time is 



^ 



R = LOGICAL EXPRESSION (IC RELAY) 
H= LOGICAL EXPRESSION (HOLD RELAY) 



For example, R = L^ • Lg. The drivers are still used 

as logical and power buffers. If a resistor rather than 
a capacitor is used for the feedback element of an inte- 
grator amplifier, then the mode control relays can 
still be used. The notation is 



U> OR [€> OR P> 



R = LOGICAL EXPRESSION 
H = LOGICAL EXPRESSION 



Thus, a summer can be made to act as a sampler in a 
sampled data system by controlling the hold relay. 

MEMORY CHAIN 

A memory chain program is 




r 



I I 

J L 



to t, t, 



Note that the complementary integrator has two IC in- 
puts. The second is obtained by patching a resistor 
between the output of the normal integrator and the U 
terminal of the C- integrator. 

AUTOMATIC RESCALING 

Sometimes it is necessary to rescale a problem during 
the computer solution in order to meet the accuracy 
requirement imposed on the simulation. This occurs 
because of the wide range of variables involved. One 
good example of this situation can be found in the sim- 
ulation of a nuclear reactor start-up, where the neu- 
tron level may rise through six decades of power. A 
technique is presented here for automatic rescaling. 
Although only a simple example is used for illustra- 
tion, the technique is applicable to the reactor start- 
up simulation and many other similar problems. 

Suppose it is necessary to generate a function, f (t), 
where 



f(t) = K^e 



en > o 



at 



This circuit will store a sequence of values of f (t). The 
first integrator requires a delayed command to obtain 
the final value of f when the program which generates 
f goes into Hold. The chain can be any length desired. 



K = amplitude scale factor, 
n 



and as the function exponentially increases and reaches 
some limit, f (t) , it is required to perform ampli- 

IXIaX 
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tude rescaling and proceed with the simulation. The 
following program automatically fulfills these require- 
ments . 




COMPUTE INTERVAL 
GATE OF LOGICAL 
CONTROL PROGRAM 



-100 



The integrator required for the generation of f(t) is a 
normal integrator. The compute interval is controlled 
by a comparator. The function switch is in the Up 
position only at the start of the simulation and the ini- 
tial IC of the integrator is entered from potentiometer 
1. The switch could be replaced by a fvinction relay 
energized by the FR bus. The 10 logic program (Fig. 1 
or equivalent) goes through its cycle and waits in the 
Compute mode until an external signal commands it to 
go into the Hold mode. Note that the comparator out- 
put is used instead of the counter output to activate 
Gate C (Fig. 1). This signal occurs when the output 
of the integrator reaches a limit set on potentiometer 
2 and the comparator is triggered. A complementary 
integrator with a delayed command learns this value 
and feeds it to the U terminal of the normal integrator 
through potentiometer 3 and a resistor R. 

As the logic goes to the IC mode, a new initial condi- 
tion is established and the entire process is repeated. 

The amplitude scale factor K is given by the follow- 
ing equation 



We see that time can be used as a dummy integration 
variable. Suppose g(x) is the solution of 



M(d)g(x) = b(t). 



That is 



g(x) = J M ^ (D)b(t)dt. 

Then f |g(x)| can also be generated, using time as a 
dummy variable, from the simultaneous solution of 

.g(x) , 
f(g(x)) = J L ^ (D)c(t)dt 

g(x) =) M ^ (D)b(t)dt. 
The program for the generation of f (g(x) ) is, 



\-h^^ "'"'l. 2, 3- 



where R is in megohms, 6 is the setting of potentio- 
meter 3 and n is the number of scale changes. 

FUNCTIONS OF A DEPENDENT VARIABLE 

This technique is not highly accurate, but is useful if 
precision is not a requirement. Suppose we want to 
generate 

f(x), x = x(t) 
If f is the solution of an ordinary differential equation 

L(D)f(t)= c(t) 
Where L(D) is a differential operator, then 




COMPUTE GATE 



-X 




N>-f(g(Xi)) 



-1. 



f(x) = J L"-^(D)c(t)dt 



The comparator is energized when g (t) = x, so that 
t = g(x). Therefore, the output of the f circuit at this 
time will be 

f(g(x2)). 
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Example: f(x) = Ke"* , -b 4 x 4 b 



FR-Coil 




t-RATE 
100 o Q 



■O— 



\B> 



> 



COMPUTE GATE 




U 



FR 

-o- 



O 

K 



]n> 



K 
100 



-o -100 




The comparator will activate when t = x. The output 
of integrator Nl is Ke" so that the output of integra- 
tor N2 is Ke"^. K these lO integrators compute at 
high speed, say with 0. 001 /ifd capacitors and elect- 
ronic mode control, the output of N2 is a closed stair- 



case approximation to Ke for slowly varying x(t). 
This ID program can be used as a subroutine for a 
non-IO main program. 

The same general comments apply to the following 
examples. 



Example: f(x) = x , x ^ 



t-RATE 



0.6 



-100 




COMPUTE GATE 



k X 3 



2 2 3 3 
Example: f(x) = Aq + Aj x +-27-X +"37^ , x ^ 0. 



100 



-100 



-O 
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PERIODIC FUNCTION GENERATION 

Suppose a non-IO program requires aperiodic func- 
tion. We can use an 10 integrator and a function gen- 
erator to provide the periodic function. The function 
sub-program is 



t-RATE 

-O— 



^ 



DFG 




f (t) 



Here the N- integrator is xmder regular lO control. The 
main program will see the output of the DFG as a per- 
iodic f\mction because its non-IO integrators are in 
Hold during the time the N- integrator is inlC and Hold. 
In fact, the Hold can be eliminated if there is no other 
10 calculation to be performed. The non-IO integra- 
tor must have the logical control 



R=F-RT 
H=F-H 



so that it will be in Hold, not IC, when the lO integrator 
is in reset. 

NEWTON'S INTERPOLATION FORMULA 



f. ^ f(x.) 
1 i' 



N='o 



•- -■ 1 o 



f [Xo,Xj,X2] 
and generally 



f [x^,X2] - f [x^,xj 



^2-^o 



L o 1 nj 

^PrV • •' '^n] -f [Xo,Xp. . ,x^,j 



X - X 

n o 



We have 

f(x) = f [xj + (x-x^)f [x^,x] (1) 

f [x^, x] = f [Xq, X J + (x-Xj)f [x^, Xj, xj (2) 

and generally 



Substituting (2) in (1) 

f(x) = f y + (x-x^)f [x^,xj 
+ (x-x^)(x-Xj)f L,x^,x| 



and 



f(x) = f [x^j . (x-x^)f [x^,xj .••• 

+ (x-x^)»»»(x-x^_j)f [xo"',Xj^J + E(x) 



where 



E(x) = (x-x^)'«»(x-x^)f [^x^*«»,x^,xj . 



If f has polynomial form, then for sufficiently large n, 
E=0. A good approximation for f, in any event, is ob- 
tained by assuming E=0. The result is called Newton's 
interpolation formula. This formula can be used to 
provide a continuous representation of a fimction de- 
fined at n points. The intervals x., x.-l need not be 

equal. However, to simplify our derivation we wiU 
assume that these intervals are equal. It should be 
remembered that the continuous memory to be devel- 
oped can be used for unequal sample intervals with dif- 
ferent values for coefficients in the computer pro- 
gram. We let X =0 and assume all x.-x. - = a. 



Then 



' [V"'^„] 



, n ' on 
n!a 



and 



f(x)=f^.2] 



^ (E-lA 



m=l 



m!a 



m 



m-1 

.n (x-x.) 
i=o 1 



so that 



' o n /x 



V o ^ Z-/ m ! i=o ^a 

m=l 



e-i) (3) 



For a specific number, k, of Epf^ a representation of 
f (^ can be obtained for (3). For n=3, we have 



X <E-l)\x X 



(E-lA 



° ©(?-!)(?- 2). 



6 V ^a ' ^a 
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collecting terms 



f(-) = f + 



(E-l)\ 2(E-lA 



6 J 



(-) 



"(E-l)\ 3(E-l)^f 



L 2 

6 ^a^ 



^]§ 






.(-f3-4f2-5f^-2yfj 



+ (fg - 3f2 . 3f ^ - f^) fj 



CONTINUOUS MEMORY 



In the analog computer we can retain a continuous rep- 
resentation, f, of f by sampling J at n points (not 
necessarily uniformly spaced), f can be generated 
later on command with the interpolation formula. The 
computer program which implements this operation 
from four imiformly sampled values of f is shown 
below. 




Comparators can be used in conjunction with one ramp 
integrator to de-energize the R- relays of amplifiers 1, 
2, 3, 4 in the proper time sequence. The integrators 
are switched from the initial condition mode to the 
compute mode when it is desired to "playback" f(t). 
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K, instead of r(t), it is desired to playback f (T-t), the 
summers and switches are added to the integrator 
chain as shown below. 




The procedure is the same as before except that the 
switches are energized at t=T. 

Although only the circuit for four uniform samples of 
f is shown, the technique can be extended to the case 
of n non-uniform samples. The pot coefficients can be 
calculated previously with a digital computer or desk 
calculator. 
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CHAPTER 16 
SAMPLED-DATA SYSTEM SIMULATION 



Sampled- data systems result from the utilization of 
telemetry, digital control, multiplexed controllers 
and other devices of a similar nature in a control 
system. The analog computer is a useful tool for the 
study, evaluation and synthesis of both linear and non- 
linear sampled-data systems. 

THE SAMPLER 

The single-rate sampler is a switch whose periodic 
duty cycle is shown below. 



At 



CLOSED- 
OPEN 



:i 




tr 



to+T 



to+2T 



-^t 



2ir 
where CO = -=. In order that the sidebands do not 
s T 

overlap, we must have 



Ff = 0, [col ^ ^ 



(The operator F is the Fourier transform. ) 

If there is sideband overlap it is not possible (with 
techniques considered here) to recover f(t) from f*(t). 
Of course, in practice, it is usually not necessary to 
recover f (t) exactly. Consequently, judgment is re- 
quired to determine the extent of allowable sideband 
overlap. This is one question which is conveniently 
answered with the analog computer. 

FUNCTION RECOVERY 

An ideal filter for fvinction recovery has the Fourier 
transform 



It f(t) is sampled by this switch, the output function 
is denoted f*(t). 



T(<^) = A(c<j)e 



i<l>ico) 



A{u}) = 1, Icol^i 



f(t) 




f*(t) 



f*(t) 



DlN 



cir 



= 0, \co\>co^ 
<j){(i)) = -k CO. 
Thus, if F(co) is the Fourier spectrum of f (t), and 
F(co) =0, |co| ^ CO 



In most sampled-data systems. At is small com- 
pared to T. In classical analysis At is assumed to 
approach zero with finite area, so that f*(t) = f (o), 
f(T), — . The introduction of this assumption may 
or may not significantly alter system behavior. When 
the analog computer is used for analysis, the sampler 
can be simulated exactly. The H- relays for the inte- 
grator amplifiers can be used to represent samplers . 
The duty cycle of these relays can be automatically 
implemented with the digital logic. 

The sampling operation produces sidebands. 



then the Fourier transform of the recovered function, 
F'(co), is 



so that 



F'(a;) = T('^)F*(cu) = F(co) e 



f(t) = ^ f (t-k). 



-jkw 



T(co) represents an ideal low-pass filter and a "cardi- 
nal data hold" if k = 0. 



Lf* = i, i; F(s + jno.^) , 
n= -oo 

(L indicates Laplace transform) 



The impuse response of the ideal low-pass filter is 



T(t) = 



o^ Sin [0^^ (k - t)] 
~ co^ (k - t) • 
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The filter is not physically realizable due to the pre- 
cursor (response for t < 0). Ideally ^j^ < w, < ('^g ~ '^i) 

where w ^ 2 w. . w shouldbe sufficiently greater than 
s Is 

2 w so that an adequate practical filter can be found,. 

The foregoing applies to piecewise continuous func- 
tions. Special filter techniques such as clamping 
can be used for discontinuous functions 1. Hold cir- 
cuits are useful approximations of the ideal low-pass 
filter. 



■j^T 



T(jaj) = Te 



Sin (f-) 



T 



= G(a;)e 



j<^(w) 



GH = T 



Sin (if-) 



(^) 



Zero-order hold filter: 

Let the zero-order filter satisfy the relationship 

f(nT+t) = f(nT), O^t < T. 



That is, the hold circuit senses the value of f at 
t = nT and holds this value until a new value is 
sensed at t = nT + T. Suppose f(t) = UQ(t) so that 

F(s) = ^ . 



Then 



F*(s) =J2^ 
n=0 



■nsT 



1 - e 



-sT 



^H= -^ 



CARDINAL HOLD 




A computer circuit which simulates a practical zero- 
order hold filter is shown below. Integrator modes 
are controlled by the lO logic . 




F*(s) 



[^] 



If the transfer function for the zero- order hold cir- 
cuit is T(s), then 



T(s)F*(s) = 



Thus 



T(s) = 



1-e 



-sT 



sF*(s) s 



T(ja;) = 



1-e 



-jwT 



J'^ 



T(j«) = Te 



-J£UT 

2 


jwT -jcoT 
e 2 -e 2 J 




2H^) 



M. Schwartz, "Information Transmission, Modulation, 
and Noise, " McGraw-Hill, 1959. 



If an ideal hold filter is to be simulated, then the main 
computation mvist be interrupted while the memory 
loads. 

First-order hold filter: 

For the first- order hold filter, we have the predictor 

f(nT +t) = f (nT) + ^ [f(nT) - f(nT - T)] , 0^ t < T. 



If f(nT) = 1 for all n, then the output o(t) from 
the filter is 



o(t) = 1 + Y ' 0^t<T 
= 1 , t ^ T. 

Thus, if the transfer function for the filter is T(s), 
then 

0(s) = T(s) F*(s) 

1 



F*(s) 



1-e 



■sT 
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0(s) = (7 + ^ ) (1 - e"^"^) 



Ts 



T(s) = (1 + ^ ) (1 - e-^"^) 



It can be shown that if T(jaj) is expressed 



T(ja;) = G(a;)e 



J^M 



introduction of fictitious samplers or the use of the 
modified Z-transform. The analog computer is an 
effective tool for evaluation of ripple since the con- 
tinuous time response of all elements of the system 
is available for observation. 

PULSE TRANSFER FUNCTION 

In what follows, it will be assumed that all samplers 
are synchronized and have the same period T . In 
order to find the pulse transfer function of an element 
of a sampled-data system, consider the simple sys- 
tem below 



then 



G(co) 



= ^w +. 



Sinf 



t2 



wT 
2 



<^(co) = - ojT + tan " CO 



G(w) 




CARDINAL FILTER 



ci>3 2at2 ^<<^S 

K-th order hold filter: 

The predictor for the k-th order hold filter is 

k 
f (nT + 1) = 2^ ^ a"" f (nT) , o 4 t < 
m = 

where A is an operator such that 



nT 



r(t)- 



We have 



r*{\) 



GCs) 



C(t) 



•C*(t) 



C(s) = G(s) R*(s) 

00 

C*(s) = ^ ^ Cis +jna;g) 



n= -00 



C*(s) = Y ^ G(s+jna;g)R*(s+jna;g) 
n= -°o 

Note that R*(s + jnw ) is periodic with period w . 
s s 

(Consider the sideband structure. ) 
Thus 
R*(s + jno) ) = R*(s), for all n 

1 °° 
C*(s) = R*(s) ^ ^ G(s +jna;g)...Equ. 1 

n = - 00 
C*(s) = G*(s) R*(s) 



Af(nT) = ^ f(nT - T). 



Practical low-pass filters: 

Recovery can also be accomplished with a variety of 
low-pass filters . The analog simulation of the sampled- 
data system usually will provide a means of deciding 
the filtering requirements . Various types of filters 
may be tried in the simulation. 

RIPPLE 

Various responses in the sampled-data system may 
contain components that are either completely or 
partially hidden by the sampling process. These are 
called ripple and may have a bearing on system 
stability. In classical analysis, the system response 
between sampling instants is examined either by the 



From the relationship between the Laplace and Z- 
transforms, it follows 



^^^' ~ R(z) 



Where G(z) is the pulse transfer function and is the 
Z-transform of g(t), the impulse response. For con- 
venience, we will use F(s), F(z) to denote the Laplace 
and Z -transforms of f(t), respectively, although F(s) 
and F(z) have different functional form in general. 

PULSE TRANSFER FUNCTION OF VARIOUS LOOPS 

For convenience, we define 



zlF(s)f ^ Z|L-^F(s)} . 
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The system 



(see equation (1) above). Thus 



r(t)- 



-rfit) 



G,(s) 



Gels) 



C(t)^ 



>(^(t) 



§|] = Z |G^(s)G2(s)( . 

It is convenient to denote Z {Gj(s)G2(s)[ by G^GgCz). 
It is important to note that G^GgOz) ^ Gj(z)G2(z). 



satisfies 



C(s) = Gj(s)G2(s)R*(s) 
C*(s) = [Gj(s)G2(s)R*(s)J * 
C*(s) = [g^(s)G2(s)] *R*(s) 



The computer block programs for the simulation of 
various simple sampled-data systems together with 
the closed- loop pulse transfer function are shown in 
the examples below. In each case the sampling fvinc- 
tion is implemented with mode relays imder lO logic 



s + a 



H = - 

s 



control. In all cases G = 

(The symbol / for an integrator means it is not imder 
lO control. The normal lO logical program is used 
for the sampling amplifiers). 



Example: 



r(t)- 



S&Uiiif;; 



r(t)— ^ 




ret) 

C(z) = GR(z) 



C*(t) 



Example: 
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Example: 



Example: 



r(t) 




•cf(t) 



C(z) 



G(z)R(z) 
l+G(z)H(z) 



r{t) 



rf> 



R=F-RT 
H=F-H 






<1\ 



R=F-RT 
H=F-H 



B> 



«-cf\t) 



r(t)- 



^ 









H 





•C*(t) 



^^^' l+GH(z) 



r (!)■ 



On 



rf> K^ 



<D 



R=F-RT 
H=F-H 



B> 



-Clt) 
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DESIGN TECHNIQUES 



Root- locus: 

Root- locus plotting techniques apply without modifi- 
cation in the z-plane. The region of stability is 
|z| < 1. The interpretation of the locus is somewhat 
different. The family of circles 

Izl = a = aT 



R(s) 




G,(s) 



G(s) 



(where G (s) is the compensating function) is given 
by T(z). Then since 



represent constant damping, a, (not damping ratio) 
for the response 



f* = L® "^ cosojtj * . 

The family of radii 

Arg z = e = wT 

represent constant frequency for the above response. 
Thus, C, the damping ratio is given by 



; = Sin 



-1 



J~2 2 

A/a + (jl) 



or equivalently 



C = Sin 



-1 



xTZ 2 

Vo- + G 



Pole- zero cancellation techniques: 



Suppose the desired pulse transfer function for the 
controller 





T(z) = 


C(z) ^c^^^) 




R(z) 1 + G^G(z) 


we have 








G^G(z) = 


T(z) 
1 - T(z) • 



This technique is useful when G is realizable. 

Digital Compensation: 

A digital controller can be used instead of a com- 
pensation network above. This controller can be 
simulated with the analog computer if it is suffi- 
ciently simple. If not, then a hybrid system may be 
used with the analog computer simulating the remain- 
ing sampled- data system. 

MULTIR ATE- ASYNCHRONOUS SYSTEMS 

Multirate and asynchronous sampled- data systems are 
at best difficult to analyze with classical techniques. 
The analog computer, with its asynchronous control 
structure finds powerful application for such problems . 
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CHAPTER 17 
PARTIAL DIFFERENTIAL EQUATIONS - PART 2 



PROPAGATION OPERATOR 



There is a class of equations which can be solved with 
the similarity transformation 



For (1) above 

u +u = -^r \(-co+lnA). 

VV XX ZTT J 



,(R) = (T'-^a'^T) 0(8R) 



where T is a suitably defined integral transform over 
aR (boundary of the region, R) and propagation into R 

is along the x coordinate of R. That is, the surfaces 

x^- constant are "concentric" to 9R. (T" A^T) is called 
the propagation operator and A the propagation function. 

-1 X 
That T AT exists, can be shown if R is homogeneous 

in the boundary coordinates. That is, usually R can 

be non- homogeneous only in x . We will use this ap- 
proach to find Green's functions for the solution of dif- 
ferential equations where we know the required inte- 
gral transforms over 9R. 

First we derive A for a few simple cases for the pur- 
pose of illustration. 



Consider the cases: 



yy 



yy 



yy 



= -u 



XX 



= u 



XX 



(1) 

(2) 
(3) 



We will use the Fourier transform and its inverse . We 
claim (for propagation along the x coordinate) 



u(x, y) - T"-^A^G(w) 



A^ e^*^ G(a;) do) = 

for arbitrary G. A sufficient condition for the vanish- 
ing of the integral is 

1 2, 2 

In A = a; 



A = e 



and 



u(x, y) 



T-^e-|"I^G(a;) 



u(0,y) = T"-'G(a,) 



so that 



G(co) = (T)u(0,y) 



Consequently 



u(x,y) = (T'^e" "''I'' T) u(0, y) 

and (T~ e '''^ I ^) is the desired propagation operator. 

For (2) we have 

oo jwy 

4- f (-'^ - lnA)e GM dw = 

so that for arbitrary G 



Then, 



= i f A^G(co)eJ"yda; 



,2.x^,, U^ 



yy 



XX 



i f -c^a'^ G(co)eJ"y dc 



f J ln2AA^G(a;)eJ'^dc 



-^ f In A a"" G(a;)e^'^ do; 



A = e 



and the operator is 



For (3) we have 



T-^ e"'^ ^ T . 



2 , 2 , 

-cj = In A 



In A = ± j CO 



i ^ 



■]w 
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Let 



so that 



u(x,y) = T"^ /^ Gj H + t"^ 4 GgCc^) 



,-1 Jwx 



-1 ->x^_2^ 



= T"%-'"^Gj(c.)+T^ e •'"'" G-Co,) 



Then 



u^(x,y) = T'Ve^'"'' G^H - T" Ve'^'^'^GgM 
u(0,y) = T"-^Gj(a;) + T""^ GgM 

u (0,y) = T"^ja;Gj(co) - T""^ jw G2(a;) 



GjH + GgM = (T)u(0,y) 



GiM - s^") =f:3<T>^x^°'y) 



G^(a;) =1 [(T)u(0,y) + J^ (T) u^(0, y)] 
GgCc^) =-|- [(T)u(0,y) - -^(T)u^(0,y)] 

T"^^(T)ju(0,y) 

+ (T"^^^(T)ju(0,y) 

u(x, y) = I T" coscox (T) j u(0, y) 

* (t-1 51^ (T)U^ (0, y) . 



(EI A)- m ^ = 0. 



Let's examine one more case 

^ (EI^)+m 
9y 8y 9x 

Let EI, m be constants and define 



^ EI 



then 



u + k u = 

yyyy xx 



w^ + k In^A = 



2 2 



9 J CO . cj 






^ 



and we see from the last case 



VF 



u(x,y) = /t"^ cos(^=^-^ )(T)j u(0,y) 



2 

. / w Xv 

+ / T ^[ ^ ] (T) I u^(0, y) 



CO 



SUMMARY FOR POLAR COORD'S (r, e) . 



2n 



(T) f (e) = / f(^) e"-'"^ d^ = F(n) 




T-^F(n^ = -^L FMe^""^ = f(e) 



27r 



n=-<» 



For V d> = we have 



u(r,e) = (T"^r"l "l(T))u(O,0) 



2 2 

For V +k<^ =Owe have 



/ 1 J (kr) \ 
u(r,e) = \T-H/^](T)j u(0,e) 



HIGHER DIMENSIONS 



We can derive similar results in higher dimension with 
the use of suitably defined transforms on 9R. 
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GREEN'S FUNCTIONS 

We proceed to use the propagation operator to obtain 
Green's functions. 

Case I: u + u =0. 
XX yy 

u(x,y) = (T-VI''l''T)u(0,y) 

OO 

J u(0,;)e-J"^dt 



i(Ju(0,Udc) f 



e-Hx3Jco(y-Od„ 



Now 



J e-i'^1^ ^ J'^y-^^do, = 



/ eJ^^y-^) e-'^^ dc 



r +COX jw(y-£,) , 
+ J e e"'^-' ' dc 



1 1 

+ 



2x 



so that 



x+j^-y) x-ja-y) x^(y-0^ 



-« X +(y-^) 



Case n: u = u 

yy X 

u(x,y) =^ Je-'^^j'^ydo; 



.( /u(0, Ue-^'^dc 



= 4; /^(O'^)d^ 



/e-'^ x^ja,(y-U^^ 



Now 
00 2 



/ 



^-oj x^ja<y-U 



2 

/e"'^'^cosa;(y-^)daj 



00 2 
+ j /e"'" '^sincoCy-O dw 



2 

/ 





= 2 /e"'^ ^cosa;(y-0 do; 



-(y-U 

1 4x 

e 



2V7rx 



so that 



00 -(y^o 

, \ 1 f 4x u(0, ;) 

u(x,y) = J e ' ' ' 

4irV7r -"» Vx 



d^. 



Case ni: u = u 
yy XX 



u(x, y) = (t' coswxt) u(0, y) 
.(t-1^t) u^(0,y) 

Let f(x, y) = (t"^ coswx TJ u(0, y) 

g(x,y) = (t-1 5iffi^ t) u^(0,y) 



Then 



OO 

f = -=- I coswx e ydco 



27r 



Now 



= i /u(0,Ud^ 

-00 

.//cosc.xeJ'^^y-^^dc.\ 
fcoscjx e^^^'^'doj = f coscox cosa;(y-c) dcu 

_oo — °o 

fcoscox e"''^^y"^^da)= n [b (x-y+^) + 8 (-x-y+£,)J 
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and 



Finally, 



oo 

f(x,y) =Y fu{0,0 [5(x-y+c) 

_ oo 

+ 5(-x-y+^)] d^ 
f (x, y) = y [u(0, y-x) + u(0, y+x)] 



Also 



- . 1 r sinwx ^jwy , 
g(x,y)=-o- J-^j-e "do; 



27r -^ ttJ 



Ju^(0,Oe'^''^dC 



4 fu(0,Odl 7^ [cosoKy-^) 

^TT "^ X ^ ^ I- 



+ ] Sin 



co(y-;)] dc 



for y-;>0: 



C sincjx 



cosoJiy-t,) dw = 7r , y - r, < X 

= 27r , y - C = X 
= 0, y - ^ > X 



for y-^<0: 



I sinojx 



cos(y-£,) daJ= tt, £,-y<x 

-oo 

»= 27r , ^ - y = X 
= 0, ^ - y >x 

for the first case (y-C>0) 

(= TT, y>^>y-x 
for the second case (y-t,<0) 

\= TT, y + x>t,>y 
so in order to meet both conditions 

f= TT, y + x>t,>y-x 



and 



1 y+x 

g(x,y) = ^j u^(0,C)dC 



u(x,y) = I [u(0,y-x) + u(0,y+x)] 

1 y+x 
. i J u (0, Udt . 

•'y-x 



Case IV: 



Consider u = u in the region O^y^Tr, O^x 
yy X ^ 



Tf(y) = f f (^) sin n^d^ = F(n) 



T""^ F(n) =- 2, F(n)sinny 
^ ' TT n=l 



u(x, y) = t" -^ a"^ G(n) = ^ ^h^^^^"" ""^^ ^^""^ 

u =- t l^{-n^) (sin ny) G(n) 

yy 7r n=i 



u = - 2. A^ logA (sin ny) G(n) 
X TT n=i 



We have 



logA = -n2 



A = e 



. 2 
u(x,y) =T''e"" "^ G(n) 



u(0,y) = T"-^G(n) 
G(n) = (T) u (0, y) 



1 2 

u(x,y) = (T"^e""''T)u(0,y) 



, V 2 v^ -n X . 
u(x, y) = - V e sin ny 

n=l 
. / ( u(0, (sin ni) dO 






(0, UdC 



/ V e"" ^ (sin ny) sin nC 
\ n«l i 
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or to write out a few terms 



EXAMPLES 



2 -X. 



■4x . 



u(x, y) = - 1 e'^'^siny f u(0, t,) (sinc)dc + e" ''sin 2y 



/fu(0,C) (sin20dc'\ 



■9x 



+ e sin 3yJ u(0, t,) (sin 3^)d^+ . . 




All of these techniques can be implemented using re- 
petitive operation. The idea is to evaluate at high speed 
the equations for fixed values for the coordinate para- 
meters. Then by changing these parameters slowly 
the solutions can be found over the region, R, of def- 
inition for the equation. 



u + u = 0: 
XX yy 



Let t = C - y 



-"OX +{y-0 



-oo X +t 



X +t 



00 

/ 



^ , ^^jo.y-t) dt 

x+t 



The program is : 
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The FF's act as an up-coiinter with successive trigger 
inputs to FFl. When the machine is in reset the counter 
output is zero and the integrators are in the standby 
mode. (The RO selector switch is in external). When 
the RO pushbutton is activated, the N- integrators go 
to compute. As the counter is incremented, the mode 
sequence is 



COUNl 


'ER 
1 


MODE 


3 2 


FF 








C 





1 


H 


1 





IC 


1 


1 


H 


1 





C 


1 


1 


H 



Note that (FF3) = when the counter has reached 101. 
When this occurs there are no more trigger inputs to 
FFl and the program stops. Before this, the trigger 
for FFl is provided by the clock (L2) when the machine 
is not in compute and by the comparator (Cj) when the 
machine is in compute. There are only two compute 
intervals. The function relay (FR) is energized (+) 
during the first compute interval and de-energized (-) 
during the second compute interval. 



The first and second integrals are calculated during 
the first and second compute intervals, respectively. 
The comparator is energized when y ± t reaches its 
maximum absolute value for each calculation. It then 
terminates the compute interval. The final output from 
the accumulator is the fimction u(x, y) for the partic- 
ular X, y chosen for the run. 



yy 



= u 



2yfiTX 



dt, 



Let t = 4 - y. 



00 — 
u(x, y) = ^— f e^"" u(0, Q d; 



u(x, y) 




y+t)dt 



The an2dog program is 



-100 



^00 




GATE I 



-100 



Both the Ic^ic program and the principle of operation 
of the cong)lete program are the same as for the pre- 
vious example. 
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CHAPTER 18 
CORRELATION ANALYSIS 



PROBABILITY THEORY 

Sample space : 

A sample space, S, is the collection of all outcomes of 
an e:q)eriment. 

Probability: 

P(S) is a function defined over S. P maps S into the 
interval ^ P(S) ^1. Pis not necessarily imivalent. 
The integral (or sum) over the range of P is unity. 
In case S is finite, P is said to ass^n a probability to 
each point of S. In case S is infinite, P is said to 
define a probability density over S. Suppose Sa is a 



nested collection of subsets of S such that S^ C 

for all i. The extremum of Sa= S and the infemumof 
S^is empty. The measure of the rai^e of P for the 
monotone sequence Sa is the probability distribution 
for S with respect to the collection S^ . 



EXPECTATION AND VARIANCE 

The expectation of S, E(S), is the measure of SP(S). 
The variance, V(S); is the measure of 

S^P(S) - E^(S) . 
Example: 
Suppose S is given by - «> < x < « and 



1 
P(x) = (27r) ^ exp 



(x - ^Y 



Then 



00 

E(x) = f X P(x) dx = /i 

-oo 
oo 

V(x) = f x^P(x)dx-/x^ = ( 



Let the S be defined by -<» < x < x , then 



STATIONARY RANDOM PROCESSES 

If the statistics of a process are time- independent, 
the process, is said to be stationary. 

Suppose f (t) is a random function defined for 
-oo < t < «>. 

Let the domain of f be denoted by D. We define an 
ordered countable set JD | so that 



D. : X. :^ t < X. ^ , X. , - X. > A 

11 1+1 ' 1+1 1 



UD. = D 
1 



If f . = f(D.), then the collection f is an ensemble. 
1 r' a 



U 



f is stationary and X is sufficiently large, the ensem- 
ble will beer godic. (That is, the statistical properties 
of the members of the ensemble are identical. ) Hence- 
forth, all ensembles considered will be assumed er- 
godic. 

FIRST AND SECOND PROBABILITY DENSITIES 

Let S be the sample space consisting of all possible 
values for the ensemble )f.(t. + A)l where A< x . 

P(S) is said to be the first probability density. P(S) 
is independent of A, since f is stationary. The ensem- 
ble is ergodic so that 



E(S) = 



Lim 

5— oo 2 6 



J fwMt 



-6 
6 



V(S) 



-^^^ YE J f' (t)dt - E(S) 



-6 



Let S^ be the sample space consisting of all possible 
pairs of values (f.(t. + 2^, f. (t. + A + T)) for the en- 
semble where A + T < X. P_,(S ) is said to be the 
second probability density. P_(S ) depends only on 
T, not on A. The ensemble is ergodic so that 



,^(x^ = j P(x)dx 



E^(S^) - ^^ 



o 

2^ Jf(t)f(t+T)dt 



is the distribution fimction. 



(It can be shown that this equation is valid for T > X .) 
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AUTOCORRELATION FUNCTION 

E (S ), above, is the autocorrelation function for f 
and is denoted by cp^^iT). That is, 



%(T) = l^^ 27 \ m(t . T)dt 

-6 



The autocorrelation function has the properties: 

1) <^jj(T) = <^jj (-T) 

2) (t)r. (o) is the average power of f. 
3) 



%(T) 



^ % (o) 



Proof: 



6 

i, 

6 



kaWh,":" fS j l'«l ■ l'(t*T)|dt 



Lim 

6-*oo 2 



-6 
f(t+T)|2- [|f(t)| - |f(t-fT)|]2Jdt 



. ^ / V Lim 
^f^ '^ 6— «> 



^ J [l^^*)| 



lf(t+T)|] 



dt 



4) f € P »- 0-j e P^ (is periodic with period v) 

Proof: f(t + T) = f(t + T +1^) 



5) <^,o contains no phase information about f . 

6) For every f there exists a imique (j)^; the con- 
verse is not true. 

7) If f & P for any v, then Litti 0« ( T ) » 



T^OOP 



8) If f can be represented by 



f = Y^ A Sin(w t + <t> ) 
/ ^ n ^ n ^n' 



n = -00 



then 



Pft (T) = X) <^nn 



u 

1 f 2 

where <p = Lim «r \ A Sin (^ t + * ) 
^nn ^ 26 J n n n 

-o 
• (sin (co^t + T + (^^) dt j 



Proof: 



^^^ i^ (f \^^^"n* " "^n^ 



-6 



A Sin(a) t + T + 6 )dt = 0, m^^n 
m m ^m' 



9) If f, f ' e bd, then 0,,., (T) = - <j),, (T) 



'f'f 



Proof: 



ci> (T) = Lim ^ t f'(t)f'(t + T)dt 



-6 



= Lim ^ [f(t)f'(t+T)] 



26 

5-* 00 

6 

3fj(t)f"(t + T)dt 



—2 Lim 2^ 
dT 6— «> 



r fj(t)f2(t+T)dt 



CROSS-CORRELATION FUNCTION 



Suppose f, g are two stationary random functions, f, 
g can be decomposed (as above) into two ergodic en- 
sembles f , g . The union of these ensembles is a 

or °a 

composite ensemble of pairs (f, g)^. Let S^ be the 

sample space consisting of all possible pairs of values 
(f.(t. + a), g. (t. + A + T)) for the composite ensemble, 

where A + T < A (a. is the ensemble duration). Since 
the ensemble is ergodic 



o 
E (S ) = Lim ^ I f (t)g(t + T)dt 



6- 



-6 



n = - 00 



E_, is the cross-correlation fxmction for the pair (f , g) 
and is denoted by <^f (T). 
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It can be shown that <^ (T) is valid for T > \. The 
cross-correlation function has the properties: 



1) Generally ^^^(T) ^ <^fg(-T) 



Consider g(t) = 



^ t < 



= 0, |t| > € 

2) 0fg(T) = cp^i-T) 

3) Max <p. (T) does not necessarily occur for T=0. 



4) Uf (T)| ^ V%^oKg(o) 



Proof: 



<^ (T) ^ Lim 25 



^ ( f^t)g2(t+T)dt ^ 



6 

JL f 

26 J 



Lim -^ \ f^(t)dt 

'-^^ -'6 



oo 

^^^ h ] 

5-»oo •» 



Lim ^ \ g (t+T) dt 

6-* 00 



AUTOCORRELATION FOR A SUM OF FUNCTIONS 



Henceforth we will denote 



Lim 

6-*°° 



6 
2^ J 



f(t)dt 



f(t)=fj(t)+f2(t) 



<^jj(T) =Mean |[f^(t) +f2(t)] [fj(t + T) +f2(t + T)]| 
= M|fj(t)f^(t + T)} + M |f^(t)f2(t + T)J 

+ M{f2(t)f^(t+T)} +M{f2(t)f2(t+T)| 

Thus, 



n iji^ ijig I2I1 12^2 

The general rule for a linear combination can be de- 
duced easily. 



EXAMPLES OF AUTOCORRELATION FUNCTIONS 
Example 1 (sinewave) : 

Let f = ASin ( cot + ip), then 

<^,, = Mean { A Sin( wt + ^) Sin( cot + wT + ^) | 



Since f e P^ , 



coA 



a +-7T 



m 2ir 



f 



Sin ( cot + 1(>) Sin( cot + ojt + 1^) dt 



„ aw + ^ + 27r 

m 27r 

dco + jp 



J 



Sin(^) Sin (t, + coT) dC 



2ir 



Pff 



47r J 



coscoT dt= -g- cos coT 



Example 2 (general random function) : 

Suppose f is defined by 

a) f(n6 +t) = f(n(5), o ^ t < 6 

b) P [f(n6) = x] =P(x) 

c) E [f(n5)] = o 



byMean|f(t)} , or just MJf} wherever convenient. Let ^j^^^ 



00 
E [mm + T)] = ( [xP(x)] fxd 

-00 

00 

W [xP(x)] [yP(y)] 



dx + 



dxdy, |T < (5 



00 
= U [xP(x)] [yP(y)J dxdy, |t| 



> 6 



00 

9 



But W xyP(x)P(y) = o, so that 
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% = (1- 



^1 ) ( x2p(x)dx = V(x) (1- 



Example 3 (binary noise) : 
),T-6 Binary noise is a special case of (2), above, where 



o, T > 6 



P(x)=-2 u^ (x-x^) +1 u^ (x+x^) . 




-s 



That is f takes on values of x , -x with probabilities 
of|. 



Then V(x) = X and 
o 



% ~ ^o ^^' 




Example 4 (periodic binary noise) : 

Suppose f is defined as in (3), above, except that 
Hi) = i{t + ^) for large p. Then (p.. is periodic. 





-8 



POWER SPECTRAL DENSITY 

The operator F denotes the Fourier transform, so that 
Fg*^G(a)). 



* M = Lim oT 



r X 

J 

L-X 



1 2 



g(C)e-^'^^d; 



* (w) = F<^ =Lim I e'^'^'^dT 



6-* oo 



-6 



Lim 2^ j g(t)g(t + T)dt 



g(t)dt 



(-) = Lim ±\{ 



Lim 



dT 



* i^) is said to be the power spectral density. We 
have 



<^ (T) = f * (a;)e^'^'^dT 
gg J gg 



WHITE NOISE 

White noise has constant power spectral density of ar- 
bitrary magnitude (say A). Thus, the autocorrelation 
fimction for white noise is 



</) = Au.(T) 



gg 
where u- is the imit impulse. 



* C^) = Lim 



1A 
J 



g(t) e^'*'* dt 



Lim I g(C)e"^'^dC 
6-*°° J^ 



PERIODIC WHITE NOISE 

If f satisfies the requirements for white noise for some 
interval 6, where 6 is large, and if f(t+6-) =f(t)rthen 



/>ff(T) = A^ Uj (T-n6). 



n=-<» 
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FINITE CORRELATION FUNCTION 



If (1) is to be implemented directly, the program 
is 



If f, g are defined only for o ^t < <», then 

1 ^ 
cf> = Lim f Jf(t)g(t+T)dt. 



% 



5-* 00 



(f) can be approximated by a sum: 

N 

*, ^ Mean 2: f (m At) g(m At + T) 
^g m=i 



(1) 



In this case At must be small enough that the high 
frequency components off, g are adequately represent- 
ed (see Chapter 16). N must be large enough that the 
mean of the sum is a good approximation to 0^ . 



CORRELATION PROGRAMS 

Simple computer programs for (1) require that 

At > T . (2) 

To obtain the best approximation. At = T + g, where 
e is the minimum time required for sampling in the 
computer program. (2) imposes a limitation on the 
magnitude of T. That is if T is large, then so is At. 
Thus, for large T, relatively few samples are made 
off, g. This restriction is not imduly severe, because 
(p, cannot be determined accurately for large T to- 
gether with finite histories for f, g. 




Integrator 1 samples f and then after an interval, T, 
integrator 2 samples g. Integrator 3 is in hold vmtil 
both samples are made and then it samples the product 
f(t.) g(t. + T). The output of integrator 3 is averaged 

by aMean Program (Chapter 3) which in turn generates 
(j), . If T is to be specified accurately, the integrator 

capacitors must be as small as possible to minimize 
the effect of the IC time constant. 

The total sampling time can be reduced by replacing 
integrator 1 with a quantizer whose input is f together 
with dither (see reference), f will be represented by 
a two-bit word (output of the quantizer). 

The table below defines the quantizing function 





BINARY 


DECIMAL 


OUTPUT 


OUTPUT 


OUTPUT 


2^x<-l 


00 


-2 


l^x < 


01 


-1 


O^x < 1 


10 


+1 


l^x ^ 2 


11 


+2 
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Except for scaling the program is 



+100 




(£> 



DIVIDER 



(P) 



COUNTER 



ES MEANS 

ELECTRONIC SWITCH I— J 



T| 



I| 



F-RT'> " +V 



COUNTER 




T2|l2 T 



I— COUNTER 



T^T£ 



^ 



R-H' 
V 



T fU 




The Dither Program is found in Chapter 10 and the Mean Program in Chapter 3. 



The correlation parameter, T, and the reset interval 
for integrator 1 are given by (IOT3 + lOTg + T^) P 

and (10I„ + lOIg + IjP, respectively, where P is 

the period of the output of the divider. 



Reference: Korn, G. A. , "Hybrid-Computer Tech- 
niques for Measuring Statistics from 
Quantized Data", Simulation, Aprill965. 
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CHAPTER 19 
ON-LINE DATA ANALYSIS PROGRAMS 



An analog computer can be used to advantage for on- 
line data analysis. Programs are given in this section 
which correspond to the more frequent operations on 
incoming data. The analog computer is used for both 
processing (analysis) and editing. The principles de- 
veloped in other chapters can be used for processing. 
The programs here are for simple editing. 



ERROR FREQUENCY 

Suppose an error is defined to occur whenever 



f(t) - g(t)| > € 



The program which measures error frequency is 



-:=c> 



ABSOLUTE 
VALUE 



r 



-€ 





A 




'- R 

— T 
S 


1 


COUNTER 




V ^ lo 




COUNTER 




-Jo 




COUNTER 



ERROR PEAK VALUE 

Suppose, with the above definition of error, it is re- 
quired to find a sequence of peak values of error, 



Max 
Max 



f - 
f - 



The program is 




+ 100 
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All the programs indicated will be foxind in various 
chapters of this manual. The integrators in this pro- 
gram are under iterative control from the C, R logical 
control inputs . Pot 1 is set to a small value to provide 
positive zero derivative detection. 

TIME BETWEEN EVENTS 

Suppose two events, E- , Eg are defined by 



E^: f(t) - g(t) = e. 



f(t) - g(t) s 



Eg: h(t) - m(t) = eg 
h(t) - m(t) s - 



A program which will determine a sequence of time 
lapses between E., Eg (not Eg, E^) is 



/^^ 



-100 



-100 




-m 



<R-H' 
0V < < F-RT' 
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CHAPTER 20 
SYSTEM OPTIMIZATION 



System optimization is a procedure for generating a 
response, yg, of a known experimental system so that 
it matches, as closely as possible, some desired res- 
ponse y^j. yd may come from tape or may be the res- 
ponse of a known but more complex system. In what 
follows we will assume that either an external device 
generates y^ or else we know the computer program 
for it. A measure of the accuracy of fit of ys needs to 
be established. This is called the Criterion Fimction 
(C F) . The C F has different forms depending on whether 
dynamic or static optimization is required. The exact 
form of the CF in either case is a matter of technical 
judgment and will depend on the particular situation. 
The CF for the dynamic problem will need to take ac- 
coimt of the transient behavior of both yg and y^j over 
some interval of time. Typical choices for the CF are 

STATIC 



F = I (ys - yd)' 



DYNAMIC 



f = C [y^W - y,(t)] 



dt. 



(Note that the CF is defined so that it is non- negative) 

The optimization problem amounts to finding maxima 
orminin/aoftheCF. Since yg is a function of its para- 
meters, say Cj, . . . , Cji, so is the CF. Thus the CF 
can be viewed as a surface over the parameter plane. 
In general the C F will have several maxima and minima. 
We will restrict our attention to finding the nearest min- 
imum from some initial point in the parameter plane. 
Once techniques are developed for the minimum prob- 
lem, they can be applied with modifications, which will 
be obvious, to the maximum problem. 

We will develop two methods for finding a minimum: 
steepest descent and relaxation. The static and dyna- 
mic optimization problems will be treated separately. 

STATIC PROBLEM 

The family of curves f = constant on the surface for f, 
the CF, are constant elevation contour lines. Differ- 
entiating 



df 



n 

j=i 






dC. = Vf 
] 



dC = 



We see that vf is perpendicular to dC and dC is a dif- 
ferential vector which is tangential to the space curve 
f = constant. The projection of the curve on the para- 
meter plane is identical geometrically to the space 
curve because it is at constant elevation over the plane. 
The projection of vf on the parameter plane is perpen- 
dicular to the curve. It is evident that Vf corresponds 
to the s lope of the surface f and vf = at the minimum . 
To get to the minimum from some point on the surface, 
then, all we need to do is follow - vf. (If we wanted to 
find the maximum we would follow + vf . ) The pro- 
jection of the path on the parameter plane is 



STARTING 
POINT 




CONSTANT (TYR) 



If the system is static (algebraic) we can change the 
parameters of yg while foUowii^ - Vf without intro- 
ducing transients in the system. Now if 



^ = T (ys - yd) 



then 



Vf = (y^ -yd)^ys- 



If we let e = yg - y^, then C can be determined by 
(since we follow - vf) 



C = -e 



jvy^dt. 



or in scalar form 



f^^s 
^1 =.-^ia^^*- 



where C is the parameter vector — i. e. , 

C = [Cj, C2, • . •, Cn] (Remember that the CF is a 

fimction of the parameters of yg . ) 



= -.e J. 



dC 



dt. 
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In practice, it is usually advisable to use gain factors, 
K, , so that 



'1 = -Vi 



ac. 



dt 



C = -K e 

n n 



— ^ 

Jac 



dt 



K is chosen so that the procedure converges but at the 
same time is stabje. When vf = 0, the integrals will 
be stationary and C will be the parameter vector (point 
on the parameter plane) corresponding to the minimum 
of f . Vyg is generated by differentiating the generating 
equation for y with respect to the parameters C . . 

DYNAMIC PROBLEM 

Except under special circumstances, vf cannot be gen- 
erated on an instantaneous basis. Even if it could, the 
parameters could not be changed during the determina- 
tion of f. Since the calculation of f is carried out over 
a finite time interval, varying the parameters would 
lead to a result which is not the CF. 

Some definitions for f will permit Vfto be computed as 
often as is f . For example, if 



' - T s! (^s - ^d)'^' 



then 



vf = 



Jo ^ 



y^) vy^dt. 



As soon as we know the components of yg, we can use 
a steepest- descent method. After each solution a new 
position in the parameter plane is chosen by taking a 
finite step along -vf . Now Vyg can be generated by the 
parameter- influence equation method. The differential 
equation for each component of Vyg is obtained by dif- 
ferentiation of the generating equation for yg. For 
example, if 



9C, 



^2 ay 

^1 o ^dC„ 
dt"^ ^ 



) -*- c 



IS 

2 ac. 



-^s' 



dy. 



9yc 



A r "s(o) . _ "^s (o) 
dt ^ aCg ^ ' aCg 



0. 



These parameter- influence equations can be solved 
with an analog program. Their non- homogeneous terms 
are obtained from outputs of the generating program 
for yg . It will be noticed that they have identical form 
except for the non- homogeneous terms. Due to equip- 
ment limitations, it is usually necessary to solve them 
one at a time, switching in the correct non- homogeneous 
inputs from the y^ program in a suitable way. This 
leads to the use ol the relaxation method. 

With the relaxation method we change only one para- 
meter, Cj, until the component, (vf)j, corresponding 
to the parameter is zero or nearly so. The pro- 
cedure is applied to all the parameters on a cyclic 
basis. An example of the resulting path in a two 
dimensional parameter plane is 




* Ci 



This method obviously is ideal for use with the para- 
meter- influence equations . Nothii^ needs to be changed 
in the influence equation program except the non-homo- 
geneous term. 



Ci y. + C„y^ = 0, y (o) = 0, y (o) = A, 



'1 's 



'2 ^s 



<3yg ayg 

the equations for ttt", ^^ can be obtained by dif- 
ferentiating first with respect to C^ and then C,. 



ac. 



d^ ^yg ^^s 

•5— I — ^) + C — ^ 
dt2 ^^^1 2 ac^ 



dt ^ 9C, 



■^s' 



9y^(i) 



In the event it is not convenient to use the parameter- 
influence method, successive values of f can be cal- 
culated for a sequence of values of Cj. This is con- 
tinued until two values of f bracket (vf)j = 0. Then 
the procedure is applied varying C^^j and so on. This 
will also lead to the minimum for f. In either case the 
subscript indexing and switching of non- homogeneous 
terms is done with a digital logic program. In the 
latter case, successive values of f are retained by a 
two or three word memory chain and comparisons made 
by comparators. 

With either method the step size must be determined 
by a gain factor which ensures proper convergence. 
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CHAPTER 21 
MEDICAL APPLICATIONS 

Originally published by: Beckman Instruments, Inc. 
Author: Hiroshi Haxa 



ANALOG SIMULATION OF VENTRICULAR PUMP 
ACTION 

Due to the complexity and non-linearity involved, the 
building of an analog model of blood circulation is one 
of the challenging problems in biomedical science. 
This example is by no means intended to cover the 
subject, but illustrates how iterative programming is 
conveniently used in simulating ventricular pump 
action. 



The ventricle has an inlet valve (tricuspid valve) and 
an outlet valve (pulmonary valve). A cardiac cycle 
consists of two periods called the diastole and the sys- 
tole, as shown in Figure 1. During the diastole, the 
heart muscle relaxes and blood flows into the ventri- 
cle as the inlet valve opens. At the completion of the 
diastole, the systole begins. During this period the 
heart muscle contracts and as the ventricular pressure 
exceeds the arterial pressure, the outlet valve opens, 
expelling blood from the ventricle. An electric analog 
of the ventricle is shown in Figure 2. The analogous 
quantities are indicated below. 



Blood Pressure 



Blood Flow 



Blood Volume 



nth CYCLE 



DIASTOLE SYSTOLE 



Voltage 



Current 



Charge 



(n+l) th CYCLE 



DIASTOLE SYSTOLE 

m> 



FIGURE 1 



The contraction and relaxation of the muscular wall of 
the ventricle is simulated by decrease and increase, 
respectively, of capacitance Cj^ in some prescribed 

manner. If we let E = ^ = elastance of the ventricle, 

the pump action is assumed to be described by the 
following equations: 



[Ed]n = 

hi- 



-t 

(A-E^_l)e%B 



[^djn 



n n ' 



* = f(q5 
n n' 



(1) 



(2) 



where 

Ej L = Elastance during nth diastole 

[e 1 = Elastance during n th systole 

-J 

E = Elastance at the end of n th diastole 
n 

E^ 1 = Elastance at the end of (n-l)th systole 
n- i 

Fq J jj = Blood volume during n th diastole 

* = Slope of [e I ^ to be determined as a 

function of q^^ , the blood volume at the 
end of n th diastole. 

A, B, and t are positive constants. 

Since the static work done by the heart ata particular 
value of q^ is known experimentally (Starling's Curve), 
* must be found as a function of q^^ such that experi- 
mental results are duplicated. The work done during 
the systole is 



W. 



Tsystole 



where 



V2 = arterial pressure 

1^2 = blood flow out of ventricle 

k = constant 

and consequently for a certain q^^ , "^^ must be such 
that FWgl ^= [Wgl Starling. In other words, itera- 
tion with respect to ^^^ is required. 

The passive circuit shown in Figure 2, equations (1) 
and (2), and an iteration circuit are implemented as 
shown in Figure 3. Note that two iterative cycles are 
required to simulate one cardiac cycle. The distinc- 
tion between the diastole and systole is made by a func- 
tion relay (Relay KO) which in turn is controlled by 
one of the memory integrators in a chain of four. The 
action of this circuitry is shown by the timing diagram 
in Figure 4. Since the compute intervals for the dias- 
tole and systole differ, the comparator which controls 
the compute interval of the iterative control program 
must, necessarily, require two different reference 
voltages which are switched alternately by relay KO. 
The function of the M2 - Ml memory chain is to make 
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ELECTRIC ANALOG OF A VENTRICLE 




VENTRICLE 



Vo 
V| 

Rqi 

R,2 
C| 

q 



VENOUS PRESSURE 
; VENTRICULAR PRESSURE 
ARTERIAL PRESSURE 
INLET VALVE RESISTANCE 
OUTLET VALVE RESISTANCE 

VENTRICLE CAPACITANCE = -^, WHERE E = ELASTANCE 
BLOOD VOLUME 

WORK DONE BY THE VENTRICLE DURING SYSTOLE » 
kj^vg i|2 dt, WHERE k = CONSTANT 



FIGURE 2 



the final value of q (blood volume) at the end of the nth 
diastole available throughout the nth systole. (M de- 
notes complementary memory. ) Note that relay KG 
prevents this memory chain from learning the end 
systolic blood volume. 



In order to determine the requiredfunctional relation- 
ships between q^^ and ^^^ , a simple iteration circuit is 



provided. With some ^^ at the output of M3, systolic 

work is computed, and this is compared to the output 
of FG 2 which follows Starling's curve. The error is 
repeatedly used to up-date *jj imtil the desired *jj is 

obtained. The above procedure is repeated for various 
values of Vq until enough information is obtained to set 
up FG 1. With FG 1 in the circuit, a ventricular model 
is ready to be used as a part of a circulation model 
which will not be discussed here. 
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FIGURE 3 
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CHAPTER 22 
A PRACTICAL APPROACH TO ADAPTIVE CONTROL 



Originally published in: CONTROL ENGINEERING, May 1966 
Authors: J. W. Banham Jr. and W. L. Smith 



Controllers are linear elements that are often re- 
quired to operate in nonlinear systems. Thus they 
cannot be expected to provide optimum performance 
over a wide range of system operating conditions. 
However, through a linear representation of the non- 
linear system a controller can be designedwith adap- 
tive features that does provide optimum compensa- 
tion for transient system requirements. The follow- 
ing text describes the design of a controller and its 
conversion to the "adaptive" type (one whose dynamic 
parameters are continuously established as a function 
of transients in the system). The approach described 
has been brpadboarded and tested for a shipboard sys- 
tem but the design techniques are applicable to other 
control systems. 

The maneuverability of a naval ship depends to a large 
degree upon the transient response characteristics of 
the steam generating and machinery systems. Most 
marine steam generators use forced draft; air for 
combustion is supplied to the fuel oil burners by vari- 
able speed blowers. The forced draft blowers in the 
machinery plant of a high speed combat ship are gen- 
erally powered by direct drive steam turbines and 
regulated by controlling the steam flow to the turbines. 
The rotational inertial characteristics of these ma- 
chines, over a wide range of operating speeds, contri- 
bute materially to the "sluggishness" of process re- 
sponse to changes in applied load. Designers interested 
in improving the system response must optimize the 
controller parameters of those system loops or sub- 
loops with the longer time constants. 

SYSTEM DESIGN STUDY 

In a conventional automatic combustion control sys- 
tem, controller settings are usually determined em- 
pirically at the load condition of maximum open loop 
sensitivity. Because controllers are linear elements 
and the forced draft blowers are nonlinear, the situa- 
tion is considerably less than optimum over a wide 
range of operating loads. 

a) Nonlinearities in the blower 

For any type of rotating machinery, the angular ac- 
celeration of the rotor assembly is proportional to the 
applied differential torque and inversely proportional 
to the polar moment of inertia. This linear relation- 
ship may be expressed by the simple differential equa- 
tion: 



where 



= T (t - T ) 



2 
w = angular acceleration, rad/sec 

2 
J = polar moment of inertia, Ib-ft-sec 

T = driving torque developed by the prime 
^ mover, Ib-ft 

T = load torque imposed by the blower, Ib-ft 



Several factors influence the applied torques. The 
torque developed by a steam turbine is a function of 
the turbine steam rate and its rotational speed. A 
typical turbine performance graph illustrating the func- 
tional relationship is shown in Figure 1. The para- 
bolic curve on this map represents the design load re- 
sistance for the turbine- fan unit. For this particular 
set of performance data, it was found practical to rep- 
resent the drivir^ torque by the algebraic expression 



Ta = 



smax, 



max 



G 



smax 



i'^ax 



where 



G /G = turbine steam flow, fraction of 

« ^"^^ rating 



max 



m 



= maximum locked rotor torque, Ib-ft 



a constant, Ib-ft 



cxj/w = angularvelocity, fraction of rating 

' max ^ 

The steam flow to the turbine nozzles varies with the 
steam supply pressure and temperature and the steam 
valve flow coefficient and port area. The forced draft 
blowers considered in this study were equipped with 
a V-ported regulating throttle valve for which the 
characteristic flow- lift curve is as illustrated in Fig- 
ure 2. The supply steam state is reasonably constant 
throughout the load range. The load torque varies di- 
rectly with the blower system resistance and with the 
square of the rotational speed. 



T - 
r 



\a;max/ 
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LINEAR REPRESENTATION OF THE 
NONLINEAR SYSTEM 



All of the elements of the combustion air flow control 
loop of a marine steam generating system are charac- 
terized by some nonlinearity. The following analysis 
illustrates one method of representing the forced draft 
blower nonlinearities by a set of linear functions . 

In the basic acceleration equation 



J ^ a 



T ) 



the torque terms are functions of the turbine steam 
flow, Gg and the turbine-fan speed, oj, as shown in 
thetext of this article. The accelerating torque, T, 

was found to be a function of both steam flow and speed; 
the graphical relationship appears in Figure 1 . In this 
study, the analytical relation between these variables 
was established by curve-fitting to manufacturer's 
performance curves, but can also be developed by a 
detailed derivation of momentum equations applied to 
the turbine wheel. For a fixed blower air system re- 
sistance, the load torque varies with the square of the 
rotational speed. The constant of proportionality, /3 , 
is established by this resistance. 

Steam flow to the turbine is regulated by a throttling 
valve characterized by the arbitrary function 

Gg = F(L) 

where L represents valve lift. 

Sufficient information is now available to permit the 
development of a set of transfer functions which de- 
scribes the open loop frequency response of the tur- 
bine-blower unit to small perturbations about a fixed 
operating level. The general perturbation response 
is given by 



Aw 



-u 



(AT, 



AT^)dt 



If a small disturbance is applied at a particular steady 

state initial condition (L , G , w ), the nonlinear 

o' so' o 



torque terms may be reduced to linear fimctions ex- 
pressed by the following equations : 



AT = AG (T G - vciG^^o)) 
a s max so so o 



AT = 2SCi) Aw 
r o 



Gs is a function of valve lift as shown in Figure 2; a 
tangent constructed at the point (L , G ) has aslope 
given by 



^°3 = ^ nu 



L = L 



AL 



Substitution of these linearized functions into the ac 
celeration equation and applying a Laplace operator 
yields 



Aoj 



Js [dL 



F(L) L ^ L(T G 
^ ' IL = L max so 



- "^Gg^o;^) - 2^co^Aa;j 



For any particular steady state operating level the 
terms d/dLF(L) Il^Lq, T^^G^^, mG^^w^ and 2^ c^ 

are constants ; the above equation may thus be e3q)ressed 
as a transfer fvinction: 



Aoj /„v _ dL 



F(L)[': 



max so 



G__ - mG__a; 1 



so 



1 + 



[j/2^coJl 



This transfer function defines the response of shaft 
speed to valve lift for small perturbations given initial 
conditions defined by L , G , and a; . 



so' 
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0.2 0.4 0.6 0.8 

cj/oj^Q^ Speed, fraction of rating 



I.O 



Figure 1. Forced draft blower turbine performance 
curves. 




0.2 0.4 0.6 0.8 I 

L/L max Valve lift, fraction of full stroke 



Figure 2. Turbine valve steam regulating charac- 
teristics. The arbitrary fimction G = F (L) used in 

Figures 3 and 7 is illustrated. Construction of a tan- 
gent at (G /G = 0. 30, L/L = 0. 53) illustrates 

S S IiidoX 

max 
the graphical approach to the linear representation of 
a nonlinear system. The tangent slope represents a 
valve gain of 1. 12 percent of steam flow per percent 
of valve lift at a point corresponding to 53 percent of 
lift. 



where p is a constant for a particular system resist- 
ance, lb- ft. In the design study, the coefficient p 
was permitted to vary with the effective fan discharge 
area (established by the number of burner air regis- 
ters in use in the steam generator). 

This information led to the model of the system shown 
in Figure 3. The problem at this point was to devise 
a control law which would provide optimum compensa- 
tion for the nonlinear system over a wide rai^e of 
operatii^ speeds. For ttds design, the system is in- 
itially described in the frequency domain as a contin- 
uous series of linear systems, each similar in form 
to the others but having different parameters. 

b) Linear representation of the nonlinear system 

One approach to the linear representation of a non- 
linear system is to expand the nonlinear elements of 
the system by means of the Taylor series, then dis- 
regard all but the linear terms. This produces a 
linear system which approximates the real system re- 
sponse to small perturbations about a fixed operating 
point. Because the steam admission valve character- 
istic is an essentially arbitrary fimction, a graphical 
approach was used. In this method, all of the non- 
linear terms are represented as linear functions by 
constructir^ tangents at selected operating points along 
the curves which represent the functions of the alge- 
braic variables. The slope of each tangent curve is 
taken as the gain of the element under consideration, 
Figure 2. The resulting linear model is illustrated in 
Figure 4. The development of the values of the pa- 
rameters associated with this model is detailed in the 
box on page 22-2. n 

A linear model like the one in Figure 4 leads to a set 
of transfer functions instead of a single fimction. 
Each function in the set represents the open loop fre- 
quency response of the system to a small signal si- 
nusoidal disturbance superimposed on a particular 
steady state initial condition. The validity of the 
analytical model constructed in this fashion was con- 
firmed by comparison with actual frequency response 
data obtained from the operating units subjected to 
pulse test analysis at two load conditions, (Ref. 1). 
The family of transfer functions thus obtained was 
used as the basis for a series of control system de- 
si^s t)ased on classical frequency response tech- 
niques applicable to linear systems. 

DESIGN OF THE CONTROL SYSTEM 

The following conditions were used as criteria for the 
dynamic performance of the closed- loop system: 

a. There must be no offset error in the steady 
state. 

b. No overshoot in the transient response to a 
"ramp" input disturbance is to be permitted. 

c. The time integral of the system absolute error 
is to be minimized. 

The conventional control loop employs a proportional- 
plus- integral controller which receives as its feed- 
back the output signal from a combustion air flow 
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DEVELOPMENT OF CONTROLLER DESIGN 



Conventional controllers used for regulating forced 
draft blowers are simple proportional-plus- integral 
devices characterized by fixed gain and reset rate 
over the entire operating range. The output signal 
of a linear controller of this type is given by 



where — = variable controller reset rate, rad/sec 



— = optimum controller reset rate 
o 



9 = ye 



-u 



edt 



where e 

r 
1 

T 



error between input and feedback signals 
controller proportional gain 

— = controller reset rate 



max 



atw/oj ^ 

max 



= optimum controller reset rate 

^* '^/'^max = ^ 



An adaptive feature was incorporated in such a con- 
troller by servoing the proportional gain and reset 
rate (integral crossover frequency) settings as re- 
quired to maintain aa optimum combination of closed 
loop system response and stability over the designed 
range of forced draft blower operation. These con- 
troller parameters are shown plotted against various 
loop operating levels in Figure 6. Curve fitting of 
these graphical functions produced the following alge- 
braic expressions: 



The output of the ad^tive air flow controller may 
therefore be expressed by the equation 



^ max / o J 



T „^ (oj/w ) / 

max \ max/y 



y/ = 



OJ/O) 



max 



where y' = variable controller gain 



7 = optimum controller gain at o) = oj 



max 



It is to be noted that the controller parameters have 
been e3q)ressed as fimctions of the control feedback 
signal (in this particular system, the blower rota- 
tional speed). To produce maximum loop gain during 
any transient disturbance, the controller gain and 
reset rate are automatically varied as a fimction of 
either the command signal * or the feedback signal 
(^, whichever of the two is least. An arrangement which 
produces this effect is illustrated in the computer sche- 
matic. Figure 7. Using the notation indicated on this 
diagram, the controller equation becomes 



w/oj = steady state loop operating level 



o max \ max/ 



6' = y 



WPmax) ^ "" ^o y 



edt 



max 



u^vax/y^"** 
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Figure 3. Forced draft blower block diagram. 
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Figure 4. Block diagram of forced draft blower perturbation model at (Lq , c^o^' 



transmitter. For the "adaptive" design, it was elected 
to utilize the same basic controller concept but to vary 
the proportional gain and integral cross- over frequency 
to maintain optimum performance. 

The conventional design was initiated by selecting con- 
troller parameters which would give the entire control 
loop a gain margin of six db and a phase margin of 35 
deg; these margin goals were derived from previous 
experience. This conventional procedure was repeated 
for each of the forced draft blower transfer functions 
at the various operating load levels. The result, de- 
tailed in the box on page 22-4, wasasetof empirically 
optimum controller settings which could be correlated 



with the loop input or output to develop analytic func- 
tions that would specify the adaptive nature of the 
design. 

The air flow controller thus obtained was simulated 
together with the forced draft blowers and air flow 
system on an analog computer. The transient response 
to input disturbances comparable to the effects of ship 
maneuvering were then studied. Typical results are 
shown in Figure 5. Note that a small amount of over- 
shoot occurred which violates one of the system de- 
sign criteria. This observation led to a repetition 
of the design procedure, this time using gain and phase 
margins of 9 db and 45 deg, respectively. Figure 6 
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illustrates the relationship between the controller 
parameters and the loop operating level for perform- 
ance with the revised open- loop margins. 

A revised analog computer program was then pre- 
pared; a simplified version of the basic program is 
shown in Figure 7. Transient response of this simu- 
lated system appears in the chart recordings of Fig- 
ure 8A; the corresponding response of the simulated 
adaptive control system appears for comparative pur- 
poses in Figure 8B. It is evident from these record- 
ings that the adaptive system is considerably superior 
to the conventional system. Inspection of the error 
recordings discloses a reduction of five to one in the 
time integral of absolute error. At the same time 
there is neither overshoot nor permanent offset. Thus 
the design criteria have been met. 



FROM SIMULATION OF ACTUAL TESTS 



Figure 5. System transient response with initial 
optimization criteria. 
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Figure 6. Controller parameters (servoed controller 
settings) as a fimction of loop operating level in the 
steady state. Open loop gain and phase marginsare 
9 db and 45 deg. 



Results obtained from the analog computer studies 
supported a decision to breadboard a pneumatic con- 
trol system from standard hardware. The basic com- 
ponents used in the adaptive controller were pneu- 
matic proportional and integral amplifiers with servo- 
driven gain settings. Two of these units provided the 
necessary proportional-plus- reset control mode. The 
gains of these two devices were automatically varied 
as a function of either the input to the closed air flow 
control loop or its output, whichever produced the 
greatest loop gain chiring any transient disturbance. 
In the steacfy state, these variables are equal and 
produce the same gain in the open loop. 



Tests of the loop were conducted by imposing "ramp" 
inputs which corresponded to those used in the com- 
puter simulation, thus aiding in an evaluation of con- 
trol system performance. The recorded results of 
one of these transients are shown in Figure 8C. The 
transient performance of the closed loop system 
closely paralleled that of the simulated system, see 
Figure SB. 



The air flow control loop comprises one of the open- 
loop elements of the steam generator's automatic 
conibustion control system. The principal benefit of 
the adaptive air flow control system when used as a 
part of the overall steam generator control system 
is its ability to reduce the deviation of controlled 
steam pressure and to eliminate the possibility of 
incomplete combustion (and accompanyii^ smoke) 
during maneuvering. Although hardware requirements 
for systems of this type are presently more com- 
plex, and less reliable than for conventional systems, 
adaptive design techniques become attractive when 
applied to systems constructed from encapsulated or 
modular design control circuitry. The U. S. Navy is 
currently investigating adaptive systems of the solid 
state pneumatic (pure fluid) and solid state electronic 
types for use in the control of shipboard main pro- 
pulsion and auxiliary machinery. 
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Figure 7. Simplified analog computer program used to study conventional 
and adaptive controller responses . 
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Figure 8. Transient response recordings. A- simu- 
lated conventional controller. B- simulated adaptive 
controller. C-breadboard adaptive controller. 
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APPENDIX I 



Glossary of Abbreviations 



c 


Capacitor 


c 


Compute 


c 


Clock 


CK 


Check 


D or DN Down 


DCU 


Decimal Counting Unit 


DFG 


Diode Function Generator 


DVM 


Digital Volt Meter 


ES 


Electronic Switch 


Fb 


Feedback 


FF 


Flip-flop 


F-RT 


Forward-Reset Real Time 


FR 


Function (Operational) Relay 


FS 


Function Switch 


H 


Hold 


1 or IC 


Initial Condition 


IJ 


Initial Junction 


IN 


Input 


INV 


Inverter 



10 Iterative Operation 

J Junction 

K Relay 

L Logical Variable 

M Meter 

NC Normally Closed 

NO Normally Open 

Output 

P or POT Coef. Potentiometer 

R Reset 

R Resistor 

R-H Reverse-Hold 

RO Repetitive Operation 

RTS Reset-Toggle-Set 

S Set 

SJ Summing Junction 

T Toggle 

TL Trunk Line 

U Up 



APPENDIX II 



The symbology used throughout this handbook is in conformity with the objectives 
of Simulation Councils, Inc., namely, to encourage the use of uniform analog com- 
puter graphics. This article was prepared by SCI, and appeared in the March 1966 
issue of SIMULATION. Permission has been granted to reprint this article. 



Presenting... 

Uniform graphics for SIMULATION 



The symbols and the methods of laying out analog 
and hybrid computer diagrams presented here are 
advocated to alleviate the confusion caused by the 
uncoordinated invention of new symbols and dia- 
gramming practices. The increasing use of hybrid 
techniques and equipment has aggravated an already 
bad situation to the point that it is often no longer 
possible for one worker in our field to read another's 
diagram. Usually this is because symbols are devised 
and diagrams are drawn to include details peculiar 
to a particular kind of equipment. Such a wiring, or 
"patching," diagram is of course necessary for set- 
ting up and checking out an actual simulation, but 
hardware-peculiar details are only confusing to those 
with other kinds of equipment. With few exceptions, 
the use of a simplified signal-flow diagram to illus- 
trate technical articles is much more effective. 

With the foregoing in mind an SCI committee 
composed of 

GEORGE BURGIN 

JOE HUSSEY 

HANS JORGENSEN 

GRANINO KORN 

JOHN McLEOD 

selected the symbols and offers the following sugges- 
tions for their use. Primary considerations were cur- 
rent usage, clarity, and simplicity. We devised no new 
symbols and, unless there were overriding indica- 
tions to the contrary, we adopted those already in 
widest use. Clarity and simplicity, we believe, will be 
enhanced by the choice of unique shapes to repre- 
sent different components, and the elimination of all 
unnecessary details in diagrams. 

There was no intent on the part of our committee 
to set up standards for the industry. However, all 
diagrams appearing in SIMULATION will be pre- 
pared according to the committee's recommenda- 
tions (as they may be modified from time to time), 
and we hope that these recommendations will prove 
attractive to others. Suggestions for modifications 
and additions are solicited. 



General rules 

The following methods and symbols are recom- 
mended for the illustration of technical articles pre- 
pared for publication. Unless the purpose of the 
article is to describe the use of a particular kind of 
equipment, and the hardware details are pertinent 
to the subject, such illustrations should not be "hard- 
ware-peculiar." In other words, the objective should 
be to show signal flow, rather than "patching" 
details. 

The primary, or overall, system diagram should 
show only the essential signal flow. Where it is neces- 
sary to show details, separate diagrams should be 
made and referenced to the primary diagram by en- 
closing the detailed area of the primary diagram in 
dotted lines with suitable notation. 

The direction of signal flow should be indicated 
by arrowheads except where the shape of the sym- 
bols makes the direction of flow obvious. Primary 
signal flow (with the exception of feedback loops) 
should be from left to right, and, if practical, each 
"line" of symbols should be made to read like the 
mathematical relation it represents. 

The choice of whether to end a line and label it 
(preferably with the symbol of the variable that the 
signal represents) when it reaches the right-hand 
side of the diagram, and then indicate its continua- 
tion with the same label as it enters again at the 
left-hand side, instead of drawing in the connection, 
should be made on the basis of clarity; if a line re- 
turning the signal from right to left will cross many 
other lines and be hard to follow, it should not be 
drawn in. 

if a diagram involves a number of identical cir- 
cuits, only one should be shown in detail, while the 
others should be indicated by boxes with appropriate 
notation. 

Components should not be numbered ynless they 
are referred to by number in the text. 

All amplifier gains should be shown just outside 
the amplifier at the point where the input enters. 
Unity gains should not be labeled. 



Always apply the test of clarity and simplicity. Ask 
yourself: "Is this the most understandable way to 
diagram this for those unfamiliar with the hardware, 
and less familiar with the subject, than I?" 
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ELEMENT 



ISignai flow 



FUNCTION 



Analog Binary 

*- — 



Digital word Mechanical 

— «^ 




NOTES 



Distinction as to kind of signal is made only at 
beginning and end of line, and then only where 
there are two or more kinds of signals to be shown 
on one diagram. Arrowheads should be omitted 
only if direction of signal flow is obvious. 



Connections 



No connection 



! Integrators 



Amplifier 
a 


10 


b 


5 


1 




1 




x=— (10a + 56 + c) 



Summing amplifier (with high-gain input) 
a 10 



High-gain amplifier 
a 10 



Inverting amplifier 



Standard 




x=— (10a + 5b-f,iC) 




> 



<=— ^(10a + 51b + c) 




-/(10b + 5c-f cOdt 



Two-mode 




RESET when R=1 
COMPUTE when R = 



No arrows; shape clearly indicates direct^ 

flow. 

No label on unity-gain input. 

"N" indicates where number of amplifier should 
be placed if, and only if, referred to in text. J 



Label high-gain input "SJ" or "G" only if perti- 
nent; the point at which c enters (at top or bot 
tom of symbol) indicates a gain of p 






d 

lowoniy 



Note that the numbers at the inputs are nowonli 
relative gains; /t, an unspecified large number 
which is a function of the particular amplifier, is 
the overall gain. 



When amplifiers are used only to change signal 
sign they should be smaller. Orientation should 
be that which makes for greatest graphical sim- 
plicity, i.e., they may be shown in vertical, or in 
feedback, paths. 

Standard integrator mode is that of basic problem. 
The point at which a enters (at top or bottom of 
symbol) indicates that it sets the initial value of x, 
so it is redundant to label it IC. 




Arrows indicating mode control should be shown 
only when two or more integrators are operating 
in different modes simultaneously. The mode- 
control input should then be labeled or coded to 
indicate the controlling mode. ~ 



Three-mode 




S R 



Track-hold 




RESET when R===1 

or R = S = 1 
HOLD when S=1, R = 
COMPUTE when R=S = 



TRACK when R = 1 
HOLD when R=0 



When used in this way the size of the integrating 
capacitor may be important, as it affects the track- 
ing lag. In such cases the value should be indi- 
cated as shown. The quantity tracked or held is 

-(a + b + c). 
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ELEMENT 



Integrators 
(continued) 




an generators Arbitary 



k 



FUNCTION 



Integrator 1 in 
Track for 5=1 

Integrator 2 in 
Track for R = 1 



(load connected) 



a = 1,fa = 
{load connected) 



x = ka 




x=/c*a 




NOTES 



Integrator 1 in ^TRl»ieK for S=«rrl; integi'at^' 2 is' 
TRACK for R = 1. In certain iterative problems y 
is required to have some value for the first itera- 
tion; thus an "initial" initial condition must be 
furnished as shown. 




The mathematical relationship shown for the 
three-terminal pot is given only to define k. It 
will not usually hold during problem solution be- 
cause if b were always zero a two-terminal pot 
could be used. If it is desirable to define k in any 
other way, its meaning should be clearly indicated. 

This should be used only if the fact that N is servo- 
set is significant to the solution of the problem, as, 
for instance, when it must be automatically reset 
after each run of an iterative problem. 



x = /(a) 



As shown, the symbol indicates a digital pot with 
integral (committed) inverting amplifier. The sym- 
bol for non-inverting digital pots should not in- 
clude the amplifier. M 

*lndicates quantized value. a 



The number N should be used only when the 
function generator is referred to by number in 
the text, or when it is desirable to explain what 
the function f is, either by a footnote, or by a 
separate graph. 



Mathematical functions (typical) 



> 



Multiplier 



+a 



X~ 



Divider 



x = va 



-X, x=- 



afe 



a = reference voltage 



In cases where the function generated can be rep-~ 
resented by a standard mathematical symbol, the 
/ should be replaced by the symbol 




These symbols should be used for all analog! 
pliers and dividers; if the kind is significant to the 
solution of the problem, it should be so stated 
in the text and/or noted on the diagram. Because 
equipment differs, the sign of the output for spe- 
cified signs of both inputs should be given, other- 
wise inversion will be assumed. 

The number N should be used only if the com- 
ponent is referenced in the text, or if more than 
one channel of a multichannel device is used. In 
the latter case the subscript, in this case /, identi- 
fies the channel. In many cases, it will be found 
that a drawing can be simplified and clarified by 
drawing the same multichannel device in more 
than one signal flow path. In this case the number 
would be the sam.e, but the subscript would be 
different for each channel. J 



Resolvers x= R cos 

Polar to rectangular Rectangular to polar y=^Rs\n6 



£ 9==arctany/x 




R=:VxM-y5 



If a resolver has additional outputs they should be 
shown only if used. 




Comparators 
With binary output With true and false (^ = 1 (a < fa) 

binary outputs Ut=0 (a > fa) 



N 



U 




N 



>* 



The symbol for a relay comparator can, of course, 
be made by combining either of these symbols 
with that of a relay. 
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Relay 





NOTES 



If the diagram is drawn as shown, the designations 
of arms and contacts are superfluous. However, 
it is often desirable, for clarity, to show the relay 
coil in one place on a diagram and the contacts 
in signal-flow paths elsewhere. In this case, N„ 
would designate the normally Open (L' = 0) con- 
tact of channel a, relay N; N^ the Arm of channel a 
of relay N, etc. 



Digital inverter 



Digital/analog switch 






The tilde (— ') is used instead of a bar to indicate 
negation. This is recommended because of the 
many other meanings of a bar over a variable. 

ON when binary signal, U, is digital "one." 



U^==a (.U=V 



Limiters 



Feed-back 



Converters 




i Diode gates 




3r ^ 




N 


>< 




+u 








— / 











Analog-to-digital 



Digital-to-analog 



5 — ^ 


y. 


X* 


Solid state 

^ 1 


1 


AND 
a 


-TN. 


h 


^ ^ 


c 


. J 




Xy 


OR 
a 




b 




c 





H 




A 


M 


NAND 
a 
b 
c 


^ 


NOR 
a 
b 
c 


i — '- 



x=— a (/<a<u) 
x=— / (a < /) 

x= — u (a > u) 



x = a (/<a<(j) 
x==:/ {a</) 
x = i; (a > u) 



x'=a 
x = a* 



x=:a 
x = 


(a>x) 
(a<x) 


AND 


x = a-6-c 


NAND 


\:^3i-b-C 


OR 


x^a^b^c 


NOR 


X =/= a -f- fa -f- c 



u= larger value at which output is limfl 
/ = smaller value at which output is limited 



•Indicates quantized value. 



Output is high (logic 1) if and only if all inpud 
are high. ;J 

Output is low (logic 0) if and only if all inputs* 
are high. 

Output is high (logic 1) if and only if at least one 
input is high [i.e., any or all of the inputs are high). 

Output is low (logic 0) if and only if at least one 
input is high (i.e., any or all of the inputs are high). 



Typical 



n^ 



VV 



As there are many kinds of flip-flops, symbolic 
representation of the operation of any but the 
simplest is not advised; if the operation is not 
obvious, an explanation in the text or a footnote 
with the diagram is recommended. 



TT 



This is Black Box number N, to be used if no sym- 
bol is given here, or in case of doubtl It can have 
any number of inputs and outputs, but the nature 
of the signals should be indicated and all signal- 
flow directions should be designated by arrow- 
heads. Here there are three analog inputs giving 
rise to one analog and one binary output. 



The function of a "black box" can often be 
made obvious by properly labeling the inputs 
and outputs. However, if it is not obvious, the 
function should be explained in the text or by a 
footnote with the diagram. If internal details are 
of interest, they should be shown in an appro- 
priately labeled auxiliary diagram. 
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